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Abstract 

The Kuramoto model is a system of ordinary differential equations for describing 
synchronization phenomena defined as a coupled phase oscillators. In this paper, a bifur- 
cation structure of the infinite dimensional Kuramoto model is investigated. For a certain 
non-selfadjoint linear operator, which defines a linear part of the Kuramoto model, the 
spectral theory on a space of generalized functions is developed with the aid of a rigged 
Hilbert space to avoid a continuous spectrum on the imaginary axis. Although the linear 
operator has an unbounded continuous spectrum on a Hilbert space, it is shown that it 
admits a spectral decomposition consisting of a countable number of eigenfunctions on 
a space of generalized functions. The semigroup generated by the linear operator is cal- 
culated by using the spectral decomposition to prove the linear stability of a steady state 
of the system. The center manifold theory is also developed on a space of generalized 
functions. It is proved that there exists a finite dimensional center manifold on a space 
of generalized functions, while a center manifold on a Hilbert space is of infinite dimen- 
sional because of the continuous spectrum on the imaginary axis. The results are applied 
to the stability and bifurcation theory of the Kuramoto model to obtain a bifurcation dia- 
gram conjectured by Kuramoto. If the coupling strength K between oscillators is smaller 
than some threshold K c , the de-synchronous state proves to be asymptotically stable, and 
if K exceeds K c , a nontrivial solution, which corresponds to the synchronization, bifur- 
cates from the de-synchronous state. 
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Fig. 1 : The order parameter of the Kuramoto model. 



1 Introduction 

Collective synchronization phenomena are observed in a variety of areas such as chemical 
reactions, engineering circuits and biological populations 071 . In order to investigate such 
phenomena, Kuramoto B5Tl proposed the system of ordinary differential equations 



where 0, = 0,-(f) e [0, 2n) is a dependent variable which denotes the phase of an z-th 
oscillator on a circle, a>i e R denotes its natural frequency, K > is a coupling strength, 
and where N is the number of oscillators. Eq. dl.ll) is derived by means of the averaging 
method from coupled dynamical systems having limit cycles, and now it is called the 
Kuramoto model. 

It is obvious that when K = 0, #,-(?) and 6j(t) rotate on a circle at different velocities 
unless (x>i is equal to coj, and it is true for sufficiently small K > 0. On the other hand, if K 
is sufficiently large, it is numerically observed that some of oscillators or all of them tend 
to rotate at the same velocity on average, which is called the synchronization 071 l42ll . 
If N is small, such a transition from de-synchronization to synchronization may be well 
revealed by means of the bifurcation theory IfTTl 1271 l28ll . However, if N is large, it is 
difficult to investigate the transition from the view point of the bifurcation theory and it is 
still far from understood. 

In order to evaluate whether synchronization occurs or not, Kuramoto introduced the 
order parameter r{t)e^^ by 



where r, ij/ e R. The order parameter gives the centroid of oscillators. It seems that if 
synchronous state is formed, r(t) takes a positive number, while if de- synchronization is 
stable, r(t) is zero on time average (see FigfJJ). Further, this is true for every t when N 
is sufficiently large so that statistical-mechanical description is applied. Based on this 
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Fig. 2: Typical bifurcation diagrams of the order parameter for the cases that (a) g(a>) is 
even and unimodal (b) g(of) is even and bimodal. Solid lines denote stable solutions and 
dotted lines denote unstable solutions. 



observation and some formal calculations, Kuramoto conjectured a bifurcation diagram 
of r(t) as follows: 

Kuramoto's conjecture 

Suppose that N — » oo and natural frequencies a^ 's are distributed according to a proba- 
bility density function g(uS). If g(co) is an even and unimodal function such that g"(0) ^ 0, 
then the bifurcation diagram of r(t) is given as Figf2](a); that is, if the coupling strength K 
is smaller than K c := 2/(ng(0)), then r(t) = is asymptotically stable. On the other hand, 
if K is larger than K c , the synchronous state emerges; there exists a positive constant r c 
such that r(t) = r c is asymptotically stable. Near the transition point K c , r c is of order 



A function g(to) is called unimodal (at a> = 0) if g(co\) > g(u> 2 ) for < a>i < oj 2 and 
g(cc»i) < g(a>2) for a>\ < u> 2 < 0. Now the value K c = 2/(ng(0)) is called the Kuramoto's 
transition point. See ll26l and 11421 for the Kuramoto's discussion. 

In the last two decades, many studies to confirm the Kuramoto's conjecture have been 
done. Significant papers of Strogatz et al. 11431 1441 investigated the linear stability of 
the trivial solution, which corresponds to the de-synchronous state r = 0. In [|43l . they 
introduced the continuous model for the Kuramoto model to describe the situation iV — » 
oo. They derived the Kuramoto's transition point K c = 2/(ng(0)) and showed that if K > 
K c , the de-synchronous state is unstable because of eigenvalues on the right half plane. 
On the other hand, when < K < K c , the linear operator T u which defines the linearized 
equation of the continuous model around the de-synchronous state, has no eigenvalues. 
However, in IPBI . they found that an analytic continuation of the resolvent (A - T\)~ l 
may have poles {resonance poles) on the left half plane, and with the aid of numerical 
simulations, they remarked a possibility that resonance poles induce exponential decay 
of the order parameter r. This claim will be rigorously proved in this paper. In [33], the 
spectra of linearized systems around other steady states, which correspond to solutions 
with positive r = r c , are investigated. They found that linear operators, which define 
linearized systems around synchronous states have continuous spectra on the imaginary 
axis. Nevertheless, they again remarked that such solutions can be asymptotically stable 
because of the resonance poles. 



0((K - K c ) 1 ' 2 ). 
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Since results of Strogatz et al. are based on the linearized analysis, effects of nonlin- 
ear terms are neglected. To investigate nonlinear dynamics, the bifurcation theory is often 
used. However, investigating the bifurcation structure near the transition point K c involves 
further difficult problems because the operator T\ has a continuous spectrum on the imag- 
inary axis, that is, a center manifold in a usual sense is of infinite dimensional. To avoid 
this difficulty, Bonilla et al. EH El and Crawford et al. O 021 EH) added a perturbation 
(noise) with the strength D > to the Kuramoto model. Then, the continuous spectrum 
moves to the left side by D, and thus the usual center manifold reduction is applicable. 
When g(to) is an even and unimodal function, they obtained the Kuramoto's bifurcation 
diagram (Fig.2 (a)), however, obviously their methods are not valid when D = 0. For 
example, in Crawford's method, an eigenfunction of T\ associated with a center subspace 
diverges as D — » because an eigenvalue on the imaginary axis is embedded in the con- 
tinuous spectrum as D — » 0. Thus the original Kuramoto's conjecture is still open. 

Despite the active interest in the case that the distribution function g(uS) is even and 
unimodal, bifurcation diagrams of r for g(oj) other than the even and unimodal case are 
not understood well. Martens et al. Il30l investigated the bifurcation diagram for a bimodal 
g(oj) which consists of two Lorentzian distributions. In particular, they found that stable 
synchronous states can coexist with stable de-synchronous states if K is slightly smaller 
than K c (see Figf2](b)). Their analysis depends on extensive symmetries of the Kuramoto 
model found by Ott and Antonsen ll35l[36ll (see also BTlO and on the special form of g{oS), 
however, such a diagram seems to be common for any bimodal distributions. 

In this paper, the stability, spectral and bifurcation theory of the continuous model of 
the Kuramoto model will be developed to prove the Kuramoto's conjecture. In particular, 
the spectrum and the semigroup of a linear operator Ti on the weighted Lebesgue space 
L 2 (R, g(oj) dco), which is obtained by linearizing the continuous model around the de- 
synchronous state, will be investigated in detail. At first, we derive the transition point 
(bifurcation point) K c for any distribution function g{oS). As a corollary, the Kuramoto's 
transition point K c = 2/(ng(0)) is obtained if g(co) is an even and unimodal function. 
When K > K c , it is proved that the de-synchronous state is unstable because the operator 
Ti has eigenvalues on the right half plane. 

On the other hand, when < K < K c , the operator T\ has no eigenvalues and the 
continuous spectrum lies on the imaginary axis. Thus the stability of the de-synchronous 
state is nontrivial. Despite this fact, under appropriate assumptions for g(to), the order 
parameter proves to decay exponentially to zero as t — » oo because of the existence of 
resonance poles on the left half plane, as was expected by Strogatz et al. MM- To prove 
it, the notion of spectrum is extended. Roughly speaking, the spectrum is the set of 
singularities of the resolvent (A - T\Y l . However, if g{aj) has an analytic continuation, the 
resolvent proves to have an analytic continuation, and it has singularities on the second 
Riemann sheet, which are called the resonance poles. By using the Laplace inversion 
formula for a semigroup, we will prove that the resonance poles induce exponential decay 
of the order parameter. It suggests that in general, linear stability of a trivial solution of 
a linear equation on an infinite dimensional space is determined by not only the spectrum 
of the linear operator but also its resonance poles. 

Next purpose is to investigate a bifurcation from K = K c . To handle the continuous 
spectrum on the imaginary axis, a spectral theory of the resonance poles is developed with 
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the aid of a rigged Hilbert space (Gelfand triplet). A rigged Hilbert space consists of three 
spaces X c H c X': a space X of test functions, a Hilbert space H (in our problem, this 
is the weighted Lebesgue space L 2 (R, g(a>)doj)) and the dual space X' of X (a space of 
continuous linear functionals on X called generalized functions). A suitable choice of X 
depends on g{oS). In this paper, two cases are considered: (i) g(to) is analytic and decays 
faster than any exponential functions (e.g. Gaussian distribution), (ii) g(co) is a rational 
function (e.g. Lorentzian distribution g(co) = 1/(tt(1 +oj 2 ))). For the case (i), X := Exp + is 
a space of holomorphic functions <f>(z) defined near the real axis and the upper half plane 
such that sup Im( . )> _ £ \(p(z)\e~^ is finite for some s > and /3 > 0. For the case (ii), X := H + 
is a space of bounded holomorphic functions on the real axis and the upper half plane. For 
both cases, it is proved that the operator T\ admits the spectral decomposition on X' by 
regarding the resonance poles as eigenvalues. Corresponding eigenvectors in X' are called 
the generalized eigenfunctions associated with resonance poles. Despite the fact that T\ is 
not a selfadjoint operator and it has the continuous spectrum, the spectral decomposition 
of T\ consists of a countable number of generalized eigenfunctions: roughly speaking, 
any element (f> in X is decomposed as (f> = 2," PnW) " A?> where fi n e X' is a generalized 
eigenfunction. Further, it is shown that for the case (ii), the decomposition is reduced to a 
finite sum because of a certain degeneracy of the space X = H + . We further investigate the 
semigroup generated by T\ and the projection to the eigenspace in the space of generalized 
functions X' . It is proved that the semigroup e Tl ' behaves as e Ti 'cp = Y,7=o ^"'Unity) • for 
any (p e X, where /t„'s are resonance poles. This completely determines the dynamics of 
the linearized Kuramoto model. In particular, when < K < K c , all resonance poles lie 
on the left half plane: Re(/l„) < 0, which proves the linear stability of the de-synchronous 
state. When K = K c , there are resonance poles on the imaginary axis. We define a 
generalized center subspace E c on X' to be a space spanned by generalized eigenfunctions 
associated with resonance poles on the imaginary axis. It is remarkable that though the 
center subspace in a usual sense is of infinite dimensional because of the continuous 
spectrum on the imaginary axis, the dimension of the generalized center subspace on X' 
is finite in general. The projection operator to the generalized center subspace will be 
investigated in detail. 

Note that the spectral decomposition based on a rigged Hilbert space was originally 
proposed by Gelfand et al. lfT8ll29ll . They proposed a spectral decomposition of a self- 
adjoint operator by using a system of generalized eigenfunctions, however, it involves 
an integral; that is, eigenfunctions are uncountable. Our results are quite different from 
Gelfand's one in that our operator T\ is not selfadjoint and its spectral decomposition 
consists of a countable eigenfunctions. 

Finally, we apply the center manifold reduction to the infinite dimensional Kuramoto 
model by extending it to an evolution equation on X' . Since the generalized center sub- 
space is of finite dimensional, a corresponding center manifold on X' seems to be a finite 
dimensional manifold. However, there are no existence theorems of center manifolds on 
X' because X' is not a Banach space. To prove the existence of a center manifold, we in- 
troduce a topology on X in a technical way so that the dual space X' becomes a complete 
metric space. With this topology, X' becomes a topological vector space called Montel 
space, which is obtained as a projective limit of Banach spaces. This topology proves to 
have a very convenient property that every weakly convergent series in X' is also conver- 



6 



gent with respect the metric. By using this topology and the spectral decomposition, the 
existence of a finite dimensional center manifold for the Kuramoto model will be proved. 
The dynamics on the center manifold will be derived when g(co) is an even and unimodal 
function. In this case, the center manifold on X' is of one dimensional, and we can show 
that the synchronous solution (a solution such that r > 0) emerges through the pitchfork 
bifurcation, which proves the Kuramoto's conjecture. A bifurcation structure depends 
on the shape of g(to). When g(to) is an even and bimodal function, it is shown that the 
synchronous solution appears through a subcritical bifurcation as is shown in Figj2](b). 

This paper is organized as follows: In Sec. 2, the continuous model for the Kuramoto 
model is introduced and its basic properties are reviewed. In Sec. 3, Kuramoto's transition 
point K c is derived and it is proved that if K > K c , the de-synchronous state is unsta- 
ble because of eigenvalues on the right half plane. In Sec.4, the linear stability of the 
de-synchronous state is investigated. We will show that when < K < K c , the order 
parameter decays exponentially to zero as t — » oo because of the existence of resonance 
poles. In Sec. 5, the spectral theory of resonance poles on a rigged Hilbert space is de- 
veloped. We investigate properties of the operator T u the semigroup, eigenfunctions, 
projections by means of the rigged Hilbert space. In Sec. 6, the nonlinear stability of the 
de-synchronous state is proved as an application of the spectral decomposition on the 
rigged Hilbert space. It is shown that when < K < K c , the order parameter tends to 
zero as t — » oo without neglecting the nonlinear term. The center manifold theory will 
be developed in Sec. 7. Sec. 7.1 to Sec. 7.4 are devoted to the proof of the existence of a 
center manifold on the dual space X' . In Sec.7.5, the dynamics on the center manifold is 
derived, and the Kuramoto's conjecture is solved. 



2 Continuous model 

In this section, we introduce a continuous model of the Kuramoto model and show a few 
properties of it. 

For the Af-dimensional Kuramoto model (11.11) . taking the continuous limit N — » oo, 
we obtain the continuous model of the Kuramoto model, which is an evolution equation 
of a probability measure p, = p,(6, to) on S 1 = [0, In) parameterized by t e R and weR, 
defined as 

77(0:= f f e^ Pt (6,(o)g(to)d6da>, (2-1) 
JrJo 

_ Po (G,(o) = h(e), 

where h(8) is an initial condition and g(to) is a given probability density function for 
natural frequencies. We are assuming that the initial condition h(6) is independent of to. 
This assumption corresponds to the assumption for the discrete model (11.11) that initial 
values {8j(0)} N j =l and natural frequencies {cOj} N . =l are independently distributed, and is a 
physically natural assumption often used in literature. However, we will also consider to- 
dependent initial conditions h(9, to), a probability measure on S 1 parameterized by to, for 
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mathematical reasons, in Sec.7. Roughly speaking, p t (8, co) denotes a probability that an 
oscillator having a natural frequency co is placed at a position 6 (for example, see 01 [Hfl 
for how to derive Eq. d2.ll) ). Since h and p t are measures on S \ they should be denoted 
as dh(9) and dp,(9, ■ ), however, we use the present notation for simplicity. The 77(f) is 



a continuous version of (11.21) . and we also call it the order parameter. 77(f) denotes the 



complex conjugate of 77(f). We can prove that Eq. (l2.1l) is a proper continuous model in 



the sense that the order parameter (11.21) of the TV-dimensional Kuramoto model converges 
to 77(f) as N — » 00 under some assumptions, see Chiba lUOll . The purpose in this paper is 
to investigate the dynamics of Eq. (l2.1l) . 

A few properties of Eq. (l2.1l) are in order. It is easy to prove the low of conservation 
of mass: 

I I p t (e,oS)g(io)d6d(o= I I h(8)g(co)d9dco = 1. (2.2) 
JrJo JrJo 



By using the characteristic curve method, Eq. (l2.1l) is formally integrated as follows: Con- 
sider the equation 

^ = co + ^—irjiDe-^-lme^), xe[0,27r), (2.3) 
dt 2^T\ 

which defines a characteristic curve. Let x = x(t, s; 9, cS) be a solution of Eq. (l2.3l) satisfy- 
ing the initial condition x(s, s; 6, to) = 9 at an initial time s. Then, along the characteristic 



curve, Eq. fl2.ll) is integrated to yield 



p t (8, co) = h(x(0, f; 8, co)) exp[| JW)*-" + ^e^^ds}. (2.4) 



By using Eq. (l2.4l) . it is easy to show the equality 



pzn pin 

I a(9,co)p t (9,co)d9= I a(x(t,0;9,co),co)h(9)d9, (2.5) 
Jo Jo 



for any measurable function a(9, co). In particular, the order parameter 77(f) are rewritten 
as 

rj(t)= f f e' rix( '' ' eM) g(co)h(9)d9dco. (2.6) 
JrJo 

Substituting it into Eqs. fl2.3l) and (12.41) . we obtain 



d n n2n 



—x(t, s;9,co) = co + 
dt 



K f f sm(x(t,0-9',co')-x(t,s;9,co))g(co')h(9')d9'dco', (2.7) 
JrJo 



and 



p t (9, co) = h(x(0, f; 9, co)) x 

<Ln 



Xt p pin 

ds- J J cos(x(s,0-8\co')-x(s,t-9,co))h(9')g(co')d9'dco'],(2.S) 
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respectively. They define a system of integro-ordinary differential equations which is 
equivalent to Eq. (12. II) . Even if h(6) is not differentiable, we consider Eq. (l2.8l) to be a weak 
solution of Eq. (l2. II) . Indeed, even if h and p, are not differentiable, the quantity (12.51) is 
differentiable with respect to t when a(9, to) is differentiable. It is natural to consider the 
dynamics of weak solutions because p t is a probability measure and we are interested in 
the dynamics of its moments, in particular the order parameter. In IfTOll . the existence 
and uniqueness of weak solutions of Eq. (l2. II) is proved. In this paper, the existence and 
uniqueness of solutions of a perturbed system of (12.11) . which includes (12.11) itself as a 
particular case, will be proved in Appendix B. The perturbed system is introduced in 
Sec.7 to prove the existence of a center manifold. 



3 Transition point formula and the linear instability 



A trivial solution of the continuous model (12.11) . which is independent of and t, is given 
by the uniform distribution p t (6, to) = 1 /(2n). In this case, 7]{t) = 0. This solution is called 
the incoherent state or the de- synchronous state. In this section and the next section, we 
investigate the linear stability of the de-synchronous state. The nonlinear stability will be 
discussed in Sec. 6. 
Let 

Zj(t,af):= e" fZIj9 p t (e,(o)de= e^ x(imoi) h(9)de (3.1) 
Jo Jo 

be the Fourier coefficients of p t (6,to). Then, Z (t, to) = 1 and Zj satisfy the differential 

equations 

— - = V^lwZj + -7/(0 - -r]{t)Z 2 , (3.2) 
dt 2 2 

and 

■Ji = j^najZ j + ^(n(t)Z hl -W)Zj +l ), (3.3) 

for j = 2,3, • • • . The order parameter 77(f) is the integral of Z\{t, to) with the weight 
g(to). The de-synchronous state corresponds to the trivial solution Z 7 = for j = 
1,2, Eq. (l3.1l) shows \Zj(t, to)\ < 1 and thus Zj(t,to) is in the weighted Lebesgue 
space L 2 (R, g(u))du)) for every t : 



= f \H 

JR 



As evolution equations on L 2 (R, g((o)dco), the above equations are linearized around the 
origin as 

^(VZTAuf^Z,, ,3.4) 



and 



dZj 1 — 

-l = j y[ZlMZj, (3.5) 
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for j = 2, 3, • • • , where M : qipS) h-> a>q(a>) is the multiplication operator on L 2 (R, g{oS)dco) 
and !P is the projection on L 2 (R,g(co)du)) defined to be 

Pq(a>) = I q(co)g(aj)dco. (3.6) 
Jr 

If we put Pq(co) = 1, y is also expressed as Vq(oS) = (q,P ), where ( , ) is the inner 
product on L 2 (R, g(oj)da>) defined as 



(<7i,<?2):= \ q\(co)q 2 (co)g{uS)dto. (3.7) 
Jr 

Note that the order parameter is given as rj(t) = PZx = (Z\,Pq). To determine the linear 
stability of the de-synchronous state and the order parameter, we have to investigate the 

spectrum and the semigroup of the operator T x : = V-LM + —P. 

Remark. We need not assume that the Fourier series ^ 00 Zj(t,to)e^ je converges to 
p t (8, co) in any sense. It is known that there is a one-to-one correspondence between a 
measure on S 1 and its Fourier coefficients (see Shohat and Tamarkin [40]). Thus the 
dynamics of {Zj(t, co)}™^ uniquely determines the dynamics of p t {6, to), and vice versa. In 
particular, since a weak solution of the initial value problem (12.11) is unique (Chiba IfTOl ). 
so is Eqs.d32l.d33h. 



3.1 Analysis of the operator V-LM 

Before investigating the operator T\, we give a few properties of the multiplication op- 
erator M : q(co) i— > coq(u>) on L 2 (R, g(u))du>). The domain D(M) of M is dense in 
L 2 (R, g(u))daj). It is well known that its spectrum is given by <x( Ai) = supp(g) c R, 
where supp(g) is a support of the function g. Thus the spectrum of V-lAf is 

cr( V-fM) = V^T • supp(g) = { V^U | A e supp(g)} c V^IR. (3.8) 

The semigroup e^ Mt generated by V-TA1 is given as e^ Mt q(oS) = e^^'qioj). In 
particular, we obtain 

(e^ lMt q l ,q 2 )= f e^Kw^M^, (3.9) 
Jr 

for any q x , q 2 e L 2 (R , g{oS)doS). This is the Fourier transform of the function qi(u))q 2 (ti>)g(to) 
Thus if qi(co)q 2 (u))g(a>) is real analytic on R and has an analytic continuation to the upper 
half plane, then (e^ Mt qi,q 2 ) decays exponentially as t — * oo, while if qi(co)q 2 (co)g(a>) is 
C , then it decays as 0(1 /f) (see Vilenkin Il48l ). This means that e^ Mt does not decay in 
L 2 (R,g(co)du)), however, it decays to zero in a suitable weak topology. A weak topology 
will play an important role in this paper. These facts are summarized as follows: 

Proposition 3.1. A solution of the equation (13.51) with an initial value q(co) e L 2 (R, g(u))dco) 
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is given by Zj(t, to) = e j ^ Mt q(to) = e j ^ 0Jt q(to). The quantity (e j ^ Mt qi,q 2 ) decays ex- 
ponentially to zero as t — » oo if g(to), q\{co) and qi{co) have analytic continuations to the 
upper half plane. 

This proposition suggests that analyticity of g(co) and initial conditions also plays an 
important role for analysis of the operator T\ . The resolvent (A- V-lAfT 1 of the operator 
V-TA1 is calculated as 

((A- ^PiMy l q u q2)= f l -— qi (u)li^jg(aj)da>. (3.10) 

JrA- y-lco 

We define the function D(A) to be 

D(A) = ((A - V^TAQ-^o^o) = f 1 -^g(co)dco (3.11) 

JrA - y-lco 

(recall that Po(co) = 1). It is holomorphic in C\cr(V-lA1) and will be used in later 
calculations. 

3.2 Eigenvalues of the operator T\ and the transition point formula 

The domain of T x = V^IM + § P is given by D(M) n D(P) = D(M), which is dense 
in L 2 (R, g(co)dco). Since M is selfadjoint and since P is bounded, T\ is a closed operator 
Il22ll . Let q{T\ ) be the resolvent set of 7\ and cr(T\) = C\g(T\) the spectrum. Since T\ is 
closed, Ti has no residual spectrum. Let <t p {T\) and cr c (ri) be the point spectrum (the set 
of eigenvalues) and the continuous spectrum of T\, respectively. 

Proposition 3.2. (i) Eigenvalues A of Ti are given as roots of 

D(A) = \, A 6 C\cr(V-TM). (3.12) 
(ii) The continuous spectrum of 7^ is given by 

a c {T x ) = cr( V^IM) = V^l • supp(g). (3.13) 

Proof, (i) Suppose that A e 0-^(7^ )\cr( V^lAf). Then, there exists x e L 2 (R, g(to)dto) 
such that 

,1* = ( V^IM + y!P)x, X ±0. 

Since A tr( V-bVl), (/I - V-lAfT 1 is defined and the above is rewritten as 

x = (A - V-iMT^Px 

= |(x,p )U- V^lMr'PoM. 
11 



By taking the inner product with P (a>), we obtain 

1 = - ttMT l P Q ,P Q ) = *D(A). (3.14) 

This proves that roots of Eq. (l3. 121) is in cr p (Ti)\cr( V-lAl). The corresponding eigen- 
vector is given by x = (A - ^\M)' l P {oS) = I /(A - V^M- If A e V^lR, x <£ 
L 2 (R, g(to)dco). Thus there are no eigenvalues on the imaginary axis. In particular, there 
are no eigenvalues on cr( V-TyVl). 

(ii) This follows from the fact that the essential spectrum is stable under the bounded per- 
turbation (see Kato [22]) : the essential spectrum of T\ is the same as <x( V-TA1). Since 
there are no eigenvalues on cr( V-lAI), it coincides with the continuous spectrum. ■ 

Our next task is to calculate roots of Eq. (l3.12l) to obtain eigenvalues of T\ = V-TyVt + 
jP. By putting A = x + V-Ty with x,y eR, Eq. fl3.12l) is rewritten as 

x 2 
x 2 + {co-y) 2 K 
oj-y 



y Jr 



;g{U))d(x) = 0. 



. x 2 + (a) - y) 2 ' 
The next lemma is easily obtained. 
Lemma 3.3. 

(i) A satisfies Re(A) > for any K > if it exists. 

(ii) If K > is sufficiently large, there exists at least one eigenvalue A near infinity. 

(iii) If K > is sufficiently small, there are no eigenvalues. 

Proof. Part (i) of the lemma immediately follows from the first equation of Eq. fl3.15l) : 
Since the right hand side is positive, x in the left had side has to be positive. To prove part 
(ii) of the lemma, note that if \A\ is large, Eq. fl3.12l) is expanded as 

Thus the Rouche's theorem proves that Eq. (l3. 121) has a root A ~ K/2 if K > is suffi- 
ciently large. To prove part (iii) of the lemma, we see that the left hand side of the first 
equation of Eq. (l3. 151) is bounded for any xj G R. To do so, let G{oS) be the primitive 
function of g(a>) and fix 6 > small. The left hand side of the first equation of Eq. fl3.15l) 
is calculated as 



xg(a))dco 
R x 2 + (a>-y) 2 



L 

r°° xg(co)duj C y 6 xg{(jS)d(x> C y+6 xg(a>)du) 

J y+S x 2 + (oj-y) 2 J_„ x 2 + (co-y) 2 J y _ 

r°° xg(co)duj r yS xg(o))da> 

Jy+gx 2 + (to-y) 2 J_ a 



x 2 + (to - y) 2 

6 



x 2 + (to - y) 2 

x* + & K ~ w ' u0; " 5)) + X- 5 (x 2 + (co-y) 2 ) 2 ' 



C y+S 2x(co - y) 
(G(y + 6) - G(y - 6)) + y -^G(co)du. 

Jv-S 
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The first three terms in the right hand side above are bounded for any x, y e R. By the 
mean value theorem, there exists a number £ such that the last term is estimated as 



2x(to - y) 



f ( 2 17 y ~^2H G{ 

Jys (x 2 + (to- y) 1 ) 1 



2xco 

(x 2 + CO 2 ) 2 



(G(y + to) - G(y - to))dto 



- f 

Jo 

= (G(y + 0) - G(y - 0)) + Wy + *) - G(y - 5)) £ j^jf**- 

(3.16) 

Since G is continuous, the above is calculated as 

x x 



(G(y + S)-G(y-6)) 



yX 2 + £ 2 JC 2 +(5 2 / 

If £ ^ 0, this is bounded for any x, y e R. If £ = 0, Eq. (l3.16l) yields 

I (5T^ (G(y + w) " ° (y - ^ = <G(y + * " ° (y - ™ [WT^ 

Since G(to) is monotonically increasing, we obtain 

G(y + to)- G(y - to) = G(y + 6)- G(y - 6) 

for < to < 6. In particular, putting to = gives G(y + 6) - G(y - 6) = 0. Thus 
G{y + to) - Giy - to) = for < to < 5. This proves that the quantity (13.161) is zero. Now 
we have proved that the left hand side of the first equation of Eq. (l3.15l) is bounded for any 
jcjgR, although the right hand side diverges as K — » +0. Thus Eq. (l3.12l) has no roots if 
K > is sufficiently small. ■ 

Lemma 3.3 shows that if K > is sufficiently large, the trivial solution Z\ = of the 
equation dZJdt = T\Z\ is unstable because of eigenvalues with positive real parts. Our 
purpose in this section is to determine the bifurcation point K c such that if K < K c , the 
operator T\ has no eigenvalues, while if K exceeds K c , eigenvalues appear on the right 
half plane (K c should be positive because of Lemma 3.3 (iii)). To calculate eigenvalues 
A = A(K) explicitly is difficult in general. However, since zeros of a holomorphic function 
do not vanish because of the argument principle, A(K) disappears if and only if it is ab- 
sorbed into the continuous spectrum cr( V-TA1), on which D(A) is not holomorphic, as K 
decreases. This fact suggests that to determine K c , it is sufficient to investigate Eq. (l3. 121) 
or Eq. (l3.15l) near the imaginary axis. Thus take the limit x — » +0 in Eq. (l3.15l) : 

2 

;g(to)dtO = — , 



lim f 

C to — y 
{ x^+o J R x 2 + (to - y) 2 



x 2 + (to-y) 2 ° K 

(3.17) 
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Fig. 3: A schematic view of behavior of roots A of Eq. (l3.12l) when K decreases. Thick 
lines denote the continuous spectrum. As K decreases, eigenvalues A\,Az, ■ • ■ converge to 
V-Tyi, V-Ty2, • • • and disappear at some K = Ki,K 2 , • • • , respectively. 



These equations determine Kj and y ; such that one of the eigenvalues A = Aj(K) converges 
to yf-Tyj as K — » K } ■ + (see Fig©. To calculate them, we need the next lemma. 

In this paper, the probability density function g(co) is called piecewise continuous 
when it satisfies following: If g(co) is discontinuous at coq, then g(coo + 0) and g(a>o - 0) 
exist, and for any s > 0, there exists 6 > such that when -6 < co - a>o < 0, \g(of) - 
g(oj - 0)| < e and when < a> - a> < 8, \g{co) - g(co + 0)| < e. If g{oS) is piecewise 
continuous except for a measure zero set, then we also call it piecewise continuous. Note 
that when g((o - 0) = g(o» + 0) # g(a>o), we can modify g(a>) on a point {co } so that it 
becomes continuous at a> . Thus if g is discontinuous at co , we can assume without loss 
of generality that g(a>o - 0) ^ g(a> + 0). 

Lemma 3.4. (i) Suppose that g(a>) is piecewise continuous. If one of the eigenvalues 
A(K) satisfies A(K) —> V-Tyy as K — > Kj. Then, g(a>) is continuous at co = yj. 
(ii) If g(co) is continuous at to = y, then 

lim n f TTT ^g(co)dco = 7rg(y). (3.18) 

x^+o J R x z + (co - y) 1 

Proof. To prove (i), suppose that g(co) is discontinuous at co = without loss of gener- 
ality. We will show that a root A(K) of Eq. (l3.12l) can not approach to zero for any K > 0. 
Put g(+0) = h + , g(-0) = h- and h + ± h_. For any s > 0, there exists 6 > such that if 
-5 < co < 0, then \g{co) -h-\ < e and if < to < 8, then \g(to) - h+\ < s. For Eq. (l3.12l) . 
we suppose \A\ = \x + V-Ty| < 5 and y > 0. The case y < is treated in a similar manner. 
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We calculate D(A) as 



D(A) 



A- V-l. 



-dio + 



CO 



A - V-T 



-dio 



CO 



+ 



r 6 1 

=- (g(co) - h + + h + ) da> + 

Jo A - y-lio 



r° 1 

— - (g(co) - h_ + h_) dio 

J-sA- y-lio 



Js A - yhAco Js A + V-l 

r 6 1 

+ =- (g(cu) - h + ) da) + 

Jo A - y-lco 



-dio 



CO 



r 6 i 

=- (g(-co) - h_) dio 

Jo A + y-lio 



+ 



1 



(h+ — hJ) dio + h 



-f 

Jo 



dio 



+ h 



f 

Jo 



dco 



lo A - V-lw V Jo A - V-Tw Jo A+ \[^\oj 

Since \A\ < 6, there exists a positive number M, which is independent of A, such that 

g(±co) 



. (3.19) 



f 

J s 



A+ V 3 ! 



-dio < M. 



10 



Thus \D(A)\ is estimated as 
\D(A)\ > \h 



Jo \A- y-lco\ Jo 



dio 



h- log 



\a- V-T< 
a+ V 3 !^ 



- s 



\A - V-M Jo \A + V-M 



dio 



2M. 



Since y > 0, \A - V-lw| < \A + V-lw| when <x> > 0. This shows that 



\D(A)\>(\h + -h4-2s) 



r 6 dio _ 

Jo \A - V-M 



/?_ log 



a- V-T<sJ 



-2M. 



(3.20) 



The right hand side tends to infinity as A — » if 2s < \h+ - h-\. This proves that as K 
decreases, a root of Eq. (l3.12l) can not approach to zero for positive K. 

(ii) The formula Eq. (13. 1 81) is famous and given in Ahlfors 0. ■ 

Recall that the second equation of Eq. (13. 171) determines imaginary parts to which A(K) 
converges as Re(A(K)) — » +0. Let y u y 2 , • • • be roots of the second equation of Eq. (l3.17l) . 
at which g(io) is continuous. Since g(co) is continuous at to = yj, substituting it into the 
first equation of Eq. fl3.17l) yields 



1 ng(yj) 



, J = 1,2, 



(3.21) 



which gives the value such that Re(/l(^)) — » as K —> K } ■ + 0. Now we obtain the next 
theorem. 
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Theorem 3.5. Let y\,yi, • • • be roots of the second equation of Eq. (l3.17l) . Suppose that g 
is continuous at y\, yi, ■ ■ ■ (this is true when g is piecewise continuous). Put 



K c := MKj = . (3.22) 

j n sup jg(yj) 

If < K < K c , the operator T\ has no eigenvalues, while if K exceeds K c , eigenvalues of 
T\ appear on the right half plane. In this case, the trivial solution Z\ = of Eq. (l3.4l) is 
unstable. 

In general, there exists Kjr such that T\ has eigenvalues when K c < K < but they 
disappear again at K = fC c , i.e. the stability of the trivial solution Z\ = may change 
many times (see Example 3.9). Such Kf } is one of the values K/s. However, if g{to) is an 
even and unimodal function, it is easy to prove that T\ has eigenvalues on the right half 
plane for any K > K c , and they are real as is shown in Mirollo and Strogatz [|32l . Indeed, 
the second equation of Eq. (l3.15l) is calculated as 

f a> — y I to 

= L^(co- y y 8(co)daj = i -r—i(z(y + ^-z(y-^- 

If g is even, y = is a root of this equation. If g is unimodal, g(y + co)- g(y -co) > when 
y < 0, co > and g(y + to) - g(y - to) < when y > 0, to > 0. Hence, y = is a unique 
root. This implies that an eigenvalue should be on the real axis, and (K, y) = (K c , 0) is a 
unique solution of Eq. (13 .171) . As a corollary, we obtain the transition point (bifurcation 
point to the synchronous state) conjectured by Kuramoto [|26l : 

Corollary 3.6 (Kuramoto's transition point). Suppose that the probability density func- 
tion g(co) is even, unimodal and continuous at y = 0. Then, K c defined as above is given 
by 

K c = (3.23) 
ng(0) 

When K > K c , the solution Z\ = of Eq. (l3.4l) is unstable. In particular, the order 
parameter -q(t) = (Z\,Pq) is linearly unstable. 

Example 3.7. If g(co) is given by 

f 1 (-1 < to< -1/2) 
g(co) = \ 1/2 (l/4<w<5/4) (3.24) 
{ (otherwise), 

(see Figl4l(a)). then roots of the second equation of Eq. (l3.17l) are y 0.706, -0.044, 0.567. 

Thus K c is given as 

2 2 
K = _ _ 

c ng(-Q.lQ6) n 
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(a) 



(b) 



> 
















-1 -1/2 


1/4 


5/4 < 




-1 -1/2 




1/2 1 co 



Fig. 4: Probability density functions g(co) for (a) Example 3.7 and (b) Example 3.8. 



More generally, suppose that g(co) is of the form 



( c\ (a\ < co < b\) 

C2 {a,2 < co < b-i) 

: : (3.25) 

c„ (a n <co< b n ) 

(otherwise), 



where c ; + and ci\ <b\ < a 2 < ■ ■ ■ <b n are constants. In this case, we can prove that the 
second equation of Eq. (l3.17l) has roots yj on the intervals (cij, bj) for j = 1, • • • , n. Thus 
K c is given by 

K c = . (3.26) 

^maxi< 7 <„ cj 

Indeed, the left hand side of the second equation of Eq. (l3.17l) is calculated as 

(y-aj 



I co — y 

lim — — —g(co)dco = ) c. log 

^jRX 2 + (co-y) 2SK U 



bj-y 



This tends to -oo as y —* a.j + and to +oo as y —* bj - 0. Now the intermediate-value 
theorem proves the desired result. 

Example 3.8. If g(co) is given as Fig|4] (b), the second equation of Eq. d3.17l) has roots 
y = and y ~ +0.7459. Thus K c is approximately given by 

2 2 
Kc = ^(0.7459) = 7T - 1.9671 ' (3 ' 27) 

It is notable that the value 0.7459 is slightly different from 0.75 which gives the maxi- 
mum of g{co). This result gives a counter-example for the Kuramoto's argument K c = 
2/ \n max g(coJ) for bimodal distributions II261I . Such a counter-example is also given by 
Martens et al. 



Example 3.9. We give an example that the stability of Z\ =0 changes many times. 
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Suppose that g(a>) is given by 



g(to) = 



-O.lw + 3.8 

3 

O.lw + 3.8 

7 





(-3 < co < -2) 
(-2<co< -1) 
(-1 <co< 1) 
(1 < co < 2) 
(2 < co < 3) 
(otherwise), 



(3.28) 



see Figj5](a). Behavior of eigenvalues A(K) is represented as FigJ5](b), in which arrows 
indicate the direction that K increases. Solutions of the second equation of Eq. (l3. 171) are 
approximately given by 



y ±l = +2.1032, y ±2 = +1.6882, y ±3 = +1.4046, y 4 = 0. 



This means that Z\ = is unstable if and only if 

2 



0.0909 



K r < K < K 



(2) 



0.1604 



and 



0.1616 



< K < oo. 



(3.29) 



(3.30) 



(3.31) 



Thus Zi = recovers the stability in the small interval 0.1604 < K < 0.1616. Indeed, 
an eigenvalue emerges at V-lyi when K = K c , it goes to the right side as K increases 
for a moment, but it again approaches the imaginary axis and disappears at V-Ty2 when 
K = Kr c '. After that, another eigenvalue emerges at V-Ty3 when^f = ^3. This eigenvalue 
remains on the right half plane for any K > (by using the theory of resonance poles 
developed in Sec.4, we can show that the eigenvalue disappearing at V-ly2 is connected 
to the eigenvalue emerging from V-IV3 in the second Riemann sheet of the resolvent of 
Ti, see Fig©. 



Remark. Although we suppose in this paper that the distribution g(co) for natural fre- 
quencies is a probability density function, it is interesting to consider a more general case 
that natural frequencies are distributed according to a probability measure g, not a func- 
tion. In this case, we can prove that K c = if and only if there is a point y e R such 
that g({y}) £ 0. The proof is done in a similar way to that of Lemma 3.3 (iii), in which 
g(u>)dto and G{co) are replaced by dg and g((-oo, co]), respectively. If g({y}) + 0, then 
G(y + 0) ^ G(y - 0). As a result, the first term in the right hand side of Eq. (l3.16l) is 
estimated as 0(1 /x) as x — » +0. This and Rouche's theorem prove that Eq. (l3. 121) has a 
root A for any small K > satisfying Re(A) ~ 0(K) as K — » +0. Thus the operator T\ has 
eigenvalues on the right half plane for any K > 0. This means that if a lot of oscillators 
share the same natural frequency, they may synchronize for any small K > 0. 
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(a) 



(b) 



g(eS) 




-3 -2 - 



2 



3 



CO 



3 



Fig. 5: (a) A probability density function g(a>) for Example 3.9. (b) A schematic view of 
behavior of eigenvalues A(K) for Example 3.9. 

4 Linear stability theory 

Theorem 3.5 shows that K c is the least bifurcation point and the trivial solution Zi = of 
Eq. (l3.4l) is unstable if K is slightly larger than K c . If < K < K c , there are no eigenvalues 
and the continuous spectrum of T\ lies on the imaginary axis: cr(Ti) = cr( V-TA4). In 
this section, we investigate the dynamics of Eq. (l3.4l) for < K < K c . We will see that the 
order parameter 77(f) may decay exponentially even if the spectrum lies on the imaginary 
axis because of the existence of resonance poles. 

4.1 Resonance poles 

Since V-TA4 has the semigroup e^ Mt and since V is bounded, the operator T\ = 
V-TA4 + —V also generates the semigroup e T[t (Kato |[22T0 on L 2 (R,g(u>)da>). A so- 
lution of Eq. (l3.4l) with an initial value (p(co) e L 2 (R, g(u))daj) is given by e Txt (j){oj). The 
semigroup e T,t is calculated by using the Laplace inversion formula 



for t > 0, where x > is chosen so that the contour (see Figj7](a)) is to the right of the 
spectrum of T\ (Hille and Phillips ETI . Yosida Il49l ). The resolvent (A -Ti)~ l is given as 




(4.1) 



follows. 



Lemma 4.1. For any <p(oS), i^(a») e L 2 



(R, g(u))du)), the equality 



((A-T.T^tfr) 

((A-yTlMr l (f>,ij/) + 



1 - KD(A)/2 



K/2 



((A - yTlMr^, P )((A - yTlMT'Po, <A)(4.2) 
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holds. 
Proof. Put 

R(A)(f> := (A - ri)"V = 0* - V-FM - 

which yields 



(A - V-iM)R(A)<f> = <f> + jPR(A)<f> 

= <f,+ j(R(A)<f>,P )P . 

This is rearranged as 

R(A)<f> = (A- "nPIAO'V + j^W' p o)0* - V^lAD-'Po. (4.3) 
By taking the inner product with Po, we obtain 

(R(A)cf>,Po) = ((A - V=TAO-V,^o) + |(«(A)0,P O )O(A). 



This provides 



(pw,p ) = (u- V^TMrv^o). 



Substituting it into Eq. (l4.3l) . we obtain Lemma 4.1. 

Eq. (l4.1l) and Lemma 4.1 show that (e T[t (f>, iff) is given by 

Y r>x+ V-Ty 

(e r 'V,<A) = bm = e At ((A - V-LM)~V, 

y^°° 2ny-lJx-^y x 

+ l-KD(A)/2 ((A ~ A/ZTA1) " 1 ^' Po)((A ~ V^AO -1 ^, if,j)dA. (4.4) 

In particular, the order parameter 77(f) = (Z\ , P ) for the linearized system (13.41) with the 
initial condition Z y {§, co) = (p(to) is given by 7]{t) = (e Ti '(p, P ). 

One of the effective ways to calculate the integral above is to use the residue theorem. 
Recall that the resolvent (A - Pi) -1 is holomorphic on C\cr(Ti). When < K < K c , 
T\ has no eigenvalues and the continuous spectrum lies on the imaginary axis : <t(Pi) = 
cr( V^Al) = • supp(g). Thus the integrand e At ((A - T x y 1 ^, </0 in Eq.(l4~4l) is holo- 
morphic on the right half plane and may not be holomorphic on cr(T\). However, under 
assumptions below, we can show that the integrand has an analytic continuation through 
the line <r(T]) from right to left. Then, the analytic continuation may have poles on the 
left half plane (on the second Riemann sheet of the resolvent), which are called reso- 
nance poles (381. The resonance pole A affects the integral in Eq. (l4.4l) through the residue 
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theorem (see Fig|7](b)). In this manner, the order parameter ?](t) can decay with the expo- 
nential rate e Re(/i)f . Such an exponential decay caused by resonance poles is well known in 
the theory of Schrodinger operators [38], and for the Kuramoto model, it is investigated 
numerically by Strogatz et al. H4ll and Balmforth et al. flU. 

At first, we construct an analytic continuation of the function F (A) := ((A-T^y 1 ^, if/). 
In the next lemma, we consider {{A — T{)~ l <f), iff) instead of {{A — TyY l <f>, iff) so that the com- 
plex conjugate of iff(co) does not appear in the right hand side of Eq. (l4.5l) (by the definition 
of the inner product, ( • , iff) involves the conjugate of if/, while ( • , iff) does not). 

Lemma 4.2. Suppose that the probability density function g(co) and functions (f>(u)), iff(co) 
are real analytic on R and they have meromorphic continuations to the upper half plane. 
Then the function F (A) := ((A - Ti)' 1 ^, iff) defined on the right half plane has the mero- 
morphic continuation Fi(A) to the left half plane, which is given by 

F^A) = ((A- V r lMrV,^) + 2^(-V r Ti)(A(-V r T^(-V r Ti) 

+ — ^QU, <f>]Q[A, if/}, 

1 - KD(A)/2 - nKg(- yTlA) 

where Q[A, cf>] is defined to be 



(4.5) 



Q[A, cf>] = ((A - V^lMrV, Po) + 2ng(- V^U)</>(- V^U). 



(4.6) 



Note that Q[A, ■ ] defines a linear functional for each A e C. Actually, we will define 
a suitable function space in Sec. 5 so that Q[A, ■ ] becomes a continuous linear functional 
(generalized function). 

Proof. Define a function F(A) to be 

f(a) = { ((A ~ ^- [M y l ^t ) ( Re ^ > °)> 

1 ((A - V^lMrV, iff) + 2tt0(- V-U)<A(- V r T^)g(- V-U) (Re(^) < 0). 

(4.7) 

By the formula (13.181) . we obtain 



lim ((A - ^P\MY ] d>,if/) - lim ((A - V-fAl) -1 ^ <A) 

Re(i)->+0 ReUH-0 

= 2ncf>(lm(A)) ■ if/(lm(A)) ■ g(lm(A)), 



(4.8) 



which proves that limRey^+o F(A) = limRe^^.o F(A). Therefore, if we show that F(A) 
is continuous on the imaginary axis, then F(A) is meromorphic on C by the Schwarz's 
principle of reflection. To see this, put <p{oj)if/{oS)g{oS) = q{co). By the formula (13.181) . 



lim 



= lim 



r°° i 
J-cA - V 3 ! 

f 

U — o 



q(co)du> 



x ^+° J-oo x 2 + (oj- y) 2 
nq(y)-n^V(y), 



a> - y 



x 2 + (to - y) 2 



q(a))dco 
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where A = x + yf—ly and V(y) is the Hilbert transform of q defined by 



VO0 = P-v.- f -q(y-t)dt, (4.9) 
n J-co t 

see Chap. VI of Stein and Weiss ATI . Since q(y) is Lipschitz continuous, so is V(y) 
(Thm.106 of Titchmarsh Il45ll ). This proves that lim v ^ +0 / (A - V-Ta>) 1 q(oS)dto is con- 
tinuous in y. Therefore, F(A) is meromorphic on C. Now we have obtained the meromor- 
phic continuation of ((A - V-lAfT i <p, if/) from right to left. Applying this to Eq. (l4.2l) . we 



obtain the meromorphic continuation of F (A) as Eq. (l4.5l) . ■ 
Eq. (l4.5l) is rewritten as 

FlW . m — ( 

1 - KD(A)/2 - nKg(- V-U)\ 
(2/ K - D(A))((A - V^lMr>,^) + ((A - V^lMrV^o) ■ (U ~ ^lMy l if/,P ) 
+ 2ng(- V-U)(-|(£(- V-U)0<- V-U) - D(A)(f>(- V-U)iK- V-U) 

-((A - V-fAirV,^) + ((A - V-LMrV, At) • if/(- V-U) 
+(U - V^IAO'V. A,) • 0(- V^IA))!. (4.10) 

This expression shows that poles of g are removable. Therefore, poles of F\(A) on the left 
half plane and the imaginary axis are given as roots of the equation 

D(A) + 2ng(- V^U) = |, Re(A) < 0, 

K 2 ( 4 - U ) 

lim D(A) = lim (D(A) + 2tt^(- V-U)) = -, Re(A) = 0, 

, Re(A)^+0 Re(A)->-0 v ' K 

and poles of the functions 0(- V-U) and </r(- V-T/l)- To avoid dynamics caused by a 
special choice of (p and if/, in what follows, we will assume that continuations of (p and if/ 
have no poles. 

Definition 4.3. Roots of Eq. (l4. Ill) on the left half plane and the imaginary axis are called 
resonance poles of the operator T\ . 

Since the left hand side of Eq. (l4. Ill) is an analytic continuation of that of Eq. (l3.12l) . 
at least one of the resonance poles is obtained as a continuation of an eigenvalue A(K) 
coming from the right half plane when K decreases from K c (see Fig. |6] and a comment 
in Example 3.9). However, if g(A) has an essential singularity, there exist infinitely many 
resonance poles in general, which are not obtained as continuations of eigenvalues. 

We want to calculate the Laplace inversion formula (14.41) by deforming the contour as 
FigJ7](b), and pick up the residues at resonance poles. We should show that the integral 
along the arc C 4 converges to zero as the radius tends to infinity. For this purpose, we 
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Kc<K 



0<K<K c 



Fig. 6: As is discussed in Sec. 3, an eigenvalue A(K) disappears from the original complex 
plane at K = K c . But it still exists as a resonance pole on the second Riemann sheet of the 
resolvent. 



have to make some assumptions for growth rates of <p(A) and \]/(A) as \A\ — > oo. Since 
suitable assumptions depend on the growth rate of g(A), we calculate the Laplace inversion 
formula by dividing into two cases: In Sec. 4. 2, g(co) is assumed to decay faster than any 
exponential functions, such as the Gaussian distribution. In Sec.4.3, we consider the case 
that g(co) is a rational function. 

Remark. In Lemma 4.2 and Def.4.3, g(co) is assumed to be analytic on the real axis. 
If g(co) is piecewise analytic on the real axis such as Examples 3.7 to 3.9, we can show 
that an analytic continuation of the resolvent has logarithmic singularities: if g(co) is not 
analytic at co = y e R, then V-Ty is a logarithmic singularity of F (A). In this case, there 
exist infinitely many Riemann sheets and resonance poles exist on every Riemann sheets 
in general. Resonance poles on the n-th Riemann sheet are given as roots of the equation 



In this paper, we do not deal with such problems. 

4.2 Gaussian-type case 

In this subsection, we suppose that the probability density function g(co) for natural fre- 
quencies satisfies that 

(Al) g(a>) is real analytic on R and has a meromorphic continuation g(A) to the upper 
half plane. 

(A2) for any y > 0, there exists a positive number C such that 



D(A) + 2nng{- V-U) = -, Re(i) < 0. 



\g(co)\ < Ce~ 



coeR. 



(4.12) 



(A3) there exists a positive number 6 such that in the angular domains 



|argO*)| < 6, |argU) - n\ < 8, 



(4.13) 



near the real axis, g(A) satisfies Eq. (l4. 121) for any y > 0. 



23 



Note that since g(a>) is real- valued on R, (Al) implies that g(a>) also has a meromor- 
phic continuation to the lower half plane due to the Schwarz's principle of reflection. The 
assumption (A3) implies that there are no poles of g in the angular domains (14.131) . These 
assumptions are satisfied, for example, when g is of the form g(co) = C'e P(w) with a poly- 
nomial P(a>) = -a/" 1 + a\co 2m ~ l + ■ ■ ■ + aim and some constant C", such as the Gaussian 
distribution. Obviously (A3) includes (A2), however, they are divided to emphasize that 
(A3) is an assumption for an analytic continuation of g. The assumption (A2) holds while 
the assumption (A3) fails, for instance, if g(co) decays as exp(-exptt> 2 ) on R. 

To make assumptions for (p and iff, we prepare a certain function space. Let Exp + (J3, ri) 
be the set of holomorphic functions on the region C„ := {z e C | Im(z) > -I In] such that 
the norm 

11011^:= sup e-^\<p(z)\ (4.14) 

Im(V)>-l/n 

is finite. With this norm, Exp + (/?, ri) is a Banach space. Let Exp + (/3) be their inductive 
limit with respect to n = 1, 2, • • • 

Exp + (/?) = limExp + (# ri) = (J Exp + (£, ri). (4.15) 

;i>l n>\ 

Thus Exp + (/3) is the set of holomorphic functions near the upper half plane that can grow 
at most the rate e^ 1 " 1 . Next, define Exp + to be their inductive limit with respect to = 
0,1,2,--- 

Exp + = HmExp + (/8) = \jExp + (j3). (4.16) 

p>o /?>o 

Thus Exp + is the set of holomorphic functions near the upper half plane that can grow 
at most exponentially (see also Fig.©; (f>(z) in Exp + satisfies \\<p\\^ n < oo for some fi, n, 
and such /3 and n can depend on <p. Topological properties of Exp + will be discussed in 
Sec. 5. 2. In this section, the topology on Exp + is not used. 

Finally, we make the assumption for the parameter K. Although we are interested in a 
positive K from a physical viewpoint, any complex number is considered in the proof of 
the next theorem. In the theory of functions, if a function f(A) converges toweC along 
a path y tending to oo, w and y are called an asymptotic value and an asymptotic path, 
respectively. Let "K c C U {oo} be the set of numbers K such that 2/K is an asymptotic 
value of the function 2ng(A). Since D(A) — > oo as \A\ — > (see Claim below), 2/K is 
an asymptotic value of 2ng(A) if and only if it is an asymptotic value of the function 
D(A) + 2ng{- V-Ll). Define a subset TCo c 7C as follows: K e <K if and only if there is 
an asymptotic path y included in the angular domain n/2 + 6 < arg(/l) < 3^/2 - 6 such 
that 

D(A) + 2ng(- V^Ll) = | + o(-j-), \A\ oo, (4.17) 

holds along y for every n = 0, 1, 2, • • • , where 6 > is the constant appeared in the 
assumption (A3). That is, if the function D(A) + 2ng(- V-Ll) converges to a value 2/K 
faster than any polynomials along a path y in the angular domain, then K e ( Kq. Put 
%\ = 7Co U {0} U {oo}. It is well known as the Denjoy-Ahlfors theorem that if g(A) is 
an entire function of order p, then the number of asymptotic values is at most 2p ll50ll . 
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In this case, TCi is a finite set. If g(A) is an entire function of infinite order or g(A) is 
a meromorphic function, the set of asymptotic values 'K can have the cardinality of the 
continuum (see Gross lfT9l and Eremenko lfT6l for such examples; see also Bergweiler and 
Eremenko [0 for a condition that a meromorphic function has at most finite asymptotic 
values). However, since the condition (14.171) and the condition that g(A) is an integrable 
real-valued function on R are too strong, it seems that "K\ = {0, oo} in most cases (actually, 
the author does not know an example of a probability density function g(a>) satisfying 
(Al) to (A3) whose cardinality of "K\ is larger than two). 
The main theorem in this section is stated as follows. 

Theorem 4.4. Suppose that K £ "K\ and g(co) satisfies the assumptions (Al) to (A3). For 
any (f>,if/ e Exp + , there exists a positive number t such that the semigroup e T,t satisfies 
the equality 

M oo 

(e Tlt <f>,^) = J^SAt,^,^ + Y j R n [t,(f>,Me A « t , (4.18) 

for t > t , where ■ • ■ ,£m are eigenvalues of T\ on the right half plane, S n [t, <p, if/]e^ nt 
are corresponding residues of F (A)e At , and where Aq,Ai, - ■ ■ are resonance poles of T x 
such that \A \ < \Ai\ < • • • , and R n [t, cp, if/]e A "' are corresponding residues of Fi(A)e At . In 
particular when < K < K c , it is written as 

oo 

(e T ^,^) = Y,Rn[t,(f>A]e A " t . (4.19) 

If A n is a pole of order k, then R n is a polynomial of degree at most k - 1 in t. When 
< K < K c , all resonance poles lie on the left half plane. In particular, the order parameter 
77(f) = (e T[t (f>,Po) for the linearized system (13.41) decays to zero exponentially as t — > 00. 

To prove Eq. (l4.18l) . we need the next estimate. 

Claim. For any (f>, ij/ £ Exp + , ((A - V-lAI) -1 ^ if/) tends to zero as \A\ —> 00. 

Proof. This is a well known property of the resolvent when A is included in the resolvent 
set. Indeed, by the definition, we have 

((A - V^lAtrV^) = f -^(p(a>)if/(a>)g((o)d(o 

JrA- yf-lco 

C x 1 — C oj — y 
= J*** + (o, - yf qia))dU) + ^ Lc* + (a> - yY q{(i))da) > 

where we put A = x + V-Tv and q(oS) = (f)(a))if/(oj)g(a>). Since <p,ifj e Exp + , there exist 
Ci, C 2 ,/3!,/3 2 > such that |0(z)| < C x eP M , \if/(z)\ < C 2 e^. Thus we obtain 

\((A - V^lMrV,^)| <dC 2 f l -^e^ +/3M g(aj)daj. 

Jr\A - ^/-\tL>\ 

By the assumption (A2) for g(aj), this integral exists when A i V-1R. When A £ 
V-TR, l/\A - V— Tco»| — » as \A\ —* 00 uniformly in a> e R, which proves that \((A - 
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^^^^My l (p,^J/)\ — > 0. In particular, an upper bound of \((A - ^[^\MY l <p, if/)\ is deter- 
mined by C\ , C 2 ,fii and /3 2 . 

To estimate the case A e V-lR, take the limit x — » +0 to yield 

lim ((/I - V-rMr V, <A) = ^(y) - * V-lV(y), 

where V(y) is defined by Eq.(l4T9l>. Due to the assumption (A2), |g(y)| < C 1 C 2 e (/Jl+/32)|v| g(y) 
tends to zero as y —> +oo. Next, let us prove V(y) —* as y —> +oo. Due to the assumption 
(A3), \q(A)\ tends to zero as \A\ — * oo in the angular domain (14.131) . This fact and Cauchy's 
theorem allow us to rewrite V(y) as 



V(y) 



i r°i 

* Jo f 



( ? (y-0-9(y + 0)^ 



1 

X Jo 



-(^(y-f)-9(y + 0)* + - I -(9(y-0-9(y + 0)*, 



where a is a sufficiently small positive number. By the assumption (A3), the second term 
above tends to zero as y —* +oo. For the first term, there exists a number /ceC such that 



f^'l r y ~ ia da 

-(q(y-t)-q(y + t))dt = 2 -f{y + Kt)dt. (4.20) 

Jo * Jo ay 

Since q(z) is holomorphic on the upper half plane except for poles of g(z), Cauchy's 
formula is applied to provide 

^(Z ) = —^= f —^— I q(z)dz = J- f "4^0 + flg^VS 

with zo = y + Kt. By the assumption (A3), q(y + Kt + ae^ 6 ) — » as y — » +oo uniformly 
in < 6> < 27r. Therefore, Eq. (l4.20l) tends to zero as y — » +oo. This completes the proof 
of Claim. ■ 
From the above proof, we see that ((A- V-1A!) _ V> i/0 has an upper bound determined 
by constants C\,C%,P\ and /? 2 . 

Proof of Thm.4.4. Let us prove Eq. (l4.18l) . Take a sequence {r n }™ =Q of positive numbers 
as follows: if 2/K is not an asymptotic value of the function D(A) + 2ng(- V-T/l), then 
there exist a positive constant A and a sequence {^„)~ =0 with r„ — > oo such that 



i_^ g( _V^T^)-|D(^) 



(4.21) 



for /t = r n e^ 6 , n/2 + 6 < 6 < 3n/2 - 6. Even if 2/K is an asymptotic value of the 
function D(A) + 2ng(- V-T/l), because of Eq. (l4.17l) . there exist a positive constant A, an 
integer m > and a sequence {r„}^ with r„ — » oo such that 



l-7rZ"g(-V=TA)-|D(A) 



> 



(4.22) 
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holds for A = r„e^ e , n/2 + S<8< 3n/2 - 6. 

Recall that eigenvalues of T\ are roots of the equation D{A) = 2/K. Since D(A) — » 
as \A\ — » oo, the number of eigenvalues on the right half plane is finite, say £o,--- ,£m- 
Take a positive number d > so that Re(£ ; ) < d for j = 0, • • • , M. With these d and r n , 
take paths C\ to C 6 as are shown in FigJ7](b): 

Ci = {d+ V-Ty | -r n <y< r n ), 

C 2 = {x+ V^Ir„ | < x < d}, 

C 3 = {r n e ^ \n/2<6< n/2 + 6}, 

C 4 = {r n e^ B \n/2 + S<6< 3n/2 - 6}, 

and C5 and are defined in a similar way to C3 and C%, respectively. We put C ( ' !) = 
2 j=i Cj. 

Let Aq,Ai,--- , A#( n ) be resonance poles inside the closed curve C (,!) . By the definition 
of r n , there are no resonance poles on the curve C (,l) . By the residue theorem, we have 

_ ' m m) \ r r 

2n V^T Yj S £*> 4>> ^ + Z R J [t ' ^ = I e At F {A)dA + I e M F x (X)dA, 

when r n is sufficiently large so that C (n) includes all eigenvalues £o,--- ,%m- Since the 
integral J c e At Fo(A)dA/(2ny^l) converges to {e Tlt <p, \f/) as n — > 00, 



#00 

(e r < >, <A) = Y! S;[f, 0, <A]^ + lim Y /?,•[*, 0, <^ 

7=0 

i=u m f 

V-l »^°° Jc 



V-l Jc 



2^^ 



c 6 +c 2 2n V-l 



e*Fi(/l)<M. (4.23) 



C3+C4+C5 
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It is easy to verify that the integrals along C 2 , C 3 , C 5 , C 6 tend to zero as n — » oo because 
of the assumption (A3). For example, the integral along C 2 is estimated as 



e At F (A)dA 



e (^V=ir„) fFo(x+ VZT r „)djc 

< |F (jc + yf-Lr n )\dx, 

Jo 



where Fo is given as (14.21) . By using Claim above, the integral along C 2 proves to be zero 
as n — » oo (r„ — » oo). The integrals along C 3 ,C5,C6 are estimated in a similar manner. 
The integral along C 4 is estimated as 



\ e At F l {A)dA 
Jc 4 



J~3tt/2-<S 
r„e r " fcos0 | J F 1 (r„e VrTe )|^ 
ti72+(5 









IF^e^l 


pnl2 


< 


max 




2r n e- r " tsin6 de 




;r/2+(5<#<37r/2- 


-5 








F^e^l 


rn/2 


< 


max 




J 2r n e~ 2r " t6ln de 




jt/2+6<8<3jt/2- 


-6 


< 


max 




\Fi(r n e^ e )\ 


X ^ e -2r„tS/7T _ e -r„t 




ff/2+5<e<3ff/2- 


-6 





Since \(A - V^M)" 1 ! -» as \A\ oo, F^A) given by Eq. (l4.10l) is estimated as 



\F\{A)\ < 



p +|A +D 2 0(-yrT /1) +jD3 ( A( _yrT /1 ) +Z ) 4 0(_yzu)^(_yzr /t )| . ^(-yrr^i 

H-Tr^-V^T^-fDU)! 

where Z) to D A are some positive constants. Since (p, ip e Exp + , there exist C 1 ,C 2 ,fi l ,/3 2 > 
such that \(f>(z)\ < C x eP M , \xf/{z)\ < C 2 e ft|z| : 

D + (D i + D 2 C x e M + D 3 C 2 ^ 2W + D A C y C 2 e ( ^ M ) ■ |g(-V-U)| 

\Fi(A)\ < i — . 

\l-7TKg(-V-iA)-^D(A)\ 

n — > oo (r„ — > oo). Then, there exists D 5 > 0, 



Suppose that g(— ^i—lr n e^ e ) diverges as n —> oo ^r„ — > oo;. men, mere exists u 5 
which has an upper bound determined by constants C\ , C 2 ,/3i and J3 2 , such that \F\(r n e 
is estimated as 



upper bound 
Therefore, we obtain 

e At F { (A)dA 



F^e^l < D 5 e^ + M r ». (4.25) 

shows that there exists D 6 > , which has an 
Zi.B\ and Bi. such that 



If g(- V-I^e^") is bounded, Eq. (l4.22l) shows that there exists 
or *>«"«4 determined by constants C\, C 2 ,/3\ and /3 2 , such that 



(4.26) 



c 4 ? 
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with some D 7 > 0. Thus if t > t := n(fi\ + fa)/ (26), this integral tends to zero as n — » oo, 
which proves Eq. (l4.18l) . From the above calculation, it turns out that the integrals along 
C2 to C(, have upper bounds determined by C\ , Ci,fi\ and fi 2 once g is fixed. This fact will 
be used to prove Thm.5.10. 

In particular when < K < K c , there are no eigenvalues on the right half plane 
(Thm.3.5). Thus Eq. (l4~T8l) is reduced to Eq. (l4T9l) . 

Finally, we prove that all resonance poles are on the left half plane when < K < K c . 
Resonance poles are roots of Eq. (l4.1 II) . which is reduced to Eq. (l3. 171) on the imaginary 
axis. We have already shown that K c is the least value such that Eq. (l3 . 171) has a root y. 
Therefore, when < K < K c , there are no resonance poles on the imaginary axis. ■ 

The condition K 7<"i is used to prove that the infinite series in Eq. (l4. 1 81) converges. 
However, to prove the exponential decay of the order parameter rj(t) = (e Tl '<p, P ), we do 
not need the infinite series such as Eq. (l4.18l) . In Sec. 5. 3, we will give another expression 
of the semigroup (e T[t (p, ip) consisting of a finite sum and a certain integral, which proves 
that 77(f) decays exponentially without the assumption K £"K\. Nevertheless, for almost 
all choices of g(co), "K\ = {0, 00} and thus K <$. "K\ is not a restriction. 

Remark. Because of the assumption (A3), there are no resonance poles of T\ on the area 

{z e C I |arg(z) - n/2\ < 6, |arg(z) - 3n/2\ < 6, \z\ > r], (4.28) 

if r > is sufficiently large. Further, the left hand side of Eq. (14. Ill) is meromorphic. This 
shows that on each vertical line {A | Re(A) = a < 0}, there are only finitely many resonance 
poles. Thus Re(A n ) — » -00 as n — » 00, and we can renumber the index of resonance poles 
A n so that 

> Re(,lo) > Re(,li) > Re(/l 2 ) > • • • . 

Note that Exp + (0) is the set of bounded holomorphic functions near the upper half 
plane. From the proof above, we immediately obtain the following. 

Corollary 4.5. Suppose that K £"K\ and g(co) satisfies the assumptions (Al) to (A3). If 
(p, if/ £ Exp + (0), then Eq. (l4.18l) is true for t > 0. 

4.3 Rational case 

In this subsection, we suppose that g(u>) is a rational function. Since g(co) does not decay 
so fast as \a>\ — > 00, we should choose moderate functions for (p and iff. Let C + = {z 6 
C I Im(z) > 0} be the real axis and the upper half plane. Let H + be the set of bounded 
holomorphic functions on C+. With the norm 

II0H := sup |0( Z )|, (4.29) 

Im(.-)>0 

H + is a Banach space. 

It is remarkable that if g(u>) is a rational function, Eq. (l4. Ill) is reduced to an algebraic 
equation. Thus the number of resonance poles is finite. The proof of the following the- 
orem is similar to that of Thm.4.4 and omitted here (since F (A) and Fi(A) are rational 
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functions, estimates of integrals in Eq. (l4.23l) are too easy. We need not introduce the set 
%). 

Theorem 4.6. Suppose that < K < K c and g(a>) is a rational function. For any 
<p, ifj e H + , the semigroup e Tit satisfies the equality 

M 

for t > 0, where Aq, ■ ■ • , A M are resonance poles of T\ and R [t, <p, \]/]e Aot , • • • ,R M \t,<p, i]/]e AM 
are corresponding residues of F 1 (A)e M . In particular, when < K < K c , Ao, •• • ,A M are 
on the left half plane and the order parameter rj(t) = (e T[t (f>, P ) for the linearized system 
(13.41) decays to zero exponentially as t — » oo. 

Since the right hand side of Eq. (14.30l) is a finite sum, the semigroup e Tit looks like an 
exponential of a matrix. The reason of this fact will be revealed in Sec. 5. 4 by means of 
the theory of rigged Hilbert spaces. 

Example 4.7. If g(oj) = l/(n(l + co 2 )) is the Lorentzian distribution, a resonance pole 
is given by A = K/2 - 1 (a root of Eq. (l4.1 II) ). Therefore rj(t) decays with the exponential 
rates e {K/2 ~ l)t . Note that the transition point is K c = 2/n/g(0) = 2. 



5 Spectral theory 

We have proved that when K > K c , the de- synchronous state (rj(t) = 0) is linearly unstable 
because of eigenvalues on the right half plane, while when < K < K c , it is linearly stable 
because of resonance poles on the left half plane. Next, we want to investigate bifurcations 
at K = K c . However, a center manifold in a usual sense is of infinite dimensional because 
the continuous spectrum lies on the imaginary axis. To handle this difficulty, we develop 
a spectral theory of resonance poles based on a rigged Hilbert space. 

5.1 Rigged Hilbert space 

Let X be a topological vector space and X' its dual space. X' is a set of continuous anti- 
linear functionals on X. We use the Dirac's notation; for jx 6 X' and <p 6 X, fj.(<p) is denoted 
by (<p I n)- For any a, b e C and <p, i// e X, \x, £ 6 X' , the equalities 

(a(f> + bif/\/j.) = a(<f> \fi) + b(\ff\fj), (5.1) 
(<P\an + b& = a(<f>\fi) + b(<t>\&, (5.2) 

hold. Several topologies can be defined on the dual space X' . Two of the most usual 
topologies are the weak dual topology and the strong dual topology Il46ll . A sequence 
{fij} c X' is said to be weakly convergent to /u e X' if (cp | fij) — » (<p | //) for each ^ el; 
a sequence {fj.j} c X' is said to be strongly convergent to /u e X' if (<p \ fij) — * (<p \ /j.) 
uniformly on any bounded subset A c X. 
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Let H be a Hilbert space with the inner product (• , -) H such that X is a dense subspace 
of H. Since a Hilbert space is isomorphic to its dual space, we obtain H c X' through 

H as H'. 

Definition 5.1. If a topological vector space X is a dense subspace of a Hilbert space H 
and a topology of X is stronger than that of H, the triplet 

XcHcX' (5.3) 

is called the rigged Hilbert space or the Gelfand triplet. The canonical inclusion i : X — » 
X' is defined as follows; for t/r e X, we denote z'(i/0 by which is defined to be 

*OA)(0) = <0l<A> = (0,iAk, (5.4) 

for any e X. If the inner product on H is Hermitian, the alternative definition 

iMty) = (cf>\~iP) = (<P,M H . (5.5) 

is used so that i(ai//)(4>) = ai{i]/){(p) for a e C. After section 5.2, we will use the latter one. 
Thus if H = L 2 (R,g(u))du)), then 

z'(i/O(0) = I (p{oj)^{oj)g{oj)dco. 
Jr 



Let A : X — » X be a continuous linear operator on X. Then, the dual operator A x : 
X' -> X' defined through 

(<f>\A x fi) = (Acf>\n), e X, (i e X' (5.6) 

is also continuous on X' for both of the weak dual topology and the strong dual topology. 
We can show the equality 

A x z(<A) = W), (5.7) 

for any \p e X. 

It is easy to show that the canonical inclusion is injective if and only if X is a dense 
subspace of H, and the canonical inclusion is continuous (for both of the weak dual topol- 
ogy and the strong dual topology) if and only if a topology of X is stronger than that of H 
(see Treves ll46ll ). If X is not dense in H, two functionals on H may not be distinguished 
as functionals on X. As a result, H' <t X' in general. 

Definition 5.2. When X c H is not a dense subspace of H, the triplet (X, H, X') is called 
the degenerate rigged Hilbert space. 

For applications to the Kuramoto model, we investigate two triplets, Exp + c L 2 (R, g(co)daj) c 
Exp^, and a degenerate one (H + , L 2 (R, g(co)du)), H' + ). 
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5.2 Spectral theory on Exp + c L 2 (R, g(a>)da>) c Exp' + 

In this subsection, we suppose that g(co) satisfies the assumptions (Al) to (A3). Since 
g decays faster than any exponential functions e~^°"\ we have Exp + c L 2 (R, g(co)da>), 
and indeed, Exp + is dense in L 2 (R, g{oj)doS) and the topology of Exp + is stronger than 
that of L 2 (R, g(aj)du)) (see Prop. 5. 3 below). Thus the rigged Hilbert space Exp + c 
L 2 (R, g(co)doj) c Exp+ is well defined. Recall that Exp + (J3, n) is the Banach space of 
holomorphic functions on C„ with the norm || • \\^ n , and Exp + (/3) is their inductive limit 
with respect to n > 1. By the definition of the inductive limit, the topology of Exp + (/3) is 
defined as follows: a set U c Exp + {fi) is open if and only if U D Exp + (J3, n) is open for 
every n > 1. This implies that the inclusions Exp + (/3, n) — » Exp + (/?) are continuous for 
every n > 1. Similarly, Exp + is an inductive limit of Exp + (J3), and its topology is induced 
from that of Exp + (fi): a set U c Exp + is open if and only if U n Exp + (J3) is open for every 
f3 = 0, 1, • • • . The inclusions Exp + (J3) — » Exp + are continuous for every = 0, 1, • • • . On 
the dual space Exp' + , both of the weak dual topology and the strong dual topology can be 
introduced as was explained (see also Fig. [9]). 

A topological vector space X is called Montel if every bounded set of X is relatively 
compact. A Montel space has a convenient property that on a bounded set A of a dual 
space of a Montel space, the weak dual topology coincides with the strong dual topology. 
In particular, a weakly convergent series in a dual of a Montel space also converges with 
respect to the strong dual topology (see Treves [461). This property is very important to 
develop a theory of generalized functions. 

The topology of Exp + has following properties. 

Proposition 5.3. Exp + is a topological vector space satisfying 

(i) Exp + is a complete Montel space. 

(ii) if {(f>j}J =l is a convergent series in Exp + , there exist n > 1 and /3 > such that {(f>j}J =l c 
Exp + (yS, n) and {(f>j}J =1 converges with respect to the norm || • H^,,. 

(iii) Exp + is a dense subspace of L 2 (R, g(oj)dco). 

(iv) the topology of Exp + is stronger than that of L 2 (R, g(uj)dco). 

Proof, (i) At first, we prove that Exp + (J3) is Montel. To do so, it is sufficient to show 
that the inclusion Exp + (J3, n) — » Exp + (fi, n + 1) is a compact operator for every n (see 
Grothendieck IBH . Chap.4.3.3). To prove it, let A be a bounded set of Exp + {fi, n). There 
exists a constant C such that = sup_ eC( e~^\(p{z)\ < C for any cp e A. This means that 
the set A is locally bounded in the interior of C„. Therefore, for any sequence {<f>j}™ =l C 
A, there exists a subsequence converging to some holomorphic function if/ uniformly on 
compact subsets in C„ (Montel's theorem). In particular, the subsequence converges to if/ 
on C n+ \, and it satisfies ||0l|y3,„ + i < C and if/ e Exp + (/3, n + 1). This proves that the inclusion 
Exp + (/3, n) — » Exp + (/3, n + 1) is compact and thus Exp + (/?) is Montel. In a similar manner, 
we can prove by using the Montel's theorem that the inclusion Exp + (J3) — > Exp + (/3 + 1) is 
a compact operator for every ft = 0, 1, • • • , which proves that Exp + is also Montel. Next, 
we show that Exp + is complete. Since Exp + (J3, n) is a Banach space, in particular it is a 
DF space, their inductive limit Exp + (/?) is a DF space by virtue of Prop. 5 in Chap.4.3.3 of 
BOl . in which it is shown that an inductive limit of DF spaces is DF. The same proposition 
also shows that the inductive limit Exp + of DF spaces Exp + (/?) is a DF space. Since Exp + 
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is Montel and DF, it is complete because of Cor.2 in Chap.4.3.3 of 

(ii) It is known that if the inclusion Exp + (/3) — » Exp + (/3 + 1) is a compact operator 
for every fi = 0, 1, • • • , then, for any bounded set A c Exp + , there exists /S such that 
A c Exp + (/?) and A is bounded on Exp + (/3) (see Komatsu If23l and references therein). By 
using the same fact again, it turns out that for any bounded set A c Exp + , there exist f3 
and n such that A c Exp + (/?, n). In particular, since a convergent series {(pj}°° =l is bounded, 
there exists ft and n such that {(f>j}J =l C Exp + (/3, n) and it converges with respect to the 
topology of Exp + (/J, n). 

To prove (iii), note that L 2 (R, g(aj)du>) is obtained by the completion of the set of 
polynomials because the probability density g(u>) has all moments due to the assumption 
(A2). Obviously Exp + includes all polynomials, and thus Exp + is dense in L 2 (R, g(to)dco). 

For (iv), note that Exp + satisfies the first axiom of countability because it is defined 
through the inductive limits of Banach spaces. Therefore, to prove (iv), it is sufficient to 
show that the inclusion Exp + — » L 2 (R, g{co) doS) is sequentially continuous. Let {<pj}J =1 be 
a sequence in Exp + which converges to zero. By (ii), there exist /3 and n such that {4>j\°j = \ 
converges in the topology of Exp + (/3, n): ||0/||/},„ — » 0. Then, 

Xco 
\<f>j(uj)\ 2 g(aj)du) 
CO 

XCO 
e ma)] g(co)d(o 
oo 

XCO 
e^g(co)daj. 
CO 

By the assumption (A2) for g(oj), the right hand side exists and tends to zero as j — > oo. 
This means that the inclusion Exp + — » L 2 (R, g(co)doj) is continuous. ■ 

The topology of the dual space Exp^ has following properties. 
Proposition 5.4. 

(i) Exp^ is a complete Montel space with respect to the strong dual topology. 

(ii) Exp' + is sequentially complete with respect to the weak dual topology; that is, for a 
sequence {fij}°° =l , if (0 1 fij) converges to some complex number for every (p e Exp + as 
j — » oo, then there exists /i e Exp' + such that = (tp\ ft) and — » /i with respect to the 
strong dual topology. 

Proof, (i) It is known that the strong dual of a Montel space is Montel and complete, see 
Treves [|46l . (ii) Suppose that converges to some complex number for every 

cp e Exp + . This means that the set j is weakly bounded and is a Cauchy sequence 
with respect to the weak dual topology. As was explained before, on a bounded set of a 
dual space of a Montel space, the weak dual topology and the strong dual topology coin- 
cide with one another. Thus {Hj}°° =l is a Cauchy sequence with respect to the strong dual 
topology. Since Exp^ is complete with respect to the strong dual topology, fij converges 
to some element /u e Exp' + . In particular, (0 1 fij) converges to ((p\/u) = C^. ■ 

Next, we restrict the domain of the operator T\ = V-lAI + f P to Exp + . We simply 
denote 7"i |ex P+ by T\. We will see that T\ is quite moderate if restricted to Exp + . 
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Proposition 5.5. 

(i) The operator T\ : Exp + — » Exp + is continuous ( note that it is not continuous on 
L 2 (R,g(to)doo)). 

(ii) The operator T\ : Exp + — » Exp + generates a holomorphic semigroup e Tlt : Exp + — > 
Exp + on the positive t axis (note that it is not holomorphic on L 2 (R, g(to)dto)). 

(iii) Eigenvalues of Ti : Exp + — > Exp + are the same as those of T\ : L 2 (R, g(to)doo) — » 
L 2 (R, g(a>)<ia>) (i.e. the roots of D(A) = on Re(/l) > 0). The residual spectrum of 
Ti : Exp + — » Exp + is the imaginary axis and the left half plane. There is no continuous 
spectrum of T\ : Exp + — > Exp + . 

Proof, (i) It is easy to see by the definition that if e Exp + , then e Exp + . Let {<f>j}J =l 
be a sequence in Exp + converging to zero. By Prop. 5. 3 (ii), there exist f3 > and n > 1 
such that H0/II/?,,, — » 0. For any e > 0, ||7 , 1 y -||y 3+e ,„ is estimated as 

H^10;ll/? +e , ;1 < II V^T^-II^ + -\(<f>j, P )\ ■ \\P \\p +£ ,n 



< sup e-^ +e)M \co^j(co)\ + ^||0y|| L 2 (Rjg(w 
weC„ z 



< 110/ll/j,,, • sup e' £M \oj\ + T ll0yllL2 ( R, g M^), 

which tends to zero as j — » oo. This proves that tends to zero as j —* oo with respect 
to the topology of Exp + , and thus T\ : Exp + — » Exp + is continuous. 

(ii) We know that the operator 7^ generates the semigroup e T[t as an operator on L 2 (R, g(to)dto) 
(see Sec. 4.1). In other words, the differential equation 

d i — K 

—x(t, to) = T lX (t, to) = y-l(oxt(t, to) + -(x(t, •), Po) (5.8) 
dt 2 

has a unique solution x(t, oo) = e T[t (p(to) in L 2 (R, g(to)dto) if an initial condition (p is in 
L 2 (R, g(to)dto). We have to prove that if </> 6 Exp + , then x(t, ■) e Exp + . For this purpose, 
we integrate the above equation as 

e Tlt ftto) = e^^ftto) + ^ f e" rioj(t ~ s \e TlS (f>,Po)ds. (5.9) 

2 Jo 

From this expression, it is obvious that if 6 Exp + (/3, n), then e Tlt (p(to) e Exp + . Next, 
we show that e T[t is a holomorphic semigroup. The proof is done in the same way as 
the standard proof of the existence of holomorphic semigroups |l22l) . Let 6 be a positive 
number appeared in the assumption (A3). In the proof of Thm.4.4, we see that the Laplace 
integral along C 2 ,C^,Cs and C(, tends to zero as r n — » oo (see FigfT] (b) and Eq. (l4.23l) ). 
This means that the contour (the straight line represented as Figj7] (a)) for the Laplace 
inversion formula (14.41) can be deformed to the contour Y represented in FigJH which 
goes to infinity along the rays axg(A) = + (n/2 + 6): 

(e^frtf,) = — l — [e^X-T^^dX, (5.10) 
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Fig. 8: The contour for the Laplace inversion formula. 



for t > 0. We show that this integral exists even if f is a complex number. For this purpose, 
put A = re ± ^~ I(6+ ' T/2) e Y and t = \t\e^ with \6\ < S - e, where e > is an arbitrarily 
small number. Then 

_ | g r|(|cos(6»±(5+ff/2))i < g -r|f|sine 

Further, ((A - T{)~ l ft,ft) decays as 0(1 /\A\) as \A\ — » oo along T because of Eq. (l4.5l) . 
the assumption (A3) and Claim in the proof of Thm.4.4. This proves that the integral 
in Eq. (l5. 101) exists for t = \t\e^ e , \6\ < 8 - e, which gives an analytic continuation of 
(e Tl 'ft, ft) to the sector {t e C | arg(?) < 6}. In particular, (e Tlt ft, P ) is holomorphic in the 
sector. Now Eq. (|5.9b proves that e Tl 'ft is also holomorphic and included in Exp + in the 
sector {t e C \ arg(t) < 6}. 

(iii) Eigenvalues are calculated in the same way as Prop. 3. 2 (i). Note that corresponding 
eigenvectors 1 /(A - y-lco) are in Exp + when Re(/l) > 0. Next, from the proof of Lemma 
4.1, it turns out that the resolvent is given by 

i / k i r°° i \ 

(A-T 1 r 1 <f>= — \ft(io) + -- — — — — ft(co)g(oj)dco , (5.11) 

for ft e Exp + . Suppose that Re(A) < 0. Since 1 /(A - yf^Aco) £ Exp + , (A - T x y x ft e Exp + 
if and only if ft satisfies 

K 1 f°° 1 r— 

0(w) + ^——z ^^ft(oj)g(aj)daj = (A- i-\co)ft(co) 

for some ft e Exp + . In particular, ft has to satisfy 

(k \ k r°° i 
1_ W) U(_V=U) = -- —ft(co)g(<o)dco. (5.12) 
^ i 1 J-ooA - y-loj 
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For fixed A, define continuous linear functionals F U F 2 on Exp + to be 

Fd<f>] = (i-|dw)0(-V^Ta), 
k r°° 1 

F2W = -t -=-^juS)g{o))do), 

^ J-00 a — y—ltt) 

respectively. If F\ = F 2 on a dense subset of Exp + , then F\ = F 2 for any (p e Exp + 
because F\ and F 2 are continuous. However, it is obvious that Fi[l] ± ^[l]. This proves 
that the set of functions (p satisfying (15.121) is not dense in Exp + , and thus the domain of 
(A - T\)~ l is not dense in Exp + when Re(A) < 0. Therefore, the imaginary axis and the 
left half plane are included in the residual spectrum. Finally, suppose Re(A) > and A is 
not an eigenvalue. Then, (A - Ti)~ l (p e Exp + for any (p <= Exp + . From Eq. (l5.1 II) . it is easy 
to prove that if a sequence {(/>j}J =l tends to zero in Exp + , then \\(A - T\)~ l (pj\\^ n tends to 
zero for some /? and n, which proves that (A - T\Y { <pj — » in Exp + . Therefore, (A - T{)~ 1 
is continuous and T\ has no continuous spectrum on the right half plane. ■ 

Eigenvalues of T\ are given as roots of the equation D(A) = 2/K, Re(A) > 0, and 
corresponding eigenvectors are 

v A {co) = e Exp + . (5.13) 

A - y-\co 

If we regard it as a functional on Exp + through the canonical inclusion i : Exp + — > Exp' + , 
it acts on Exp + as 

Kvd(<f>) = {<P\va) = ((f>,v A ) = f (p{oj)v A {(jj)g{oj)dto = f 1 — ^){w)g{io)dw, (5.14) 

Jr JrA - y-lco 



for 4> e Exp + . Due to Eq. (l4.7l) . the analytic continuation of this value from right to left is 
given as 

^-^(f>((D)g((o)d(o + 2tt0(- V-U)g(- V-U). (5.15) 

IrA - V-lcd 

Motivated by this observation, let us define a linear functional n(A) e Exp' + to be 



/- 

Jr. 



(4> I fx(A)) = f L_0( w )g( w )rf w + 2;r0(- V^lA)g(- V^U), (5.16) 

JrA - y-lco 

when Re(/l) < 0, and 

<^|//(A)>= lim f — ^gic^doj, (5.17) 

•^+° Jr(jc + yhly) - y^ioj 

when A = yf-Ty, y e R. Indeed, /x(A) is continuous and thus an element of Exp' + . Let 
{<f> m }°° =1 be a sequence in Exp + converging to zero. By the definition of the topology on 
Exp + , <f> m (oS) converges to zero uniformly on compact sets on the upper half plane. Hence, 



36 



the second term <p(- V-lA)g(— yf-lA) of Eq. (l5. 161) tends to zero as m — » oo. The integral 
of the first term of Eq. (l5. 161) is decomposed into three terms J M ■ ■ ■ dco, JTJ • • -dco and 

X— M rM 
m • • • d<x>. Since m (a») — > uniformly on compact sets, the term J_ M • • • dco tends to 

J/"*oo r—M 
•• • dco exists, two terms J • • • and J _ m • • • dco become 

sufficiently small if M is sufficiently large. This proves that (<p m \fi(A)) —* as m —* oo. 

Similarly, we can verify for the case A = V-Ty that jd(A) is a continuous functional (note 

that the existence of the limit in Eq. (l5. 171) was shown in the proof of Lemma 4.2). 

We expect that /u(A) plays a similar role to the eigenvectors. Indeed, we can prove the 

following theorem. 

Theorem 5.6. Let A ,A U --- be the resonance poles of the operator T y and T* the dual 
operator of T Y defined through 

(cf>\ 7?f> = (T 1( f>\£), <f> e Exp + , f 6 Exp;. (5.18) 

Then, the equality 

T* f i(A n ) = A nM (A n ) (5.19) 

holds for n = 0, 1, 2, • • • . In this sense, A n is an eigenvalue of Tf, and fi(A n ) is an eigenvec- 
tor. In what follows, /x(A n ) is denoted by \x n and we call it the generalized eigenfunction 
associated with the resonance pole A„. 

Proof. The proof is straightforward. Suppose that Re(A n ) < 0. For any cf> e Exp + , 
(<P\T^ n ) = (T l( f>\ia n ) 

= f l -^(T 1( f>)(aj)g(co)daj + 2n(T 1( f>)(- ^lA n )g(- V^LL.) 

JrA h - y-lco 

= f — (p(to)g(to)dto + ^ f 1 — g(co)dco ■ f (p{oj)g{co)dco 

JrA h - y-lco 1 JrA u - y-lco Jr 

+2tt|^0(- V-UJ + | jr^(fi>)g(w)dwjg(- V-UJ 

= A J f l -^(f>(Lo)g(to)dto + 2ncp(- <-iA n )g{- V^H.) ) 

\JrA„ - y-lco I 

+| fj(co)g(to)dco\D(A n ) + 2ng(- V^UJ - j\ 

= A^\fi„) + ^ Jj>(co)g(co)dw{D(A n ) + 2ng(-yf^A n )-jy (5.20) 

Since A n is a resonance pole, it is a root of Eq. (l4. 111) . Thus we obtain 

(<p | T*/j, n ) = A n ((p | ju„) = (<p | A n fi n ), 
which proves the theorem. The proof for the case Re(A n ) = is done in the same way. ■ 
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Define a dual semigroup (e Tit ) x through 

(4>\(e Tit )y) = (e Tit 4>\fi). (5.21) 

for any <f> e Exp + and /u e Exp' + . 

Proposition 5.7. (i) A solution of the initial value problem 

j t Z=Tl& ^(0)=^gEx P ; (5.22) 

in Exp' + , which is continuous in an initial condition with respect to the weak and the strong 
dual topologies, is uniquely given by g(t) = (e Tlt ) x fi. 

(ii) (e Tlt ) x has eigenvalues e Aot , e A[t , • • • , where A , A u ■ ■ • are resonance poles of T x . 
Proof, (i) For any <f> e Exp + , (e T[t ) x /j. satisfies 

^<0|(e r ")V) = ^-(e Tl, 4>\fi) = (T l e Tit cf>\fi) 
at at 

= (e Tl, T l c(>\ f i) = ((t>\T x (e Tl, ) x M). 

This provides 

4(e r '0V = T x (e Tl ') x iu. (5.23) 
at 

It is easy to verify from Eq. (l5.21l) that when t = 0, (e T,t ) x is the identity. This proves that 
(e Ti ') x {i is a solution of the initial value problem (15.221) . 

To prove the uniqueness of a solution, let us show that a solution of the equation 
dcp/dt = T\(p in Exp + is unique (this fact was mentioned in the proof of Thm.5.5 (ii)). 
Indeed, let <p\ and 2 be two solutions of the equation with the same initial condition in 
Exp + . Then, Eq. (l5.9l) provides 

Mt, oj) - <h(t, co) = - J e^^iMs, ■ ) - <p 2 (s, ■ ), P )ds. 



This is estimated in L 2 (R, g(u))dco) as 



\\<Pl(t, ■ ) - fc(U ■ )\\LHR,g(<o)da>) ^ ^ I K<M S ' " ) _ 020, • ), Po)\ds 



< 
~ 2 



K 

k r*. 



his, - )-<p2(s, -^LHR^d^ds. 



Now the Gronwall inequality proves that \\(p\(t, •) - (f> 2 (t, • )llL 2 (R,g(w)tM = 0, and thus a 
solution is unique in L 2 (R, g(a>)dto). Since Exp + c L 2 (R, g(a>)du)), a solution is unique in 
Exp + . For a solution <p(t, to) e Exp + , the canonical inclusion yields 



d 
dt 



-K<t>) = m<p) = Tfm. 
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Since i : Exp + —* Exp' + is injective, it turns out that a solution of (15.221) in z'(Exp + ) is 
unique. Now we suppose that the initial value problem (15.221) has two solutions f\(t,fi) 
and f 2 (t,fj), which are continuous in \x e Exp' + . There exists a sequence {£ n }™=i c *(Exp + ) 
such that £ n — » n with respect to the weak dual topology and the strong dual topology 
because z'(Exp + ) is dense in Exp' + . Since a solution in z'(Exp + ) is unique, f\(t, g n ) = f 2 (t, £„). 
Taking the limit yields fi(t,fi) = fiit,[i), and they coincide with (e T,t ) x /j.. 
(ii) Eq. (l5.19D implies 



d 
dt 



Thus = e Ant /u n is a solution of the equation d^/dt = Tfg. Obviously e Ant fj,„ is continu- 
ous in fi n e Exp' + . By the uniqueness of a solution, we obtain e /! "'//„ = (e Tlt ) x fj. n . ■ 

If we define a semigroup e 7 *' generated by T* to be the flow of (15.221) . then Prop. 5.7 
(i) means e T *' = (e T,t ) x . Prop. 5.7 (i) also implies that a solution of the inhomogeneous 
equation 

j t £ = T^ + f(t), f(t)eExp' + , (5.24) 
which is continuous in an initial condition /i 6 Exp+, is uniquely given by 

= {e Tl 'Yii + f (e Tl{t ~ s) ) x f(s)ds. (5.25) 
Jo 

This formula will be used so often when analyzing the nonlinear system (I3.2I) . (I3.3I) . 

In what follows, we suppose that < K < K c and K TCi so that Eq. (l4.19l) is 
applicable. To develop a spectral decomposition, we make the additional assumption. 

(A4) All resonance poles are poles of first order (i.e. all resonance poles are simple roots 
ofEq.<gT[])). 

Note that this is the assumption only for g(a>). For example, the Gaussian distribution 
satisfies (A4), see Appendix A for the proof. Under the assumption (A4), R n [t, (p, ip] in 
Eq. (|4. 191 ) is independent of t. Since R n [t, </>, ip] is the residue of Fi(A) given as Eq. (l4.5l) . it 
is calculated as 

K 



where 



R„[t,<f>,lff] = —((/>\Mn>(i/f\fin), (5-26) 
LL) n 

D n := lim — !— f 1 - ^D(A) - nKg(- V^U)) (5.27) 

is a constant which is independent of (p, ip. Note that Q[A, (f>] given by Eq. (l4.6l) is just the 
definition of the functional fx(A). Thus Eq. (l4. 191) is rewritten as 

oo ^ 

(.e Tit (f), </>) = V ^-e x "\<t> | fi n )(if, | iu n ), (5.28) 

for t > to. By using the canonical inclusion i, the left hand side is rewritten as 

(e Tlt cf>,^) = (e T ^\t) = (<t>\(e T,t ) x ^). (5.29) 
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Therefore, we obtain 

(e T ") x m = ^<*>ty \ Mn ).\ n n ), (5.30) 

,j=0 ' Lu n 

for t > to- Note that the complex conjugates arise in the coefficients because functionals 
are anti-linear. Since Eq. (l5.30l) comes from Eq. (l5.28l) . the infinite series in the right hand 
side of Eq. (l5.30l) converges with respect to the weak dual topology on Exp^. However, 
since Exp + is Montel, it also converges with respect to the strong dual topology. We divide 
the infinite sum in Eq. (l5.28l) into two parts as 

M K °° K 

(e Tit <f>,i{,) = Y j —e A " t (<t>\iu n )(i{,\iu n )+ J] (5.31) 

n=0 " n=M+\ " 

where M e N is any natural number. The second part zZ7=m+i ["'] does not converge when 
< t < t Q in general. However, since (e Tlt <p, if/) is holomorphic in t > and continuous at 
t = 0, we obtain 

M K K 
(0,<A) = X — -<0|>u„><tA|ju ;i > + lim V — -e A »'((P\n n Xi{,\n n ), (5.32) 

where the second term has a meaning in the sense of an analytic continuation in t. Through 
the canonical inclusion, the above yields 

m K 

\lff) = V— <^|/i„>-|i"n>+^M[«AL (5.33) 

00 k - 

W]:=lim X — e^\fi n )-\fi n ), 

,j=M+l ' Lu n 

which gives the spectral decomposition of | 0-) e ?'(Exp + ). 

Theorem 5.8 (Spectral decomposition). Suppose that < K < K c , K <£. K.\ and g(oS) 
satisfies the assumptions (Al) to (A4). 

(i) A system of generalized eigenfunctions {[i n }™ =0 is complete in the sense that if {(p \ yu„) = 
Ofor n = 0, 1, • • •, then = 0. 

(ii) fj.o, fj,i, ■ ■ ■ are linearly independent of each other: if 2,?=o a n \/u n ) = with a n e C, then 
a n = for every n. 

(iii) Let V M be a complementary subspace of span{/i , ■ ■ ■ in Exp^, which includes 
ix j for every j = M + 1 , M + 2, • ■ ■ . Then, any | iff) e z'(Exp + ) is uniquely decomposed with 
respect to the direct sum span{/i , • • • ,/u M } © V M as Eq. (l5.33l) . and this decomposition is 
independent of the choice of the complementary subspace Vm- 

Proof, (i) If <0 1 //„> = for all n, Eq. d5T28l) provides (e r 'V,^) = for any iff e Exp + . 
Since Exp + is dense in L 2 (R, g(a))dco), we obtain e Tl '(p = for t > t . Since e T{t is 
holomorphic in t > and strongly continuous at t = 0, we obtain (f> = by taking the limit 
t +0. 
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(ii) Suppose that £,?=o a n \fj. n ) = 0. Since (e Tl ') x is continuous, 

CO CO CO 

= (.e Tlt ) x J]a n \Hn) = Y, a n(e Txt T\Hn) = J]a n e 1 » t \fi n ). 

n=0 n=0 n=0 

We can assume that 

> Re(i ) > Re(/li) > Re(i 2 ) > ■ ■ ■ , 

without loss of generality. Further, on each vertical line {/l|Re(/l) = a < 0}, there are 
only finitely many resonance poles (see Remark in Sec.4.2). Suppose that Re(/io) = • • • = 
Re(/l,t) and Re(Ak) > Re(Ak+i). Then, the above equality provides 

k oo 
n=0 n=k+\ 

Taking the limit t — » oo yields 

k 

= YimYa n e-^ m \Hn). 

f— >oo i—l 

n=0 

Since the finite set ,f*k of eigenvectors are linearly independent as in a finite- 

dimensional case, we obtain a n = for n = 0, • • • , k. The same procedure is repeated 
to prove a„ = for every n. 

(iii) Let V M and V' M be two complementary subspaces of span{/i , ■ ■ ■ ,Hm\, both of 
which include (Xm+uMm+i, •"• Let 

Exp' + = span{ju , ■ • • ,Hm) © V M = span{ju , ■ • • ,Mm) © V' M 
be two direct sums and let 

M M 

1 <A> = ^ a "l^ + V = J] a 'n\Vn) + V', V 6 V M , V 6 V' M 

n=0 n=0 

be corresponding decompositions. Then, (e T[t ) x \ if/) is given by 

M M 

(e Tit r\^) = ^a„^V„> + (e Tl O x v = J] a' n e J "'\ Li n ) + (e r ")V. 

They give decompositions of (e Tlt ) x \ iff) with respect to two direct sums 

span{(e r ") x //o, ■ • ■ , {e Tlt ) x Lx n ) © (e r ") x V M , span{(e r '0Vo, ■ • ■ , (e Tlt ) x fx n } © (e r ") x V*. 

respectively. Since (e Tlt ) x \fX n ) = e* nt \fx n ), span{(e r ' f ) x yUo, ■ • • , (e Tlt ) x ju n ) = span{// , ■ ■ ■ ,Hn), 
and the sets (e Tlt ) x VM and (e T[t ) x V' M also include [Xm+i,Hm+2, ■ * • ■ Because of part (ii) of 
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the theorem, the decomposition of (e Ti ') x \ iff) with respect to above direct sums is uniquely 
given as Eq. (15.30l) for t > t . This implies 

a n = a' n (= \ Mn )), (e Tlt ) x v = (e r ")V (= £ —^\ lff \ /in ) . \ Mn )). 

2D„ n =M+\ 2D n 

Since (e Tl ') x is continuous in t > 0, we obtain v = V . ■ 

When K > K c , the spectral decomposition involves eigenvalues and eigenvectors in 
Exp + . In particular, if all eigenvalues are not degenerate (that is, all roots of Eq. (l3.12l) are 
single roots), then Eq. (l4. 1 81) proves 

K _ _ » K 



(e T «YW) = X -W»W„> -\v„) + V — e*.^|^> • |^>, (5.34) 
for t > t , and 



n=0 „ = o 2D„ 



M ^ M' ^ 

I <A> = X — <^ I v„> • I v„> + V — (iff \ fin) ■ I/O + n M >m, (5.35) 

,i=0 n=0 

where £ , •• • ,% M are eigenvalues of T y on the right half plane, v„ := v^ n (a») e Exp + are 
corresponding eigenvectors defined by Eq. (l5.13l) . and where E n are defined to be 

E n = lim — !— (l - f £>(£)) . (5.36) 

It is easy to extend Thm.5.8 to the case that K > K c : a system of eigenvectors {v„}^l 
and generalized eigenfunctions {/i n }^ is complete, and they are linearly independent. 
If eigenvalues and resonance poles are degenerate, the right hand sides of Eq. (l5.33l) and 
(15.351) become more complex involving polynomials in t. In this paper, we do not consider 
such cases for simplicity, though it will be treated in the same way as a finite dimensional 
case because multiplicities of eigenvalues and resonance poles are finite. 

This theorem suggests the expression of the projection to the generalized eigenspace. 

Definition 5.9. Denote by n„ : Exp' + — > span{//„} c Exp+ the projection to the general- 
ized eigenspace. For | iff) e ?'(Exp + ), it is given as 



n n \^) = —(iff\/i n )-\fi n ). (5.37) 
2D n 

Unfortunately, the projection IL, is not a continuous operator . For example, put 
iff m (co) = e ~^ mcj . Then, i(iff m ) = | e^ moj ) converges to zero as m — > +00 with respect to 
the weak dual topology by virtue of the Riemann-Lebesgue lemma. However, 

{i]j m \Hn) = f ] ~^ e ~ ^^gi^dco + 2ne-' nA "g(- V=H.) (5.38) 

JrA„ - V-lo> 
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does not tend to zero. It diverges as m — » oo when Re(/l„) < 0. This means that II„ : 
Exp^ — * Exp' + is not continuous with respect to the weak dual topology. To avoid such a 
difficulty caused by "weakness" of the topology of the domain, we restrict the domain of 
II„. For constants C > 1 and a > 0, define a subset V c ,a C Exp + to be 



V c , a = {<t> 6 Exp + | |0(z)|e~ a|z| < C when Im(z) > 0}. (5.39) 

We restrict the domain of the projection to i(V c , a ) C Exp' + and denote it by Tl^'" = ry^j. 
Note that the set {e~'^ ma) }^_ above is not included in Vc, a - To discuss the continuity of 
II„' a , let us introduce the projective topology on Exp' + (see also Fig.|9]and Table 1). In the 
dual space Exp^, the weak dual topology and the strong dual topology are defined. An- 
other topology called the projective topology is defined as follows: Recall that Exp + (fl, ri) 
is a Banach space, and its strong dual Exp + (/?, ri)' is a Banach space with the norm 

\H*p n := sup |<0|^>|. (5.40) 

We introduce the projective topology on Exp + (J3)' = P|«>i Exp + (/3, ri)' as follows: U c 
Exp + (/?)' is open if and only if there exist open sets U n c Exp + (J3, ri)' such that U n n 
Exp + (/?)' = U for every n > 1. It is known that the projective topology is equivalent to 
that induced by the metric 



^ 1 ,^ 2 ):=^- 1 + || ^_^, (5.41) 



see Gelfand and Shilov IfTTTl . When the projective topology is equipped, Exp + (J3)' is called 
the projective limit of Exp + (J3, ri)' and denoted by Exp + (/3)' = lim Exp + (jS, n)'. In a similar 
manner, the projective topology on Exp' + = P|/?>o Exp + (^)' is introduced so that U c Exp' + 
is open if and only if there exist open sets Up c Exp + (J3)' such that Up D Exp^ = £/ for 
every >0. This topology coincides with the topology induced by the metric 

2? 1 +dp(ju l ,iu 2 ) 

In this manner, Exp' + equipped with the projective topology is a complete metric vector 
space. When the projective topology is equipped, Exp + is called the projective limit of 
Exp + (J3)' and denoted by Exp' + = limExp + (^)'. 

By the definition, the projective topology on Exp' + is weaker than the strong dual 
topology and stronger than the weak dual topology. Since Exp + is a Montel space, the 
weak dual topology coincides with the strong dual topology on any bounded set of Exp^.. 
This implies that the projective topology also coincides with the weak dual topology and 
the strong dual topology on any bounded set of Exp' + . In particular, a weakly convergent 
series in Exp^ also converges with respect to the metric d and the strong dual topology. 
Note that the projective topology is stronger than that induced by the norm || • \\*^ n . There- 
fore, a convergent series with respect to the metric d also converges with respect to the 
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Exp + (0) c - c Exp + (i3) c ... c Exp + c Z 2 (R, g(co)da>) c Exp; c - c Exp + (j3)'c - c Exp + (0)' 



u u u n 



U canonical n 

Exp + (i3, «) inclusion Exp + (j3, «)' 

u n 

Exp + U 03, 1) Exp + n (^, 1)' 



Fig. 9: A diagram for the rigged Hilbert space Exp + c L 2 (R, g(co)doj) c Exp^. 



Exp + (/3, n) 


Banach space: /3,,, = sup \4>{z)\e /3|z| 

Im( z )>-l/n 


Exp + (/3, n)' 


Banach space: = sup |<0|£>| 

ll<%,.=i 


Exp + (£)' = limExp + (je,/i)' 




Exp+ = lim Exp + (/3)' 





Table 1: Metric vector spaces used in Section 5. 



norm || • for every /? and n. 

Theorem 5.10. For any C > 1 and or > 0, the following holds. 

(i) On i{V c ,a) c Exp' + , the weak dual topology, the projective topology and the strong dual 
topology coincide with one another. 

(ii) The closure i(Vc, a ) of i(Yc,a) is a connected, compact metric space. 

(iii) When K $ 7("i, ifj'" : i{V c , a ) —* Exp^ is continuous for every j = 0, 1, • • • with 
respect to the projective topology and sequentially continuous with respect to the weak 
dual topology. 

(iv) For the system (12.11) . give an initial condition po(0, co) = h(6), where h is an arbitrary 
measure on S l . Then, corresponding solutions of (13.21) . (13.31) satisfy Zj(t, ■ ) e Vi,n for any 
t>0andj = 1,2, 

Proof, (i) At first, we show that the set i(V c , a ) c Exp^ is equicontinuous. For any small 
s > 0, we define a neighborhood U = U(e) c Exp + of the origin so that if cf) e U nExp + (y), 
then 

sup <T rkl |0(z)| < eD(a,y), 

lm{z)>0 

where D{a, y) is a positive number to be determined. Then, for any (f> e U n Exp + (y) and 
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<A> 6 i(V C ,a), 



K0|iA)| < f \(p(io)\ ■ \if/(co)\g(io)duj 
Jr 

Jr 

< £CD(a,y) f e (a(+y)M g(w)dw. 
Jr 

Since g(w) decays faster than exponential functions, the integral E(a, y) : = J R e (a+y)M g(to)dco 
exists. If we put D(a, y) = \/CE{a,y), we obtain \(<p\ip)\ < e for any <p e U and 
| ip) e i(V c ,a). This proves that i(V c , a ) is an equicontinuous set. In particular, i{Vc, a ) is a 
bounded set of Exp' + for both of the weak dual topology and the strong dual topology (see 
Prop. 32. 5 of Treves II4610 . Since Exp + is Montel, the weak dual topology, the projective 
topology and the strong dual topology coincide on the bounded set i(Vc, a )- Thus it is 
sufficient to prove (ii) for one of these topologies. 

(ii) Obviously V c , a c Exp + is connected (actually it is a convex set). Since the canon- 
ical inclusion i is continuous, i(Vc, a ) and i(Vc, a ) are connected. Since the strong dual 
Exp^ is Montel (Prop. 5. 4), every bounded set of Exp' + is relatively compact, which proves 
that i(Yc,a) is compact. By the projective topology, i(Vc, a ) is a metrizable space with the 
metric (15 .421) . 

(iii) We show that Hf, a is continuous at e i{V c , a ) with respect to the projective 
topology. Then, U„' a proves to be continuous on i{V c , a ) because H„' a is linear. Since 
the projective topology is metrizable, it is sufficient to show that Tl^" is sequentially 
continuous. Since every weakly convergent series in Exp+ also converges with respect to 
the projective topology, it is sufficient to show that Yl„'° is sequentially continuous with 
respect to the weak dual topology. Take a sequence {iA m }^ =1 c V c , a such that i(ifr m ) = 
| ijj m ) — » € i(Vc,a) as m — » oo with respect to the projective topology. Since the weak 
dual topology is weaker than the projective topology, it also converges with respect to 
the weak dual topology. In the proof of Thm.4.4, we have obtained Eq. (l4.23l) for each 

e Exp + . Putting ifr = \j/ m in Eq. (l4.23l) yields 

M „ #(«) „ 

(e Tlt <f>,i/, m ) = V 7T£re ijt (<P I VjX^m I Vj) + lim V ^^'(t I f^jX^m I Hj) 

3— lim \ e A, F {A)dA 1 —= lim \ e A, F x {A)dA. 

In^Pl 11 ^ Jc 2 +c 6 In ^hl Jc i+ c 4 +c 5 

In the proof of Thm.4.4, we proved that J^ 2+c e At F (A)dA and +C4+C e M F\{A)dA tend to 
zero as n — » oo. In particular, we showed that when \<p(z)\ < C\e^ and |^(z)| < C2^ 2 ' z| 
on the contour, these integrals have upper bounds determined by C\, C2,fii and /?2- Thus 
if tff m e Vc, a , these integrals have upper bounds which is independent of m, and it is easy 
to verify that they tend to zero uniformly in m. As a result, the infinite series in Eq. (l4.18l) 
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converges uniformly in m, and we obtain 



M K °° K 

lim (e Tlt <t>, tffj = lim V I Vj)(ifr m \ vj) + lim V — ^'((f> \ Mj )(^ n \ //,-> 

m—yco m— >oo i—l 2.tL ; m— »oo i—i 2.L) ; 

j=0 1 j=0 J 

M K °° K 

= Yj tf^^ I lim I ^> + X ^?r^'<^ I w> lim I 

Ztt, ; m—yoo t—l 2.L) ; m— »oo 

7=0 1 7=0 (5.43) 

Since | i// m ) — » with respect to the weak dual topology, the left hand side above is zero: 

\im(e Tit <p,^ m ) = \im(e Tit </>\~fi m ) = 0. 



Therefore, we obtain 

= V — lim <<//,„ I v y > • | v,-) + V — lim <<//,„ | // •> • | ^>. 
7 =o j=Q LU ) 

Now Thm.5.8 (ii) shows that (see also Eq. (l5.36l) below) 

lim<^ m |^> = (5.44) 

for all j. This proves that 



lim Il)' a \ ifsj = lim — (ifs m \Hj) ■ \fij) = 0, 

m—>oo J m— »oo 



with respect to the weak dual topology. This implies that ITj' a is sequentially continuous 
at e i(V c , a ) with respect to the weak dual topology. Since Exp + is Montel, a weakly 
convergent series in Exp+ also converges with respect to the projective topology. Thus 
Hj' a is continuous at e KV c ,a) for all j with respect to the projective topology. Note that 

if I (/>,„) £ i{Vc, a ), Eq. (l5.44l) may fail in general as is shown in Eq. (l5.38l) . 

(iv) To prove Zy s Vi,o, recall that Zy is defined by Eq. (l3.1l) . We want to estimate the 
analytic continuation of Z;(t, oj) with respect to co. Put X(t) = e^* M, From Eq. (l2.3l) . 
it turns out that X satisfies the equation 

at 2 2 

X(0) = e^ e . 

PutZ(0 = ^(t)e^ p(t \ ?](t) = ((fa* 1 **® with £,£,p,q £ R. Then, the above equation is 
rewritten as 

f + +*irJ = < V=lReM - Im(c))£ + f _ *^V^\ 

at at 2 2 
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which yields 

§ = -Im(w)f + - e) cos(q - p). (5.45) 
at 2 

This equation shows that if Im(a») > and |£| = 1, then d^/dt < 0. Therefore, if the initial 
condition satisfies |X(0)| < 1, then \X(t)\ < 1 for any t > and any Im(a») > 0. Thus the 
analytic continuation of Zj(t, u>) to the upper half plane is estimated as 



\Zj(t, w)\ < f \X(tWh(6)d9 < 1 , j = 1,2, •• • , 
Jr 



which means that Zft, to) e V ifi for every t > 0. ■ 

Eq. (l5.33l) gives the spectral decomposition for an element | iff) e z'(Exp + ). Now the 
spectral decomposition for £ e i(Vc, a ) is obtained as follows. For £ e i(Vc, a ), there exists 
a sequence {| iA/)}^ c i'(Vc,a) such that | ^) — » Then, from the proof of Thm.5.10 (iii), 
we obtain 

(e T,t ) x t = lim(e Tlt ) x \^i) 

j->oo 

= lim V -=e Ant (it/j\fi n ) ■ \/j, n ) 
K 



Y^ e v lim<lA , H x {5A6) 



and 



M 



£ = V — \im(i/fj\ Mn ) ■ \/u n ) + lim^M^]. (5.47) 

t—t O T\ ;— >c» l-)0O 

From Eq. (l5.47l) . it turns out that the projection IT,, : Exp^ —* Exp' + is given by 



K 



Tint = — lim<^ | • I i"n>, <A; 6 V C ,a, I */fj) ~> £ (5 AS) 



for any £ e i(V c , a ). 

Define the generalized center subspace of T[ to be 



E c = span{^„ | A n e V^R} c Exp^. (5.49) 



Proposition 5.11. E f is a finite dimensional vector space. 

Proof. We have to prove that the number of roots of Eq. (l4. Ill) on the imaginary axis is 
finite. However, it is obvious because D(A) is analytic and tends to zero as \A\ — > oo along 
the imaginary axis. ■ 
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Now we suppose that K = K c . Then, the operator T { has resonance poles A , ■ • ■ ,A M 
on the imaginary axis. Let /i , ■ • • be corresponding generalized eigenfunctions; that 
is, E c = span{// , ■ ■ ■ ,Hm)- Let be a complementary subspace of E c in Exp+ includ- 
ing ftM+i, f*M+2, • • • • In this case, Eq. (l5.33l) gives the decomposition with respect to the 
direct sum Exp^ = E c © E^. Let II f : Exp^ — » E c be the projection to E c with respect 
to the direct sum Exp^ = E c © E^. Although E^ may not be unique, ILJ if/) is uniquely 
determined for | if/) e ?'(Exp + ) because of Thm.5.8 (iii). The complementary subspace E^ 
including fi M+ i,fi M+2 , • • • is called the stable subspace. Eq. (l5.30l) shows that H s (e Tlt ) x \ if/) 
decays exponentially as t —* oo, where H s = id - H c is the projection to E^ . 

Theorem 5.12. Let /i , ■ • • ,p n be generalized eigenfunctions associated with the reso- 
nance poles on the imaginary axis. For any | if/) e z'(Exp + ), the projection to the center 
subspace IL satisfies 
(i) 

U c \if/) = U \if/) + --- + U„\if/) 

K .— — . K 



(ii) 



(iii) 



-=-(>// \fi ) ■ \Mo) + ■■■ + -=-{ift \n„) ■ \Hn)- (5-50) 
2D 2D n 



n c T*\if/) = T*n c \if/), (5.5i) 

U c (e Tlt r\^) = (e Tlt m c \^). (5.52) 



IT; = n c . (5.53) 

(iv) n c is continuous on i(V c , a ), and Eq. (15.50l) is continuously extended to any £, e i(Vc, a )- 

Proof, (i) and (iii) follows from the definition, (ii) is verified by using Eq. (l5.7l) and equal- 
ities T*fi n = A n iJL n , (e Tlt ) K fi n = e Ant n„. (iv) follows from Thm.5.10. ■ 

For Theorem 5.12, it is sufficient to assume that resonance poles on the imaginary 
axis are poles of first order instead of the assumption (A4). If resonance poles on the left 
half plane are poles of higher order, the projection to the complement E^ may become 
complex, however, the projection to E c is still given as above. For the bifurcation theory, 
it is sufficient to assume this fact. 



5.3 Spectral theory on F p c L 2 (R, g(a>)da)) c F' 

In the previous section, we have developed the spectral theory on Exp^, and the order 
parameter is completely resolved: it is decomposed into a countable set of exponential 
functions of the form e A "' . However, for the bifurcation theory, it is sufficient to construct 
more coarse decompositions. To develop the center manifold theory, it is sufficient to 
express the semigroup as the sum of a neutrally stable term and an exponentially decaying 
term. To do so, we expect that it is sufficient to investigate resonance poles near the 
imaginary axis. For this purpose, let us introduce a space F p , which consists of certain 
holomorphic functions defined near the real axis, and study the spectral theory of the 
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semigroup e Tlt on F p . On this space, we can remove the technical assumption K $ 7<"i 
because we do not use resonance poles near infinity. 

For a fixed positive number p > 0, let F P (J3, n) be the set of holomorphic functions on 
the region {z e C | - 1 /n < lm(z) < p + 1 In] such that the norm 

||0IU„:= sup \(f>{z)\e-^ (5.54) 

-l/n<Jm(z)<p+i/n 

is finite. With this norm, F p (J3, n) is a Banach space. Let be their inductive limit 

with respect to n = 1, 2, • • • 

F p (J3) := Urn F p (0, n) (= [j F p (J3, n)). (5.55) 

n>\ n>l 

Next, define F ; , to be their inductive limit with respect to/3 = 0, 1, • • • 

F p := Urn F p (fi) (= [j F p {fi)). (5.56) 

/3>0 /3>0 

Obviously Exp + (/?) c F p (/3) and Exp + c F p c L 2 (R,g(ct>)da>). The analysis of the space 
F p and the operator T\ on F p is done in a similar way to that of Exp + . By the same way 
as the proof of Thm.5.3, we can show that F p (J3) and F p are complete Montel spaces. 
The topology on F p is stronger than that of L 2 (R, g(oj)dco), and F p is a dense subspace of 
L 2 (R, g(aj)du)). Thus the rigged Hilbert space F p c L 2 (R, g(co)du)) c i 7 ^ is well defined. 
On the dual space F' , three topologies are defined; the weak dual topology, the projective 
topology and the strong dual topology. Since F p is Montel, on any bounded set of F' p , 
these topologies coincide with one another. In particular, a weakly convergent series on 
F' p also converges with respect to the projective topology and the strong dual topology. 
The projective topology is metrizable with the distance d defined by Eq. (l5.42l) (we use 
the same notation as before). 

As Prop. 5. 5, we can verify that T\ is continuous as an operator on F p , and it gener- 
ates a semigroup e T[t : F p — » F p . Unlike the semigroup on Exp + , e Tit on F p is not a 
holomorphic semigroup. However, we do not need this property because the spectral de- 
composition of e Tit on F p shown below holds for every t > 0, although the decomposition 
(14.181) on Exp + holds for t > t Q . 

The resonance poles A„ and the generalized eigenfunctions p n are defined in the same 
way as in Sec. 5. 2, although only the resonance poles which lie on the strip {z e C | — p < 
Re(z) < 0} will be considered. It is easy to verify that p n is continuous on F p , and thus 
fi n € F' p c Exp^. Prop. 5. 6 and 5.7 also hold on F p . 

To develop the spectral theory on F p , we start with the Laplace inversion formula like 
as Thm.4.4. 

Proposition 5.13. Suppose that < K < K c and g(uS) satisfies the assumptions (Al) to 
(A3). For a positive number < a < p, suppose that there are no resonance poles on the 
line Re(z) = -a. Then, for any (p,if/ e F p , the equality 

" e~ at r°° 

(e r 'V,<A) = J^Rj[t,(/),if,]eV + — p.v. e^'F^-a + <^\y)dy (5.57) 



49 





C3 1 


Ci 




> 

c 4 

-a 


* 


r 

> 




y 




r 


3 ► 






Q 





Fig. 10: The contour for the Laplace inversion formula. 



holds for t > 0, where A$, • • • ,A„ are resonance poles of T x in the strip {z \ -a < Re(z) < 0}, 
Rj[t, (f>, if/]e Ajt are corresponding residues of the function F l (A)e /lt , and where F^(A) is de- 
fined by Eq.(l431). 

The contour used in the proof of this theorem does not pass close to infinity in the 
domain n/2 + 6 < arg(z) < 3n/2 - 5. Thus we need not introduce the set "K of asymptotic 
values, and the proof is simpler than that of Thm.4.4. 

Proof. The proof is similar to that of Thm.4.4. Let d > be a sufficiently small number 
and r > a sufficiently large number. Take paths C\ to C 6 to be 

C\={d+ ^^Ay \ -r<y<r], 
C 2 = {x+ V^r|0 < x < d}, 
C 3 = {x+ V^Ir | - a < x < 0}, 
C 4 = {-a + V^ly | - r < y < r}, 

and C5 and are defined in a similar way to C3 and C2, respectively, see FigflOl 

Note that if r is sufficiently large, all resonance poles in the strip {z\ - a < Re(z) < 0} 
are included in the closed curve C = J)f=i C, because there are no resonance poles on the 
area (14.281) . In particular, the number of resonance poles in the strip is finite. 
By the residue theorem, we have 

27T^J]Rj[t,(l>,i/f]e^ t = [ e At F (A)dA+ [ e^F^dA. (5.58) 

JCi+d+C! JC3+C4+C5 

By using \(A - V-lAfT 1 ! ~ 0(1 /\A\) and the assumption (A3), it is easy to verify that the 
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integrals along C 2 , C 3 , C 5 and C 6 tend to zero as r — » oo. Thus we obtain 



(e Tl, <p,lp) = Y RAU&Me^ + \im — e M F x (A)dA 

r ^2nypl J-a-yRr 

- e~ at C 

= y RAt, (p, if/]e Ajt + lim — e^'F^-a + V^IyMy. (5.59) 
j-J r^oo 2n J_, 

We have to show that the principle value of the integral in the right hand side exists. It is 
easy to verify that 

((A - V^lMrV,^) = j J \f>(co)ifr(co)g(aj)daj + O(^) 

as \A\ — » oo outside the imaginary axis. Put A := f R cf>(a>)if/(u))g(u))du). Then, the definition 
of Fi(A) provides 

Fl( _ a + V^Tv) ~ -^=r- + 0(l/y 2 ) (5.60) 
as y — » +oo. Therefore, for a sufficiently large M > 0, the integral is estimated as 

Xco 
oo 

J-M Jm V— lv J~oo v-ly 

= f e ' riyt F l (-a+ V r Iy)rfy + 2A f ^^Jy + 0(1/M). 

J-m Jm y 



J -'CO 
sin(y?)/y<iy exists for ? > 0. This proves Prop. 5. 13. ■ 

Since F\(-a+ V-ly) is smooth iny, the integral in Eq. (l5.57l) tends to zero as t —* oo by 
virtue of the Riemann-Lebesgue lemma. When < K < K c , resonance poles Aq, ■ ■ ■ ,A n 
satisfy Re(/l 7 ) < 0. This proves that (e T,t (p, if/) converges to zero as t — » oo. In particular, 
the order parameter r/(t) = (e Tlt <p,Po) for the linearized system (13.41) tends to zero as 
t — » oo. This fact was proved in Thm.4.4 under the assumption that K i ( Ky, and now the 
assumption is removed. 

By using the canonical inclusion i : F p — » F' p , the left hand side of Eq. (l5.57l) is 

rewritten as {<f> \ (e T[t ) x if/) . On the other hand, by using the functional ju(A) defined in 
(15.161) . Fi(A) given by (14.51) is rewritten as 

Fi(A) = + f -^(<P\^X))m^X)), (5-61) 
where the function G(A) is defined to be 

G(A) = 1 - ^D(A) - nKg{- V^U), (5.62) 
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and the functional iff ■ 11(A) is defined to be 



(4> | iff ■ nW) = f l -—tffa)Mo)g(a>)d(o + 2n<K- V^U)^- V^U)g(- V^U). 

JrA- y-lco 



As a result, we obtain the spectral decomposition of (e T[t ) x on F „. 



(5.63) 

p* 

Theorem 5.14. For Eq. (l5.57l) . assume further that resonance poles Aq, ■ ■ ■ , A„ are poles of 
first order. Then, the equality 

(e r ") x l^> = J\—^ t mF J )-\tJij) + ^- l e-^ I > t \j}-iu(-a+ yTAy))dy 

•_Q 2Dj ^ J-OO 

+ e-^ yt =(i/r\ix(-a + V^ly)) • \n{-a + ^ly))dy 

47T J_«, G (_ fl+ yZT y ) 

" i r—a+ V— I"" 

= y— <^'<^h^>+ — = ^ifM^p 

r*—ii+ Y— Too i 

+ = (5.64) 

4^V r l J-«-V^Too G(/l) 

holds for any iff e F p and f > 0. In particular, since (e r ' f ) x |iA) is continuous at t = 0, 
taking the limit yields 

l^> = y — <^>-|/i;> + =lim e^-/^))^ 

2D j In V 3 ! f ^ +0 J -a- V^Too 

+ — lim e At ^=(if/\fi(A))-\n(A))dA, (5.65) 



4tt V 3 ! '^ +0 V^Too G(A) 
which gives the spectral decomposition of | if/) e i(F p ). 

If Aq, • • • ,A n are poles of higher order, the first term in the right hand side becomes 
more complex involving polynomials in t. In this paper, we avoid such a situation for 
simplicity. Since Eq. (l5.64l) comes from Eq. (l5.57l) . the integrals in the right hand side of 
Eq. (l5.64l) converge with respect to the weak dual topology on F' . However, since F p is 
Montel, they also converge with respect to the strong dual and the projective topologies. 

A similar result to Thm.5.8 is stated as follows. 

Theorem 5.15. (i) A system of generalized functions 

T := {jxj, n{A) \j = 0, • • • ,n, A = -a + V^Tv,y e R} (5.66) 

is complete in the sense that if (iff\fij) = (iff\fJ.(A)) = for j = 0, • • • ,n and A = -a + 
y/^y^y e R, then if/ = 0. 

(ii) Suppose that a functional £(A) parameterized by A e C satisfies 

/->-«+ V-Too f*- a+ V— loo 

(e T,t ) x I £(A)dA= I e A! t(A)dA e F' p . (5.67) 

J -a— V-Too J —a— V-T°° 
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If Z"=o a j\ /f/) + Xa'^vSr = f° r a i e then a; = (J = 0, • • • , n) and 

(iii) Let V be a complementary subspace of span{// , ■ • • m ^» which includes any 

functionals of the form f_"^Jzj™ £(A)dA satisfying Eq. (l5.67l) . Then, any | if/) e i(F p ) is 
uniquely decomposed with respect to the direct sum span{// , ■ ■ ■ ,fJ- n } © V as Eq. (15.65l) . 
and this decomposition is independent of the choice of V. 

Proof, (i) If (if/ 1 fxj) = (if/ 1 fj.(A)) = 0, we obtain 

r'—a+ V-T°° 

(PoK^'OV) = = e M {P \^-n(X))dA 

2n\ — l J-a-V^Too 

= — ^<«AIM^))^ = 0. 

2n\ — l J-o-V^oo 

Differentiating in ? and taking the limit t — > yield (T"P = for n = 0, 1, 2, • • • . 
Since r".Po is a polynomial of degree n in <x>, we can prove by induction that {of \ if/) = 
for n = 0, 1, 2, • • • . Since the set of polynomials is dense in F p , ((f> \ if/) = for any <f> 6 F p . 
This proves if/ = 0. Parts (ii) and (iii) are proved in the same way as the proof of Thm.5.8. 
■ 

Because of Thm.5.15, we can define the projection n ; : i(F p ) — » span{/j ; } to the 
generalized eigenspace by Eq. (l5.37l) . Now we suppose that K = K c . Then, the operator 
Ti has resonance poles on the imaginary axis, say Aq, ■ ■ ■ ,A n . Take a positive number a 
so that all other resonance poles satisfy Re(/Jy) < -a, j = n + 1, n + 2, • • • . In this case, 
Eq. (l5.65l) gives the decomposition with respect to the direct sum F' p = E c © E^-. The 
projection n c to the generalized center subspace E c is given by n c = n + • • • + n„ as 
before. The complementary subspace E^ satisfying the condition given in Thm.5.15 (iii) 
is called the stable subspace. Eq. (l5.64l) shows that U s (e T,t ) x \if/) decays faster than e~ at , 
where Tl s = id - II C is the projection to E^-. 

Unfortunately, the projection n c to the generalized center subspace E c . is not contin- 
uous on i(F p ), however, we can show that it is continuous on a subset of i(F p ) like as 
Thm.5.10 (iii). Define a subset V P (C, a) c F p (a) to be 

V P (C, a) = {(f) e F p (a) \ \(p(z)\e~ alzl < C when < Im(z) < p\. (5.68) 

Note that V c , a defined by (15.391 ) is included in V p (C,a). Then, we can prove that IL : 
i(V p (C,a)) —* E c is continuous with respect to the projective topology on F p . Indeed, 
take a sequence {if/ m }™ =l c V P {C, a) such that i{if/ m ) = | if/ m ) —* as m —* oo with respect 
to the weak dual topology. Consider the formula (15.581) for each if/ m . We can prove that 
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the integrals along C 2 to converge as r — » oo uniformly in m. This fact yields 
= ]im(4>\(e T "r^J 

m— >oo 

- K e~ at C 

= lim V —-e Ajt ((f>\iUj)(if/ m \fij) + lim lim — e^'F^-a + ^\y)dy 

7=0 3 U ~' 

^— ' ZL) i m— >oo r— »oo J m— >oo 

j'=0 

Taking the limit ? — > oo proves the equality (15.441) for j = 0, • • • , n, which means that 
ITo, • • • ,n„ are continuous on i(V p (C, a)). The detailed proof is the same as that of 
Thm.5.10 (iii) and left to the reader. This also shows that the spectral decomposition 
holds for £ 6 i(V p (C, a)) like as Eq. (l5.47l) . These facts are summarized as follows (see 
also Thm.5.12). 

Theorem 5.16. Let po, ■ • ■ ,Pn be generalized eigenfunctions associated with the reso- 
nance poles on the imaginary axis. For any | <//) e i(F p ), the projection n c to the center 
subspace satisfies 
(i) 

n c |<//> = n |iA> + --- + iW> 

g 

= ^=-{^\Hq) ■ \Mo) + • •• + ^^(t/f\Mn) ■ \Mn>- (5.69) 
2D 2D n 



(ii) 



n c r 1 x |<A) = r 1 x n c |^>, (5.70) 

n c (e Tit r\^) = (e Tlt fn c \^). (5.71) 



II] = n c . (5.72) 

(iv) n c is continuous on i(V p (C, a)) with respect to the projective topology. Thus Eq. (|5.69b 
is continuously extended to the closure i(V p (C, a)). 

We have shown that II C is continuous on i(V p (C, a)). In general, there are no constants 
M such that the inequality ||iL£||* n < M\\^ n holds for £ 6 i(V p (C, a)) because i(V p (C, a)) 
is not closed under the multiplication by a scalar. To obtain an estimate of ||Il c £||^ n , we 
need an additional assumption for g(a»): 

(A5) For fixed p > 0, g(co + V^Ip)/g(<x>) e F p . 

For example, the Gaussian distribution satisfies this assumption for any p > 0. 
Next, for C > 1 and a > 0, we define a subset W P (C, a) of F p to be 

W P (C, a) := {(A 6 F p (a) | \xj/{z + V^/iAfe)^ 1 < C, when < lm(z) < p}. (5.73) 

In other words, if/ e W P (C, a) if and only if <//(z+ ^Ip) I '<A(z) 6 V P (C, a). Let i(W p (C, a)) c 
F' p be a closure of i(W p (C, a)) = \ <// e W P (C, a)} with respect to the weak dual topol- 
ogy. If £ e i(W p (C, a)), then kg e i(W p (C, a)) for any k e C. 
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For elements in i(W p (C, a)), let us estimate the norm of the projection Ily. 

Lemma 5.17. Suppose that g{oS) satisfies (Al) to (A3) and (A5). 

(i) For each £ 6 i(F p ), ||£||^„ is bounded as n — » oo. 

(ii) For every f3 = 0, 1, • • • and n = 1,2, • • • , there exists a positive number Q(J3) such that 
the inequalities 

\\m,i < \\m,n> mi* z Qm\m,i (5.74) 

hold for £ e i(F p ) (this means that norms || • \\^ n are comparable IfTTlO . 

(iii) For /u(A) defined by Eq. (l5.16l) . the linear mapping \iff) i-> (ifs\/j.(A)) from i(W p (C,a)) 
into C is continuous with respect to the projective topology when —p < Re(A) < 0. In 
particular, if a resonance pole Aj satisfies —p < Re(/l ; ) < 0, the corresponding projection 
n ; - is continuous on i(W p (C, a)). 

(iv) For every {3 = 0, 1, • • • and n = 1, 2, • • • , there exists a positive number D Ca pj such 
that the inequality 

Iin,^„ < D c ^M\\h (5-75) 



holds for £ei(W p (C, or)). 

In Thm.5.10, we have proved that ITy is continuous on i(V p (C, a)). This lemma means 
that under the additional assumption (A5), Hj proves to be continuous on i(W p (C, a)). 
Since the set i(W p (C, a)) is closed under the multiplication by a scalar, we can obtain the 
inequality above. 

Proof, (i) \\£\\l n has an upper bound 



M\\*/3,n SUp 

Mh.n = l 



f (f>(a>)t(co)g(co)daj < f ^(oj)\g(co)dco, 
Jr Jr 



which is independent of n = 1, 2, • • • . 

(ii) The inequality \\g\\* { < ||£||^„ follows from the definition. It is easy to verify that 
the inclusion F p (0, 1) — > F p (J3,n) is continuous. Thus its dual operator from F p (fi,n)' 
into F p (0, 1)' is continuous. This shows that there exists a positive number Q(fi) such that 
\\£\\p n ^ Q(fi)\\t;\\l i- Since the norm || • ||* ; I is bounded as n -* oo, we can take Q{fi) not to 
depend onn= 1, 2, • • • . 

(iii) Let {|<A m >}^ =1 C i(W p (C, a)) be a sequence converging to zero as m — > oo with 
respect to the projective topology. By the definition of the projective topology, we have 
|| \if/ m )\\p n — » for every /3 and n. This means that (/ 1 ifs m ) — * uniformly in / e F p (J3,n) 

satisfying \\f\\^ n < C for each C > 0. Due to the part (i) of the lemma, || \if/ m )\\*p n — > 

uniformly in n = 1, 2, • • • , which shows that (/ 1 t^ m ) — > uniformly in / 6 V p (C,/3). For 
a positive number p > satisfying -p < Re(/l), (i/r m | yu(/l)) is given by 

(i/r m \fi(x» = f + 2^ m (- V-U)g(- V-U) 

JR/I - V-lw 

= f -=^ty m (u + V-fp)g(w + yf-ip)d(o 

JrA - y-l(a) + V-ip) 

l if/ m (o) + V r T/?)g(w + V 3 !/?) 



L 



rA- V r T(o;+ V 11 !/?) ^ m (w)g(w) 



if/ m (io)g(co)dio, 
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where we used the residue theorem. Putting 



j m {o>) :- — — 

A - V 3 !^ + ^h^p) iff m ((o)g((o) 

provides (\f/ m \p{A)) = (f m \ \j/ m ). By the assumptions (A5) and if/ m e W P (C, a), there exist 
C > 1 and a' > such that /„, e V p {C',a') for every m. Therefore, (tfr m \p(A)) — > 
as m — > 0. Since the projective topology is metrizable, this implies that the mapping 
h* (if/\p(A)) is continuous. 

(iv) Since II ; is continuous on i(W p (C, a)) with respect to the metric d, for any s > 0, 
there exists a number £ CjQ ,j > such that if d(^, 0) < 6 c , a ,j, then d(JI£, 0) < s for 
£ e ?'(W p (C, or)). For e > 0, take £ 6 /(Wp(C, a)) and numbers /3, n and 8 c , a ,pj such that 
ll£ll^„ ^ dc,afij- We can take /? and n sufficiently large and S c , a ,pj sufficiently small so that 
d(£, 0) < 8c,aj holds, which implies d(H£, 0) < s. By the definition of d, this yields 

2 n K 2 8 e 
If £ e /(VK /? (C, a)), then £/||£|£„ e i(W p (C, a)). Thus we obtain 



n 8c,afi,& 



3 W%,n 



< t^-, (5-76) 

flu 1 - 2"* 



for £ 6 i(W p (C,a)), which yields Eq.d577j) by putting Dc, a ^,j = 2"k/(1 - 2 n K)/5 c , a fij. 
Because of the part (ii) of the lemma, the same inequality holds for any J3 and n with 
a suitable constant D c ,a,/3j- Since the norm || • \\*^ n is bounded as n — » oo, we can take 
Dc,aj3j n °t to depend on n = 1, 2, • • • . The result is continuously extended to the closure 
i(W p (C,a)). m 

Since the projection Tly and the mapping h-» (if/ \ p(A)) is continuous on i(W p (C, a)), 
the spectral decomposition (15.651) also holds for£ e i(W p (C, a)). In particular, let {iff m }Z=\ c 
W P (C, a) be a sequence such that | if/ m ) —* g e i(W p (C, a)). Then, the projection of £ is 
given by 



n £ = lim (ifr m \ pj) ■ \ pj). (5.77) 

2Dj m - >0 ° 

Next purpose is to estimate norms of semigroups. At first, we suppose that K < K c . 
In this case, there are no resonance poles on the imaginary axis and thus n,~ = 0. 

Proposition 5.18. Suppose that < K < K c and g(oj) satisfies (Al) to (A3) and (A5). 
Take the constant p, which is used to define F p and W P (C, a), sufficiently small. For 
every /3 = 0, 1, • • • and n = 1, 2, • • • , there exist positive numbers M c , a ^ and a such that 
the inequality 

\\(e Ti ra;,n < M C , aJ3 e- jat m\;, n , j = 1,2, • • • (5.78) 



holds for £e i(WJC,a)). 
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Proof. At first, we show the proposition for j = 1. Because of the assumption (A3), there 
are no resonance poles of T\ on the area 



{z 6 C | |arg(z) - tt/21 < 6, |arg(z) - 3n/2\ < 6, \z\ > r), 

if r > is sufficiently large. Further, the left hand side of Eq. (l4.1 II) is analytic. This shows 
that resonance poles can not accumulate on the imaginary axis. Therefore, there exist 
positive constants a and a' such that all resonance poles satisfy Re(/L.) < -a' < -a. For 
the definition of F p and W P (C, a), we choose the constant p so that Re(A n ) < -p < -a'. 
Then, 

\\(e Tlt rt\\;, n _j(e {Ti+a)t ra;, n 



and (e (ri+fl) x £/||£||;„ is given by 



ie^T-^rr = — e^'\-§-- l x{A))dA 



jyr /^-fl'+V-Too 
JL_ \ ^a)tJ_ { _L 

An^TA J-a'-yRoo G{A) Wp,n 



+ f= e (A+a), =(-^—\ f i(A))-\iu(A))dA, 



for £ e i(F p ). Since -a' < -a, this quantity tends to zero as t — > oo for each £ 6 /(Fp). 
Because of Lemma 5.17, the right hand side is bounded with respect to the norm || • ||* n 
uniformly in £ e i(W p (C, a)) and t > 0. This proves that there exists a positive constant 
L c ^ such that ||(e r ' f ) x ^,„ < e~ at L c ^\W pn for f G i(W,(C, a)). Since the norm || • 
is bounded as n — » oo, we can take L Ca p not to depend on n = 1, 2, • • • . The result is 
continuously extended to the closure i(W p (C, a)). 

Next, let us consider Tj for j = 2, 3, • • • . Let z 1; z 2 , ■ • ■ be poles of g{co) on the upper 
half plane. Take p > so small that p < Im(z„) for every n. Then, Cauchy's theorem with 
the aid of the assumption (A3) proves that 

<0|(^W> = (e^^) 



= \ e jyrici, <p(oj)il/(to)g{o))doj 
Jr 



J-> y—ip+oo 
e j ^ wt cp(co)ifr(co)g(co)dco 

e j ^ I<JJt (f>(co + ^Ap)i//(co+ V^p)g(co + yPlp)dco, 

CO 



= e~ Jpt 



for any <p, i// e F p . Since g{co + V-Tp) decays faster than any exponential functions 
as to — » +oo, the integral in the right hand side above exists. Now take a sequence 
{|</r m )}^ =1 c i(W p (C, a)) converging to zero as m — > oo with respect to the projective 
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topology. Then, 

oo 



e n <p(oj+y-\p) ifr m (aj)g(co)daj. 

f m ((o)g(a)) 



Because of the assumptions, there exist C > 1 and a' > such that the functions 

lff m (OJ)g{Oj) 

are included in V p (C',a') for every m, 7 and ? > 0. This proves that (0| (e {Tj+jp)t ) x i/f m ) 
tends to zero as m — > 00 uniformly in j = 2, 3, • • • and t > 0. Therefore, the operator 
( e ( r ;+i» f )x j s continuous on z'(V7 p (C, or)) uniformly in j = 2, 3, • • • and t > 0. By the same 
way as the proof of Lemma 5.17 (iii), we can show that there exists a positive number 
N c , afi such that ||(e< W) x £|£,„ < N c , afi \W^ for f 6 /(W^C, or)). Choose the constant 
p > so that it satisfies both conditions Re(/l„) < —p < -a' and p < Im(z„) above. Then, 
putting M c , a fi = niax{L CQ , iyS , N c , a fi} yields the desired result. ■ 

Next, we suppose that K = K c . In this case, there exist resonance poles on the imagi- 
nary axis and n c ± 0. Then, we can prove the following proposition. 
Proposition 5.19. Suppose that K = K c , g(of) satisfies the assumptions (Al) to (A3) and 
(A5), and that resonance poles on the imaginary axis are poles of first order. Take the 
constant p > sufficiently small. Then, for every = 0, 1, • • • and n = 1, 2, • • • , there 
exist positive constants L c , a ^, M c ,a,p, N c , a fi and a such that the inequalities 

ll( e r ' f ) x n^ii;,„ < L c , afi u\\i n , (5.79) 

W(e Tlt TnMU < Mc, a *e- a, M\\fr, (5-80) 



hold for £ e i(W p (C, a)), and the inequality 

s AT _ „-jat\\ t 



^T^IIL < N c ^e- iat U\\* j = 2, 3, • • • , (5.81) 



holds for £ei(Wp(C, or)). 

Proof. Let A ,--- ,A n be resonance poles on the imaginary axis. For £ e i(W p (C, a)), 
(e Ti ') x Tl c t; is calculated as 

I I) - l) I / I) 

bince we obtain 

n v n n 

11 (e r ") x iu n;,„ < X 11 -=-<^> • l ^> = 2 11 \\; M < £ zwjll^- 

7=0 2L) y ;=0 7=0 
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This provides Eq. (|5.79b for £ e i(W p (C, a)) with L c>a ^ = Y/l =0 D Ct(Xt pj. The result is con- 
tinuously extended for £ e i(W p (C, a)). The proofs of Eq. (15.80l) and (15.811) are the same 
as that of Prop 5.18 with the aid of the fact that II S = id-Ii c is continuous on i(W p (C, a)). 
■ 

Note that the constant a > in Prop. 5. 19, which determines the decay rate of (e Tl ') x , 
is chosen so that all resonance poles Aj on the left half plane satisfy Re(/l y ) < -a, and 
all poles zj of g(co) on the upper half plane satisfy lm(zj) > a. If resonance poles on the 
imaginary axis are poles of higher order, ||(e Tl ') x n r £||^ ;J may diverge as t — » oo with a 
polynomial rate. 

The next proposition is used to show that the inequalities (15.78l) , (l5779l ), (15.80l) . (15.81l) 
are applicable to solutions of the Kuramoto model. 

Proposition 5.20. (i) For any C > 1, i(V 2p (C, 0)) c i(W p (3, 0)). 

(ii) Put V = Uc>i ^2 P (C, 0). Then, the generalized center subspace E c is included in i(V); 

E c c Kyj c i(Wp(3,0)) c F' p . (5.82) 



Proof, (i) If a function if/ e V 2p (C, 0) has zeros on the region < Im(z) < p,\J/ W P (C, a) 
for any C and a. To prove that = i(tfr) 6 2X^(3,0)), let us perturb the function 
i//€ V 2p (C,0). Forn = l,2,---,put 

iJ/(oj) = iff((o) + 2Ce <~ h ^ + "P. (5.83) 
Since if/ 6 V 2p (C, 0), if/(co + ^Ip) £ V P (C, 0). For < lm(to) < p, we have 
if/(co + V^p) 

j/((0) 

which implies iff e W p (3,0). It is easy to verify that 2Ce^ na ' + " p —> as n — » oo with 
respect to the weak dual topology. Therefore, i(ifr) —* i(ifi) e i(W p (3, 0)) as n —* oo for any 
C > 1. 

(ii) Put va((d) = 1 1 {A - V-T<x>). Let Aj = yf—Tyj be a resonance pole on the imaginary 
axis. By the definition, the corresponding generalized eigenfunction pj is given by 



pj = lim i(v x+ v=r y /w)) = lim | v x+ ^ (oj)), (5.84) 

JC-»+0 -T— >+0 • J 

where the limit is taken with respect to the weak dual topology. It is easy to verify that 
v A (u)) 6 V for Re(/l) > 0. This implies that pj e i(V) and thus the generalized center 
subspace E c is included in i(V). ■ 

Note that i(V) is a closed subspace of F' p . If if/ e V, i(if/) = \iff) satisfies inequalities 

(15.781) . (15.791) , (15.80l) . (15.81l) . in which the constants depend only on /?. Since E c c i(V), the 
generalized eigenfunctions in E c also satisfy the inequalities with the same constants. The 



if/(co+ V Z Tp) + 2Ce VZT " 



C + 2Ce"" Im(w) 



if/(a>) + 2Ce^~ lnu)+n P 



< < 3, 

_ 2Ce n{ P~ lm(cJ)) - C ~ 
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set i(V) has all properties for developing a bifurcation theory: it is a metric space including 
all solutions of the Kuramoto model and the generalized center subspace. The projection 
Il c is continuous on i(V). The semigroup (e T,t ) x admits the spectral decomposition on it, 
and norms of the semigroups (e Tjt ) x satisfy the appropriate inequalities. By using these 
properties, we will prove the existence of center manifolds in Section 7. 

Remark. Since the generalized eigenfunction /uj is given by Eq. (l5.84l) and since E c c 
lOO, Eq . (157771) provides 



Since UjfUj = <%u,, we obtain 

«a<-» ^(»> i «> = — «» = j * j # * (5.86) 

for i, j = 0, l,-- - ,n. This is formally written as = -p^/y, which means that 

a system of generalized eigenfunctions {yu ,/ii,- ■ ■ ,f* n } forms a complete "orthogonal" 
basis of E c . 



5.4 Spectral theory on (H+, L 2 (R, g(cj)da)), H' + ) 

In this subsection, we suppose that g{co) is a rational function. Let H+ be a Banach 
space of bounded holomorphic functions on the real axis and the upper half plane (see 
Sec.4.3). In this case, H + is not a dense subspace of L 2 (R, g(co)doj), and thus the triplet 
(H + , L 2 (R,g(u))da)), H' + ) is a degenerate rigged Hilbert space. Indeed, recall that a triplet 
is degenerate if and only if the canonical inclusion is not injective. 

Proposition 5.21. The canonical inclusion i : H + — » H' + is a finite dimensional operator 
; that is, i(H + ) c H' + is a finite dimensional vector space. 

Proof. By the definition, 

/(fA)(0) = (0 = (0,^) = f <f>(o>)tKo))g((o)dcj, 

JR 

for 6 // + . Let z 1? • • • ,z„ be poles of on the upper half plane. By the residue 
theorem, we obtain 

f (p((o)iff((o)g(a))d(D + [ (pire^iffire^gire^^yf-ire^de 
Jr Jo 

n 

= 2wV=T j)Res( Zj ), 

7=1 

where r > is a sufficiently large number and Res(zy) denotes the residue of (f>(to)ifr(a))g(co) 
at z,-. Since is a rational function which is integrable on the real axis, the degree of 
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the denominator is at least two greater than the degree of the numerator : g(a>) ~ 0(1 /\co\ 2 ) 
as \o)\ — » oo. Since ft, ft e H+ is bounded on the upper half plane, we obtain 

n 

(ft\ft) = 27rV^T^Res(z 7 ), (5.87) 

as r —* oo. This means that the action of | ft) e i(H + ) on H + is determined by the values 
of ft(cS) and its derivatives at z\, ■ ■ ■ ,z„. In particular, if the denominator of g is of degree 
M, then i(H + ) * C M . ■ 

Since i(H+) is of finite dimensional, the semigroup (e T[t ) x restricted to i(H + ) is a finite 
dimensional operator. This is the reason that Eq. (l4.30l) consists of a finite sum. In what 
follows, we suppose the assumption (A4). Then Eq. (l4.30l) is rewritten as 

M 

(e Tit cf>, if/) = ((f>\ (e T ") x «A) = V TFr e A " , { ( t> I UnW \Hn), (5.88) 

where definitions of D n and | fj. n ) are the same as those in previous sections. Now we have 
obtained the following theorem. 

Theorem 5.22. Suppose that g(a>) is a rational function with the assumption (A4). For 
any ft e H + , the equalities 

m K 

(e^'flft) = Y!— ^l^>-lM»>. (5-89) 

\ft) = y j ^(ft\fln)-\^n), (5.90) 

hold. In particular, a system of generalized eigenfunctions {/i„}^l forms a base of i(H + ). 
The projection n„ : i(H+) — * span{//„} c is defined to be 

U n \ft) = — (ft\JT n )-\ix n ), n = 0,---,M (5.91) 
2D n 

as before. Since i(H + ) is a finite dimensional vector space, n„ is continuous on the whole 
space. Note that solutions Z 1; 7^, ■ ■ ■ of the Kuramoto model are included in H + (we have 
proved that Zy 6 V\fl in Thm.5.10 (iv)). Thus the bifurcation problem of the Kuramoto 
model is reduced to the bifurcation theory on a finite dimensional space, and the usual 
center manifold theory is applicable. 



6 Nonlinear stability 

Before going to the bifurcation theory, let us consider the nonlinear stability of the de- 
synchronous state. In Sec. 4 and Sec. 5. 3, we proved that the order parameter rj(t) = is 
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linearly stable when < K < K c ; that is, the asymptotic stability of rj(t) = is proved for 



the linearized system (13.41) . For a system on an infinite dimensional space, in general, the 
linear stability does not imply the nonlinear stability. Infinitesimaily small nonlinear terms 
may change the stability of fixed points. In this section, we show that the de-synchronous 
state Zj(t) = 0(j = 1,2, • • •) (which corresponds to p t = l/2n) is locally stable with 
respect to a suitable topology when < K < K c . In particular, the order parameter proves 
to decay to zero as t — » oo without neglecting the nonlinear terms. 



Recall that the continuous model (12.11) is rewritten as Eqs. (l3.2l) . (l3.3l) by putting Zj(t, aS) = 
JT e^i e p t (6, u))d6 with the initial condition 

Zj(0,af)= I e" f=iie h(6)de := hj e C. (6.1) 
Jo 

We need not suppose that h{6) is a usual function. It may be a probability measure on S 1 . 



By using the canonical inclusion, Eqs. (l3.2l) . (l3.3l) are rewritten as a system of evolution 
equations on ]T^li F' of the form 



' iz 1 ) = r 1 x |z 1 )-f< J p |z 1 )|z 2 ), 

at 2 

j t \Zj) = T*\Zj) + ^ ( <P |Zi>| Z y _!> - (P | Z t )\ Z J+1 >) , j = 2, 3, • • • , < 6 - 2) 

k \ZjO~)) = hj\P ), 
where linear operators Tj are defined to be 

T 1( f>(co) = ( V^IM + |f )0(w) = ^iaj<p(co) + j{P 1 0>P o M, (6-3) 

and 

7>(w) = V^TjWt^M = yTljco(p(to), (6.4) 

for j = 2,3,-". TJ axe their dual operators. The main theorem in this section is stated as 
follows. 

Theorem 6.1 (local stability of the de-synchronous state). Suppose that < K < K c 
and g(o>) satisfies the assumptions (Al) to (A3) and (A5). Then, there exists a positive 
constant dp such that if the initial condition h(6) of the initial value problem (12.11) satisfies 



\H = 



2k 

" rie h(6)de 



Jo 



<6 P , j=l,2,---, (6.5) 



then the quantities 



(Zj,(p)= [ fe^ j0 <p((o)p t (9,(o)d(odO 
Jo Jr 



tend to zero as t —> oo for every (p £ F P (J3) uniformly in j = 1, 2, • • • . In particular, the 
order parameter 77(f) = (Z\,Pq) tends to zero as t — » 00. 
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This theorem means that the trivial solution Zy = of (16.21) is locally stable with 
respect to the weak dual topology on F' p . In general, dp — » as fi — > oo. One of the 
reasons is that the norm || • \\^ n goes to infinity as ft — » oo. For the case g(co) is a rational 

function, we can show the same statement : (Z ; , <p) tends to zero as t — » oo for every 
<p e H + if the initial condition satisfies (16.51) . in which dp is independent of fi. 

Proof of Thm.6.1. Since we have Prop.5.18, the proof is done in a similar manner to the 
proof of the local stability of fixed points of finite dimensional systems. Eq. (l6.2l) provides 



| Zi(t, •)) = (e T ' (t - t0) r\ Zi(? , ■)> - f f (Po I Z t (s, -))(e Tlit - s) r\ Z 2 (s, -))ds, 
I Z~(7~)> = (e T ^ t0) r\Zjh~)) + j\(Po I Z^))(e T ^- s) ) x \Z hl (s,-)) 



(6.6) 



-<P„ I Ziis, -)K^- S) ) X \ Z j+l (s, -)))ds, 



for < t < t. Since \Zj) e i(Vi, ) c i(V 2p (l,0)) c j(W p (3,0)) for every t > and 
j = 1, 2, • • • , Prop.5.18 is applied to show that there exists M 3 ^p = Mp > such that 

II I ZA^))\\; n < Mpe- ait - t0) \\ | Z^~))\\l n + y f M^-^-'^ll | ZjJ~))\\l n - IN Z^~)>||^s, 



\Zj(t r ))\\; n <Mpe-^ t - to) \\\Z j (to,-))\\ 



jK r' 

z J to 



(lllz^^-^ + lllz^,^,-))!^)^. 



(6.7) 

Take a small constant <5g > such that hj = f n e-i^ e h(9)d6 satisfies Eq. (l6.5l) . Let us 
show that there exist 1 > and Np > such that 



\Zj(t, -)Wp, n <SpNp 



(6.8) 



holds for any t > t and j = 1, 2, • • • . Indeed, || | Zj(t, -))l^„ is calculated as 

up I (p{co)Zj( 

/3,« = 1 Jr 



II I Z/f, •)>!£,„ = sup 



,•(/, aj)g(io)dcL> 



= sup 



r r « 

Jr Jo 



e^ jxme ' (o) h(e)d9g((o)dio 



Let 



n=-oo 

be the Fourier expansion of gV-r/xW);0,w)_ t^j^ we obtain 

II I Z/f, •)>!!«„ = sup V I (f>(o))a jn (t,co)h n g(a))da) 



< 8p ■ sup 



oo p 

™— _™ *^R 



(f>(a))a; n (t, u))g{o))d(x) 
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In general, for any C 1 periodic function f(9) = YiZ=-oo fn.e n6 -> the inequality 

CO 

Yj \fn\ < max |/(0)| + 2 max \f'{9)\ (6.9) 

n=-oo 

holds, where f n are Fourier coefficients of /, see Newman Il34ll for the proof. Applying 
this to the function ^(p((o)e^ jx(t - m ^ g(co)dco, we obtain 



^ I I 4>{u))a jn {t, co)g(co)dco 

< maxl f (f>(co)e^ Mtmaj) g(co)dco +2 max f cp(co)^- 
9 IJ R e J R 39 

Cauchy's theorem with the aid of the assumption (A3) provides 
V I (p(to)a jn (t, co)g(co)dco < I \(p(co)\g(co)dco 

+2 maxl f 0(to+ yfZi a)je ^Mwm^^a)lL e V=T^,o ; ^ + V=T«) (w + V^Ta)dw 
e IJ R 50 

where a is a positive number. Let us show that the second term above is bounded uni- 
formly in j = 1,2 • • • andf. PutX = X(co) = e *-Mt&M m From the proof of Thm.5.10 (iv), 
we can verify that there exist e > and f > such that when t > t , \X(to + V-la) I < 1 —s 
for any t > to- Therefore, the factor jX j ~ l included in the second term tends to zero as 
j — » oo. To estimate the derivative dX(co + V-la)/<90, note that it satisfies the equation 

d IdX , — \ / / — / — \dX , — 

— I —{to + V-la) = [ V-lu) - a - K?](t)X(to + V-la) — (<x» + V-la), 
dt\o9 j v ' 39 



which is solved as 

—(to + V^a) = X\ t=t0 • exp 



f ( V^la; -a- Kt](s)X(co + V^la)) ds 



Put X = te^P, T] = (e^ 9 with £,p,qe R. Then, £ satisfies the equation (15451) . By 
using it, we obtain 



Re 



f ( V^lco -a- Kr](s)X(co + V^la)) ds 

Jtn 



= J (-a - cos(p - qj) ds 

JtQ 



= -a — — -ds + log 

Jtn * S 



W(0 2 
l-^o) 2 
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Since < £(t) < 1 - e for t > to, the second term above is bounded in t and the first term 
tends to -oo as t — > oo. This proves that dX/d8 — » as t — » oo. Therefore, there exists 
such that Eq. (l6.8l) holds for t > t Q . Since the norm || • \ \* Bn i s bounded as n — > oo, we can 
take iVyj not to depend on n. Then, the first equation of (16 71) gives 

II | Z x (t, -))\\; M < SpMpNpe'^ + 13 / p e- a{t ~ s) \\ \ Z x (s, -))\\*pjs, (6.10) 

for t > t . Now the Gronwall inequality proves 

|| \Z~JtP))\\l n < 8 fi MpN fi e^ KM ^ /2 - a)it - to) . (6.11) 

Since Mp and Np are independent of the choice of dp, by taking 5p sufficiently small, this 
quantity proves to tend to zero as t — » oo. Substituting it into the second equation of (16.71) . 
we obtain 

\\\Zji~))\\l n < SpMpNpe~ ja(t - t0) + jK(SpMpNp) 2 f e ~^) e ^KM p N p i2-a)(s-t ) ds (6 U) 

J t 

for 7 = 2, 3, • • • . It is easy to verify that the right hand side tends to zero as t — » oo 
uniformly in j. 

Now we have proved that if the initial condition satisfies (16.51) for each /?, then || | Zj(t, -))||^ 
decays to zero as t — > oo for every j and n. By the definition of the norm || • ||^ n , this means 

that (Zj(t, •), 0) as f oo for every e F P 08). ■ 



7 Bifurcation theory 

Now we are in a position to investigate bifurcation of the Kuramoto model by using the 
center manifold reduction. Our strategy to detect bifurcation is that we use the space of 
functionals F' p instead of the spaces of usual functions F p or L 2 (R, g(co)dco) because the 
linear operator T\ admits the spectral decomposition on F' p consisting of eigenfunctions, 
while the spectral decomposition on L 2 (R, g(co)dco) involves the continuous spectrum on 
the imaginary axis; that is, a center manifold on L 2 (R,g(co)doj) is an infinite dimensional 
manifold. To avoid such a difficulty, we will seek a center manifold on F' . At first, 
we have to prove the existence of center manifolds. Standard results of the existence of 
center manifolds (see ||5l|9l|24l|47l) are not applicable to our system because the space F' p 
is not a Banach space and the projection n c to the center subspace is continuous only on a 
subspace of F' p . Thus in Sec.7.1, the existence theorem of center manifolds for our system 
and a strategy for proving it are given. The proof of the theorem is given in Sec. 7. 2 to 7.4. 
In Sec.7.5, the dynamics on the center manifold is derived and the Kuramoto's conjecture 
is solved. Readers who are interested in a practical method for obtaining a bifurcation 
structure can skip Sec.7.1 to 7.4 and go to Sec.7.5. Throughout this section, we suppose 
the assumptions (Al) to (A3) and (A5) for g(to). Existence of center manifolds for the 
case that g(to) is a rational function is trivial because the phase space i(H+) is a finite 
dimensional vector space. 
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7.1 Center manifold theorem 

Let i'(!F) be a certain metric subspace of the product space YlkLi F'p w ^ n a distance d^, 
and i(T) its closure. These spaces and the metric will be introduced in Sec. 7. 2 and 
Sec. 7. 3. Let O f be the semiflow on i(T) generated by the system (16.21) (for the existence 
of the semiflow, see Sec.7.3). For the generalized center subspace E c of T\ defined by 
(15.49b . put 

CO 

E c = E,x{0}x{0}x---c]~[f;. (7.1) 

k=i 

Let E^ = E^ x F' p x F' p x • • • be the complement of E c . The existence theorem of center 
manifolds is stated as follows. 

Theorem 7.1. Suppose that g(co) satisfies the assumptions (Al) to (A3) and (A5). There 
exist a positive number e and an open set U of the origin in c n^=i F' P sucri that 
when \K - K c \ < s , the following holds: 

1 A A 

(I) (center manifold). There exists a C mapping q : E c — » E c n z'(!F) such that the dim-E f 
dimensional C 1 manifold defined to be 

W c hc = {y + q(y)\yet c }nU (7.2) 

is <D r invariant (that is, $>t(Wf 0C ) C\U c W£ c ). This is called the local center manifold. The 
mapping q is also C 1 with respect to the parameter s .= K - K c , and q(y) ~ 0(y 2 , sy, s 1 ) 
as y, e — > 

(II) (negative semi-orbit). For every £ 6 , there exists a function w : (-oo, 0] —* i(T) 
such that w(0) = £ and O^mCs)) = u(t + s) when t > 0, s < -t. Such a w(?) is called a 
negative semi-orbit of (16.21) . As long as w(£) e £/, w(0 6 Wf . In this case, there exist 
Ci > and a small number b > such that 

rfoo("(0,0)<C 1 e fa . (7.3) 

(III) (invariant foliation). There exists a family of manifolds {M^}^ eW c c £/, parameter- 
ized by £ 6 W^ c , satisfying that 

(i) M f n W£ c = {£}, ^ = U, and M e n M? = if f ± ?. 

(ii) when O t (£) e U, $,(Mj) nt/c M^. 

(iii) 

= {uei(T)nU 

there exist a > b and C 2 > such that ^(O^w), $,(£)) < C 2 e" a '}. 

By Prop. 5.1 1, Wf is a finite dimensional invariant manifold. Part (III) of the theorem 
means that Wf is attracting with the decay rate e~ c ", where the constant a is the same 
as that in Prop. 5. 19. Further, (Ill)-(iii) means that the semiflow near Wf is eventually 
well approximated by the semiflow on Wf oc if t > is large. In particular, if (16.21) has an 
attractor N near the origin, /V is included in Wf oc . Since the topology induced by the metric 
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coincides with the strong and the weak dual topologies on any bounded set, N is attracting 
for both of the strong and the weak dual topologies. Due to the spectral decomposition 
(15.651) . any element |Zi) e i(F p ) is decomposed as 

n 

\Z l ) = Y_ i oryl Hj) + Yt , aj eCJ,e E c \ (7.4) 

i=o 

where |/i ), ■ • • , \f* n ) 6 E c are generalized eigenfunctions associated with the resonance 
poles on the imaginary axis. Then, part (I) of the theorem means that if (| Zi), | Z 2 ), • • • ) e 

Y 1 ~ 0(a 2 ,sa,s 2 ), \ Z k ) ~ 0(a 2 , ea, s 2 ), k = 2, 3, • • • , (7.5) 

as s, a , ■ ■ ■ , a n — » 0, where a = (a ,-- - ,a n ). Substituting Eq. (l7.4l) into the system 
(16.21) with the condition (17.51) . we can obtain the expression of q{y) as a function of 
s, ao, • • • , a n . The dynamics on W c loc is realized by a system of ordinary differential equa- 
tions of ao, ■ ■ ■ ,a n : 

da; 

— = fj(s; a , ■ ■ ■ , a n ), j = 0, • • • , n. (7.6) 

If (16.21) has an attractor N near the origin, N is an attractor of the system (17.61) . In this 
manner, (16.21) is reduced to a finite dimensional problem. Such a technique to investi- 
gate bifurcation is called the center manifold reduction. Part (II) implies that the center 
manifold is characterized by the property that the dynamics on it is sufficiently slow. 

Although we prove the existence of W c loc in i{T), from a physical viewpoint, especially 
we are interested in an initial condition of the form |Z ; -(0)) = hj\P ) (which corresponds to 
the initial condition of the form po(6, oS) = h(6) for the system (12.11) ). Then |Z/(0) e '(Vi.o) 
(Thm.5.10 (iv)). This means that an attractor of (16.21) which is reached from the initial 
condition |Z/0)> = hj\P ) is included in W c loc n UZi Wio)- 

Sec. 7. 2 to 7.4 are devoted to prove Thm.7.1. It is well known that a global center 
manifold uniquely exists only when a Lipschitz constant of nonlinear terms of a system 
is sufficiently small. Thus in Sec. 7. 2, we consider a perturbed system of (16.21) so that its 
Lipschitz constant becomes sufficiently small, while it coincides with the original system 
in the vicinity of the origin. Because of the perturbation, a solution may fall out of z'(Vi,o) 
and go into a larger space. Thus we will introduce the space ( J T , and show that solutions 
(Zi,Z 2 , • • •) of the perturbed system are included in T . We will prove in Sec. 7. 3 and 
Appendix B that the perturbed system generates a smooth flow to prove that the center 
manifold is smooth. Once we obtain the existence of a proper phase space, a spectral 
decomposition of the linear operator, estimates of norms of the semigroups and a smooth 
flow whose Lipschitz constant of nonlinear terms is sufficiently small, then the existence 
of the center manifold is proved in usual way with a slight modification. We demonstrate it 
in Sec. 7.4. In Sec. 7. 5, we perform the center manifold reduction: equations of ao, • • • , &„ 
are obtained and investigated. The order parameter rj(t) is defined as rj(t) = (Z\,Pq) = 
(Po I Zi). On the center manifold, it is written as 

n n 2 

(PoIZi) = ^ a~(fj(P | fij) + (Po \Yi) = J] Y®^> + 0(a2 > £a > £2) ' (7 - 7) 
7=0 7=0 
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where we use (P = 2/K c , which follows from the definition of resonance poles. 
Therefore, if a bifurcation diagram of (17.61) is obtained, a bifurcation diagram of the order 
parameter is also obtained. In this way, the Kuramoto's conjecture will be proved in 
Sec.7.5. 



7.2 Phase space of the perturbed system 

Recall that the trivial solution Zj(t) = (j = 1, 2, • • • ) corresponds to p t = 1 /2n called the 
de-synchronous state (Sec. 3). Since we are interested in bifurcations from p, = 1 /2n at 
K = K c , put p t = l/2n + p t and K = K c + e. Then, Eq. (l2. II) is rewritten as 

dt d6 K An^TA 

+—^=(n(f)e-^ - W)e^ e ) + -4=0KO«" ^ - me^p) = 0, (7.8) 
47rV r I 2^ ' 

where 

77(0= I I e^~ w p t {6,to)d6d(o = I I e^p^^dOdco. 
JrJo JrJo 

An initial condition p = ^(v> is a real- valued measure (signed measure) on S param- 
eterized by co e R satisfying j h(6, to)d6dco = 0. The first step to prove the existence 
of center manifolds is to localize the nonlinear term so that the Lipschitz constant of the 
nonlinear term becomes sufficiently small. For this purpose, let^ : [0, <x>) — » [0, 1] be a 
C°° function, and consider the perturbed continuous model of the form 

§ + Ucopt + -^- m e-^° - W)e^ W ) 
ot oQ y Any^l 

+—^=(jl(t)e-" rXB ~ me+^m + -4=^"^ - n&e^mPt) = 0. (7.9) 
An V^T 2 ^ ' 

If we put 

e^ e p t {Q,to)dQ, (7.10) 

o 



Eq. (l7.9l) yields a system of equations 

JZi / — K c Is K * \ 

-± = V^koZ, + -frj(t) + (-77(0 - -rj{t)Za)x{t), 



= j V^lwZ; + ^ (77(0^-1 - 77(OZ i+1 )i-(0, j = 2, 3, 



(7.11) 



If^ = 1, this coincides with the original system (I3.2I) . (I3.3I) . If Jt"(0 is sufficiently small, this 
perturbation makes the Lipschitz constant of the nonlinear terms of (13.21) . (13.31) sufficiently 
small (Note that when proving the existence of center manifolds, the bifurcation parameter 
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e is regarded as a dependent variable. Thus erj(t) = s{Z\ , P ) is regarded as a nonlinear 
term). A concrete definition of x will be specified in Sec.7.3. Note that Z = because of 
£ n kd, co)d6dto = 0. 

Eq.( |7.9l ) is integrated by using the characteristic curve method. The characteristic 
curve x = x(t, s; 6, to) is defined as a solution of the equation 

^ = co + Ut)e-^ lx - W)e^ x )x(t), (7.12) 

satisfying the initial condition x(s, s; 6, to) = 6 at an initial time s. Along the characteristic 
curve, (17.91 ) is integrated to yield 

p,(6,(o) = h(x(0,t;e,to),to)exp[^ f (77(5)^ ^'^^ + ^e^ x(sM )x(s)ds] 

~ Jo 

+ exp[| J (rj(T)e~ + -^ e <~^^ x (j)d T ] x 

l Kc + ^ {S \ n{s)e-^ sM +me^ sM )ds. (7.13) 

Once x(t, s; 6, to) and 77(f) are determined, this p t gives a weak solution of (17.91) . The 
existence of solutions of (17.1 II) follows from that of the integro-ODE (17.121) and (17.131) . 
which will be proved in Appendix B (see also Prop.7.3). 



In Sec. 5, we have proved that a solution of (I3.2I) . (I3.3I) is included in Vifl. This property 
may break down because of the perturbation^. Thus we define an appropriate phase space 
for (17.1 II) and prove the existence of the flow on it. Define a subspace T of the product 
space Y\J=i F p as follows: (Zi,Z 2 , • • •) e f if and only if there exists a signed measure 
h{6, to) on S parameterized by to e R such that 

(Fl) 

I h{0, to)d6 = 0. 
Jo 

(F2)Forj = 1,2,---, 

X2rr 
e^ je h(6, to)d6 e V (= [j V 2p (C, 0)). (7.14) 
C>1 

(F3) Define Zj for j = -1,-2, • • • by Eq. (l7.14l) . There exist positive constants C and y 
such that 

sup \Zj(to)\ < Ce ljlr (7.15) 

0<Im(6L))<p 

for all j e Z. 

Recall that the set V c F p was defined in Prop. 5. 20. Hence, T is the set of Fourier coeffi- 
cients {Zj}J = , of signed measures included in Y\%\ V whose growth rate in j is not so fast. 
Although Zj(to) is assumed to be holomorphic near the region < lm(to) < p, h(6, to) 
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need not have an analytic continuation to the complex a>plane. 

Example. If h is of the form h(9,to) = hi(9)h 2 (to) with h 2 (to) e V, the assumptions (F2) 
and (F3) are satisfied. In particular, |Zy(o>)| is bounded uniformly in j. Recall that from a 
physical viewpoint, we are interested in an initial condition of the form h(9, to) = h(9) (see 
Eg. (IQ>1. Let 6(9+to) be the delta function at 9 = -to. lfh(6, to) = 5{9+to)-\ /(2n), these 
assumptions are satisfied with C = l,y = p because Zj(to) = e~^i w and \Zj(to)\ = e jIm ^ } . 
If h{9,to) = 9cos(9 2 to), we can verify that \Zj(co)\ ~ 0{e j ~) for to e C as j —* +00, and 
thus (F3) is violated. 

The existence of solutions of Eq. (l7.1 II) will be prove in Sec.7.3 and Appendix B after 
X is specified. In this section, we show that if solutions exist, they are included in T . 

Proposition 7.2. For a given function^ : [0, 00) — > [0, 1] and an initial condition 
in T , suppose that a solution of (17.1 II) exists and a function x(t, s; 9, to) has an analytic 
continuation with respect to 9 and to (these facts will be verified in Sec.7.3 and Appendix 
B). Then, 

(i) the solution is included in T for any t > 0. 

(ii) for each t > 0, /? = 0, 1, • • • and n = 1, 2, • • • , || \Zfit, -))\\*p n is bounded uniformly in 
j e Z, where \Zjt, •)) = i(Zj(t, •))• 

Proof, (i) At first, let us verify (F2). To show that Zj(t, ■ ) 6 V if Z ; (0, • ) 6 V, we use the 
equality 

r 2n r 2n 

I a(9,to)p t (9,to)d9 = I a(x(t,0;9,to),to)h(9,to)d9 
Jo Jo 

a(x(t,0;G,to),to)exp[~^ (^(r)^^™ + W)e^ x(M )x(r)dr] x 

[ KC + 4 f (S) ) (rt s)e ~ % ~' V0: ""' + ~>^ e % ~ U V " : ""') dsd6, (7.16) 

for any measurable function a(9, to), which is proved by substitution of (17.131) . Note that 
ifx(f) = 1, it is reduced to Eq. (l2.5l) . From this, it turns out that Zj is expressed as 

z ^t, to) = J^e^ ljx(t '°- M h(e,to)de 

Jo Jo 2 Jo 

( Kc+ 4 f (S) ) ^ S)e ~ ^ lx(m + W)e ^ x( - m ) dsd9. (7.17) 

At first, let us show that ^V-ufoO;^) 6 y j s p rove( j m fa. same wa y as Thm.5.10 
(iv). Put X(t) = e^ 1 *^^. Then X satisfies the equation 

— = V^ltoX + - (77(0 - r](t)X 2 ) X (t), (7 lg) 

X(0) = e^~ w . 
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Putting X = ^e^ p , rj = ^e^ q with p,q eR yields 

^ = -Im(^ + f^(l - f) cos(p - q)x(t). (7.19) 
at 2 

This equation implies that if Im(a») > and £ = 1, then J^/J? < 0. Since \X(0)\ = 1, we 
obtain \X(t)\ < 1 for any t > and any Im(u>) > 0. This proves that X(t) is bounded on the 
real axis and the upper half plane: X(t) = e ^ x ^-M £ v 2p (l,0) c V. 
Next thing to do is to investigate Y(t) = g-V-LjfeOftw)^ w j 1 j c ] 1 satisfies 

_ = _VZT^--(, (0 y 2 -,(o)m (720) 

y(0) = e-^i . 
Putting y = %e yr ~ lp , T] = Ce^i with p,q eR yields 



^ = Im(w)£ + f «1 - £ 2 ) cos(p + ^f), 

^ = -Re(w) - f Z(g + -) sin(p + q) x (f). 
at 2 t; 



(7.21) 



When |Re(o»)| is sufficiently large, the averaging method is applicable to construct an 
approximate solution. Eq. fl7.21l) is averaged with respect to p to provide the averaging 
equation d^/dt = Im(a>)£, which is solved as %{f) = e lm{oj)t i;{0). Therefore, a solution of 
Eq. d7.21l) is given as 

m = gtafr* + 0(-±—\ (7.22) 
v |Re(w)K 

as |Re(o>)| — » oo. See Sanders and Verhulst 091 for the averaging method. This im- 
plies that Y is bounded as Re(a») — » +oo for each Im(a») and t. Thus g-V-f-tftO;^) g y 
Therefore, the second term in the right hand side of Eq. d7.17l) is in V; the second term is 
bounded as Re(<x>) — > +oo for each j and t. 

Next, we show that the first term in the right hand side of Eq. (l7. 171) is in V. Let 

OO 

e ^mo M = £ ajn (t t0))e +** (7 .23) 

n=-oo 

be a Fourier expansion of e V-i/*W>;0,w) 

oo 

= ^ w)Z„(0,w). (7.24) 

n=-oo 

Since the series d7.23l) converges uniformly in 6, the right hand side of (17.241) exists for 
each oj. Since e ^-^M t ^ e < (0 ) [ s holomorphic, so is a/ n (f, co). By the assumption (F2), Z„(0, co) 
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is also holomorphic. Therefore, the right hand side of (17.241) converges to a holomorphic 
function near the region < lm(co) < p. By the assumption (F3), there are positive 
constants C and y such that 



sup \Z n {Q,co)\ < Ce Wr . (7.25) 

0<lm(w)<p 



This provides the inequality 



sup \y aj n (t,aS)Z n (0,u)) < sup C \ e^ 7 \a jn (t, (o)\. (7.26) 

OsamMs^^t-i 0<ImM<p n tri 

Let us prove that the right hand side exists. Eq. (l7. 121) shows that x(t, 0;6 + 2n,co) = 
x(t, 0; 9, to) + 2n. With this property, we use Cauchy 's theorem to the function e ^M f <°'M 
along the path represented in Fig fTTT a) to yield 



a jn (t,u>) = i- f e -^ne e ^lMtm<o) dg 



c. 



■- r 



V-ln0 y/-ijx(t,0;6,aj) 



In 



c 2 



for n = 0, 1, 2, • • • and j = +1, +2, • • • , where r > can be taken arbitrarily large be- 
cause e V-T/*W>;0,to) j s analytic in # e C. By the same way as above, we can show that 
e ^=iMtfifi- V=Tr^.) is es ti m ated as 

| g V^/tfoO;0- V^Tr,^ _ g (-Im(w)f+r)j + fl/ 1 \ 

MRe(w)K' 

as |Re(w)| — > oo for each t and Im(a»). This provides 

\a jn (t, co)\ < e~ nr U-^ t+ -y + 0(—^—)) , 
\ v |Re(w)l7 

for n = 0, 1, 2, • • • . When « < 0, we take a path represented as FigflTTb). which yields 

M^)l^™( e '™ + o(^)) 

as |Re(o»)| — » oo in the same way. Therefore, we obtain 

CO 

sup V a ;n (£, <x»)Z„(0, co) 

Im(w)<p l ;J= _ oo 

sup C y e M(y ' r) [ e^^+'W + ^-Wt-r)j + Q t — * — \\ 
;im£< P £^ \ l |Re(a>)K/ 
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0<lm(&>)<p' 

< 

0<lm(w)<p n=l 



By taking r > y, it turns out that the right hand side of Eq. fl7.26l) exists and bounded as 
Re(<x>) — » +00 for each j and t. This proves that Zj(t, co) is bounded as Re(a») — » +00 and 
Z/(f , co) e V for each j and 
To verify (F3), put 

A(s,0) = exp[-y f {rj(j) e -^(jm<o) + ^7) e ^(r,o;^)^ (r) j T ] x 



Then, Zy is rewritten as 



00 w pin 

Zj(t,a>)= V a jn (t,to)Z n (0,co)+ I I e^^^A^s, 6)d6ds. (7.28) 

nfdo -'O Jo 

From the above calculation, we obtain 

00 

sup \Zj(t,co)\ < sup sup C y^y-r)(\ e ^^^^)\l + \ e ^t^^m)y\ 

0<lm(w)<p 0<lm(a))<p 0<8<2n 



+ sup sup In f lACj.^j-le^'-^K 

0<lm(w)<p 0<e<2w Jo 



which proves that Zj(t, co) satisfies (17. 151) for some C and y. 
Finally, let us verify (Fl). Note that Eq. fl7.12l) provides 

%(s,O;0M = expf-l^^r),-^' ^ + Jjfa^M)%T)dT]. 
This shows that the second term in the right hand side of Eq. fl7.16l) is rewritten as 

r' K c + Sx(s) p . ^ ^3 , _ ^ V^v(.v,0;^)\ 

Jo 47r Jo <96> v ' 

In particular, when a(6, co) = 1, this value vanishes because x(s, 0; 6, co) is periodic in 9. 
This fact and Eq. fl7T6l) yield 

/-»2;r /-»2;r 

I p t (6,co)d6= I h(6,co)d6. (7.29) 
Jo Jo 

Therefore, if an initial condition satisfies (Fl), so is p t (6, co) for any ? > 0. Now the proof 
of Prop.7.2 (i) is completed. 

To prove (ii), note that when co e R, Zj(t,co) is bounded uniformly in j e Z be- 
cause they are Fourier coefficients of a measure p, on S 1 . This shows that || \Zj(t, -))lle„ is 
bounded uniformly in jeZ. ■ 
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Fig. 1 1 : The contour for obtaining Eq. fl7.27l) . 



7.3 Localization of the semiflow 

By using the canonical inclusion, we rewrite Eq. (l7.1 II) as an evolution equation on R x 

7' 
P 



WZ\F' of the form 



dt 

d- 



-| Zj> = 7%| Z 1 ) + - (s(P I Z X )\P Q ) - K(P Q | Zi>| Z 2 >)^(0, (7.30) 
^| Z ; > = r*| Z,-) + ^ ( (K\Z\)\ Z hX ) - (P Q I Zi)\ Z j+ i))x(t), J = 2, 3, • • • , 



where | Zy) = z'(Z ; ). The trivial equation ds/dt = is added to regard s = K - K c as a 
dependent variable. The operator T 10 is defined by (16.31) . in which K is replaced by K c . 
Note that r 10 has resonance poles on the imaginary axis. 

In what follows, we denote an element of the space Rx \\™ =l F' p by z = (zo, Z\,Z2, • * ■ ), 
where zo e R and {z\,Zz, - • •) 6 n*li ^p- We also denote it as z = (Zfe)^i- A metric on 
this space is defined as follows: The dual space F p {fi,n)' of F p (J3,n) is a Banach space 
with the norm H^,, = sup| M/jn=1 The projective limit F P (J3)' = \im F P (J3, n)' is a 
complete metric space with the metric dp defined by Eq. fl5.41l) . Next, for the projective 
limit F' p = UmFpifi)' we introduce the metric by 



1 



d{f o = > - dp{u) 



(7.31) 



where {P(/3)}J =0 is a certain sequence of positive numbers such that Y!p=o 1/P(fi) con- 
verges. This defines the same projective topology as that induced by the metric (15.421) . 
The constants P(J3) will be determined in Sec. 7.4 so that d(£, 0) plays a similar role to a 
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norm. Define to be 

CO 

doo(z, z) = sup d(z k , z' k ), z = (zo, zu ■ ■ ■ ) e R X FT F' (7.32) 

if 

where d(z*;, z^) for > 1 is the distance on F' defined as above, and d(zo, z' ) = \zo - z' Q \ 
for zo, z' Q e R. Let Q be a subspace of R x n^=i F£ consisting of elements z = (zo, Zi, • • * ) 
such that sup fe0 d(zk, 0) is finite. With this metric ^ is a complete metric vector space. 
Metric vector spaces and definitions of the metrics used in this section are listed in Table 
2. 



F p (J3,n) 


Banach space: \\(p\\^ n = sup \<p(z)\e /3|z| 

-l/n<Im(z)<p+l/n 


F p (fi, n)' 


J — I 1 II L'l 1 ;|: | / I | U \ | 

Banach space: 1 1^1^ = sup \{(p\%)\ 

II<%„ = 1 


F p ffl = limF p (J3,n)' 


W'«-2j 2 -l + ||f-^ 


F' p = \imF p Q3y 




CO 

£cRx[]f; 

k=\ 


doo(z,z) = supd(zk,z k ) 

k>0 


r, kt) 


T is a subspace of nS=i F p satisfying (Fl) to (F3); 
i(!F) is its inclusion into n~ =1 F' p . R x i(T) c 


X ( ~ m) (= R x i{T)) 


D m (z,z') = *"Uo(z,z') 


CO 

Icf]l w 

m=0 


D(w,v) = supD m (w ( ~ m) ,v ( * m) ), 

m>0 

« = uf-D, ■ ■ ■ ) e n:= * ( - m) , = (4 _w) )r=o 6 xi ~ m) 



Table 2: Metric vector spaces used in Section 7. Definitions of the spaces X ( m \ X and the 
constant k will be given in Section 7.4. 



By the definition, it is easy to verify that 

dUz, +z) = d„(z t z', 0) < ^(z, 0) + dM, 0). (7.33) 

A sequence {z (m) = (zi m) )^ }^ =1 in Q converges to x = (x k ) k=Q if and only if z ( Q m) -» x on 
R and \\z k m) - Xk\\*p n —> uniformly in k > 1 for every f3 > and n > 1. On the other 
hand, since a weakly convergent series in F' p also converges with respect to the metric 

d, a sequence {z (m) = (4 )fc=oK£=i in ^ converges to x = (**)£lo if and only if z ( m) -* x 
on R and (cf> \ z ( k m) ) — > (<f> I x k ) uniformly in k > 1 for every <p e F p . It is notable that a 
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mapping / = (/o,/i, • • • ) : Q —* Q is continuous with respect to the metric doo if fa is 
sequentially continuous with respect to the weak dual topology uniformly in k. Indeed, 
let {z (m) }™ =1 be a convergent sequence with respect to doo, which satisfies doo{z {m) , x) — » 

1 oo 

'm=l 



as m —* oo. Since the weak dual topology is weaker than that induced by doo, {z (m ^°° 



weakly converges to x; that is, zi -4 x onR and (0 | z^ m) ) — » (0 1 x^) for every 6 F p . By 
the assumption, fk(z im) ) weakly converges to fa(x). Since F p is Montel, it also converges 
with respect to doo'- doo{fa{z {m) ), fa{x)) — » 0. Since a metric space satisfies the first axiom 
of countability, this means that / is continuous with respect to the metric doo- This fact 
will be used when proving that the semiflow generated by (17.301) is continuous in initial 
conditions. 

For locally convex topological vector spaces X and Y, a mapping / : X — > Y is called 
Gateaux differentiable at x G X if there exists a map df(x) : X — » F such that the map G 
defined by 

G(f , x ) = df(x)x , t e R, x G X, ^ 

G(0, jc') = 0, 

is continuous at (0, 0) G R x X. Further, if df(x) is linear and if / : X — » 7 and rf/ : 
X x X — » y are continuous, / is called C 1 or Frechet differentiable. Note that / can 
be Gateaux differentiable even if it is not continuous (for example, a linear operator is 
always Gateaux differentiable). For our situation, it is convenient to use the weak dual 
topology. As was mentioned in the previous paragraph, we can show that a mapping 
/ = (/o> fu " ' " ) : Q -* Q is Gateaux differentiable with respect to the metric doo if there 
exist maps dfa(z) such that 



lim 

f->0 



{(p\h(z + tz')- fail)) 

{(j) I df k (z)z ) 



0, (7.35) 



t 

for any <p G F p uniformly in k. Further, if dfa(z) are linear and if fa and dfa are sequentially 
continuous with respect to the weak dual topology uniformly in k, then / is C 1 with respect 
to the metric doo. 

Let i{T) be the subspace of n^=i F' p consisting of elements of the form 

(i(Zi), i(Z 2 ), ■ • • ), where (Z x , Z 2 , • • • ) e T. 

Due to Prop. 7. 2 (ii), R x i{T) c ^. Thus with the distance doo, R X is a metric vector 
subspace of Q, and the closure R x i(T) is a complete metric vector space. 

Now we specify the function x(t). Let E c be the generalized center subspace of T w . 

Put 

00 

E c = R x E c x {0} x {0} x • • • c R x Y\ F' p . (7.36) 

k=i 

Let P c : R x fltli F'p ~* E c be the projection to E c defined by 

P c = (id R , n c , 0,0,---), (7.37) 
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(z'J R is the identity on R) and P s = id-P c the projection to the complement of E c . Because 
of Lemma 5.17 and Prop. 5. 20, P c and P s are continuous on R x i(^F) and E c is included 
in the closure R x i{T). Let^r) be a C°° function such that^f) = 1 when < t < 1, 
< x(t) < 1 when 1 < t < 2, and^f) = when t > 2. Taking a small positive constant 
6 1 , we replace %(t) in (17.301) by 

m--=xlp^Yx(W)\l (V.38) 

where 

oo 

z(t) = (s, I Z^~)>, | Z^t~)), ■ ■ • ) e R x Y] F' p , (7.39) 

k=\ 

and || • || £c is a norm on E c defined as follows: An element y e E c is denoted by y = 
(y ,vi,0, •••) withv! = £"=o a,-| //,•>. Then, ||y|| £c is defined to be 

IMk = (|v | 2 + |aol 2 + --- + |a„| 2 ) 1/2 . (7.40) 
Eq. (l5.69l) shows that n c | Zi) is given as 

U c \Z l ) = Y j a j \iu j ), a j = —{Z l \n j ). (7.41) 
7=0 2£) y 

With these aj, \\P c z\\e c is given by 

\\PcZ\\ Ec =(£ 2 + \a \ 2 + --- + \a n \ 2 ) 1/2 , (7.42) 

for z = (e, |Zi>, • • • ) e R X Since E c is a finite dimensional vector space, the topol- 
ogy on E c induced by || • \\e c is equivalent to that induced by the metric d^. With this^(0, 
we can prove the existence of solutions of Eq. d7.30l) . 

Proposition 7.3. Eq. d7.30l) with %(t) given by (17.381) generates a C 1 semiflow (p t on 
R x i{T). That is, for a given initial condition z e R x i(!F)> Eq. d7.30l) has a unique 
solution denoted by (p t (z), which is C 1 in z, on R x i(*F) for any ? > 0. 

Proof. At first, we show the existence of solutions of Eq. (l7. 1 II) on T . For a given initial 
condition (Zi(0, ■),■■■) eT , there exists a signed measure h(8, oS) satisfying (17. 141) . Such 
a h is uniquely determined because there is a one to one correspondence between a mea- 
sure on S 1 and its Fourier coefficients (see Shohat and Tamarkin [40]). Thus the existence 
of a solution of Eq. (l7. Ill) follows from the existence of a solution of Eq. d7.9l ) with the 
initial condition p = K$, of). Recall that Eq. d7.9l) is rewritten as the integro-ODE (17. 1 21) 
and (17. 1 31) . A proof of the existence of solutions of (17. 1 2b . (17. 1 31) for any t > is done 
by the standard iteration method. Since the proof needs a long calculation, it is left to 
Appendix B, in which we will also prove that x(t, s; 6, a>) is analytic in 8 and oj. Once a 
solution p, of (17.9b is obtained, a solution of (17. 1 II) is given through (17. 10b . Then, Prop. 7.2 
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is applied to show that solutions of Eq. (l7. Ill) are included in T . Note that when % is given 
as (17.381) . (17.1 II) becomes an autonomous system. Therefore, solutions define a semiflow 
on T. This implies that the dual Eq. fl7.30l) generates a semiflow <p t on Rx i(T). The proof 
of smoothness of Cp t will be done in Appendix B. It is also proved by the iteration method. 
■ 

The semiflow is also denoted componentwise as 

fr(z) = (zb, $(z), #?(z), ■ • • ), z = (zo, zi, ■ ■ ■ ). (7.43) 

By virtue of the variation-of-constant formula (see Eq. fl5.25l) ). <p\ proves to be of the form 

#"(zo, zi, ■ ■ ■ ) = (e r 'Tz ; - + gfeb, Zi, • • • ), J = 1, 2, • • • , (7.44) 

where g J t are nonlinear terms. Now we introduce another localization factor. Let 6 2 > 
be a sufficiently small positive number. By using a function x(t) above, we multiply 
the function x(dco(P s z, 0)/6 2 ) to the nonlinearity g J t , and define a perturbed map (f t = 
(zo,<p},(fi,- ■ ■) to be 

^(z ,Z!, • • • ) = (e^'Tzj + g/(z ,Z!, • • • ) •4 Joo( ^' Q) ), (7.45) 
for j = 1,2, • • • . Put 

^(zo,zi,---) = g f (zo,Zi,---)-^ -jr )• 

Fix a positive number r > 0, and put 

L = (id*, (e T n x , (e T2T T, ■■■), g = (0, g\, g 2 T ,--- ). 
Then, the time t map (p T of <p t is denoted as 

<p T : R x z'CF) -4Rx «'(T), ^ T (z) = Lz + g(z). (7.46) 

This is the desired localization of the semiflow of the original system (16.2b . By Prop. 7. 3, 
g J t is a C 1 mapping on R x i'CF). Since Joo(-,0) and P v are continuous on R x z'CF), 
g : Rxz'(F) — » Rxz'(F) is also continuous on Rxz'(F). Hence, the map ip T is continuously 
extended to the map on the closure R x i(T). Unfortunately, the distance rf«,(z, 0) is not 
C 1 in z. However, on the region such that d^PsZ, 0) < 5 2 or d^PsZ, 0) > 26 2 , g is a C 1 
mapping because;^(rfcoCP.sZ, 0)/5 2 ) becomes a constant. To see that g(z) ~ 0{z 2 ) as z —* 0, 
let us simply denote Eq. fl7.30l) by d\ Zj)/dt = T*\Zj) + Gj(z), where z = z{t) is given as 

Eq. fl7.39l) . Then, g/ in Eq. fl7.44l) is written as 

g/(Zo,Zi,---)= f (e T ^TGjiz{s))ds, 
Jo 

with the initial condition z(0) = (z ,Zi, ■ ■ •)■ Note that G is C°° and G(z) ~ 0(z 2 ) as z -> 0. 
Since solutions z(0 is C 1 in initial conditions, we ob tain G (z(s)) ~ O(z(0) 2 ) as z(0) — » 0. 
This shows that g/ is of order 0(z 2 ) as z — > on R x i(T) with respect to the metric d^. 
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When \\P cZ \\ Ec < S u \r](t)\ < 1 and ^(P^O) < 6 2 , then* (\\P c z\\ E JSi) = hx(\r](t)\) = 
1 and xidoo{P \z, 0)/6 2 ) = 1 . Thus Eg . (17 .461) coincides with the time r map of the semifiow 
of the original system O). When \\P c z\\ Ec > 26 1 or d^Ps Z, 0) > 2£ 2 , then^(||P c z|| £c /<5i) 
• X(d°°(P sZ, 0)/^) = 0. In this case, g = and Eq. d7.46l) is reduced to the linear map. 
Therefore, by taking 6\ and 62 sufficiently small, the Lipschitz constant of g 

t • / \ d BO (g(z),g(z')) 

Lip(g) := sup_ — — — - — (7.47) 

«'€Rxi(r) a oo\z, z ) 

can be assumed to be sufficiently small. 

Remark. We introduced the factors for localization in two steps. The one* (\\P c z\\ E J6i) 
is multiplied to the nonlinearity of the equation (17.301) . and the other x(doa(P s z, O)/^) is 
multiplied to the nonlinearity of the semifiow (17.441) . The reason is that if we multiply 
both of them to the equation (17.301) . then the proof of the existence of solutions for (17.301) 
(which is completed in Appendix B) becomes too difficult; x(\\P c z\\e c /8i) is essentially 
a finite dimensional perturbation, although x( doo (P S Z, O)/^) includes infinite dimensional 
terms P s z. On the other hand, if we multiply both of them to the nonlinearity of the 
semifiow of the original system (16.21) . then a center manifold of the resultant perturbed 
mapping does not coincide with a center manifold of the original system (16.21) because 
the perturbed mapping is not a semifiow for any differential equations in general (i.e. 
the property (p t+s = cp t o cp s is violated because of the perturbation for the semifiow), see 
Krisztin ||2"41 for details. However, if we introduce these factors in two steps as above, a 
local center manifold of the original system is correctly obtained as follows: In Sec. 7.4, 
we will prove the existence of a center manifold for the map (17.461) . We will show that 
if 5\ is sufficiently small, the center manifold is included in the "strip" doo(P s z,0) < 62. 
Since the map (17.461) is the same as (17.441) . which is a semifiow of the system (17.301) . when 
doa(P s z,0) < 6 2 , the obtained center manifold is a center manifold of the system (17.301) . 
When ||P c z|| £( . < 61, (17.301) is reduced to the original system (16.21) . Therefore, a local 
center manifold of (16.21) is obtained as a restriction of the center manifold of (17.301) to the 
region \\P c z\\ Ec < S x , see FigfJH 



7.4 Proof of the center manifold theorem 

Let us prove that the mapping ip T defined in Eq. d7.46l) has a center manifold, which gives 
a local center manifold for the original system (16.21) . The strategy of the proof is the same 
as that in Chen, Hale and Bin [jH, in which the existence of center manifolds is proved for 
mappings on Banach spaces. At first, we need the next lemma to treat the metric d^ as a 
norm. 

Lemma 7.4. For u = (uq, u\, ■ ■ ■ ) 6 R x z'CF), suppose that there exists a positive constant 
£3 such that \\uj\\* { < S3 for j = 1, 2, • • •. If constants {P(/?)}J =0 are sufficiently large, there 
exist positive constants A c and A s = A s (6 3 ) such that the inequalities 

doa(L m P c u , 0) < A c dcc(u, 0), m = 0, +1, +2, • • • , (7.48) 

and 

dL(L m P s u , 0) < A s e~ amr d x {u, 0), m = 0, 1 , 2, • • • , (7.49) 
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Fig. 12: A center manifold (the black curve) for the map (17.461) coincides with that of the 
semiflow (17.441) . The center manifold of the semiflow (17.441) coincides with that of the 
original system (16.21) in the region ||P c z||e c < 8\. 



hold, where a > is the constant appeared in Prop. 5. 19. 

Note that since P c is a projection to the finite dimensional vector space E c , the linear 
operator L restricted to E c is equivalent to an exponential of a matrix. Thus the inverse 
L exists on E r and L~ m P r is well-defined. 



Proof. For u = (w , * • * ) e R x i{T), L' n P c u is given by 

L m P c u = Oo , (e TlomT ) x U c ui , , , • • • ). 



Since u x e i(V) c i(W p (3,0)), Eq. d5779b is applied to yield 



T mTX _ v 1 ll(e r,omT ) x ncKil£„ " 1 L fi \\ Ul \\* fin 

* ((eU } n " Ml ' 0) = h i + ||(^)xn cMl |l^ " It * i + L,|| Ml |i;„ 



where = L 3 0j g. We can assume without loss of generality that Lp>\. Then, 
Eq. d5.74t provides 

1 H"lllo.» „ \H 1 W U ^}.n 



1 I I ^ 1 I I 

2" 1 + Hi L 



and 



l(Ml , ) = V ± " < V ± " M/3 ;" = d^m.O), 



A l INIL A l (20e )n 
rffl("i,0) = > — < > — = Q(P)do(u u O). 
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By using them, d((e Tl0 ' nT ) x Il c u i , 0) is calculated as 

d(( e TwmT fIi c u x , 0) = > — — -r 1 — C -LL-L 

/?=0 

CO 

-< E 



-P(fi)l+d /3 ((e T ^r rUcUlj0) 

1 LpdpiuuO) 
o P08)l+L^( M i,O) 



P^) l+db(«i,0)' 

where we choose a sequence {P(J3)}^ =Q so that £J=o LpQi/B) / P(J3) converges. Then, 

^ TO l ^ TO' ^ TO l + db(m,0) 

y^g© /y _J_ v-1 y 1 dfg(«l,0) 



TO v tif TO 7 U TO 1 + ^(«i,0) 



y3=0 /?=() j8=0 

= A C J( M1 ,0), A^^.^-L)" 1 . (7.50) 

Thus we obtain 

doo(L m P c u,0) < sup{u ,A c d(ui,0)} < sup{A c w , A c d(u\, 0)} = AcdooO, 0), 

where we suppose that A c > 1. Note that we did not use the condition ||wy||* y < 63 for 

Eq. CZ3S> - 

Next, L m P s u is given by 

L m J>,« = (0, (e riomr ) x n sMl , (e r2 '" T ) x «2, •••)• 
Eq. (l5.80l) is used to yield 



T v A 1 IK^Tn^di* 

< - 1 Ay— 1 INR, 

- ^ 2» 1 + M^—ll^ill^ " ^ ^ 2" 1 + e-^Wu^; 

where Mp := M 3 0(8 is assumed to be larger than 1. Note that the condition y < t> 3 
yields || M; -||* jB < Q(J3)6 3 . When \\uj\\^ n < Q(J3)6 3 , putting A' p {8 3 ) = (1 + Q(J3)6 3 ) provides 

1 , 1 

<A' p (5 3 y 



\ + e-^\\ Ul \\i n ~ ^ -l + Hmll^' 
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uniformly in m = 0, 1 , • • • . Therefore, we obtain 

d fi ((e Tlomr ) x U sUl ,0) < A'^Mpe-^dpiuuO). 

By the same way as in Eq. d7.50l) . we can verify that there exists a constant A s (5y) > 
such that 

d((e Twmr ) x IL s u u 0) < A s (6 3 )e- amT d(u u 0). 

In this calculation, constants P{fi) are chosen sufficiently large as before. Similarly, 
Eq. d5.81b shows that 

d((e T ' mT ) x Uj,0) < A s (S 3 )e- amT d(Uj,0), 

holds for every j = 2, 3, • • • . Note that the constant A s can be taken so that it is indepen- 
dent of j because the constant N c , a ,i3 in Eq. (l5.81l) is independent of j. Thus doo(L" l P s u, 0) 
satisfies Eq. d7.49b . ■ 

If \\PcU\\ Ec > 26 x or dUPsU, 0) > 26 2 , then g{u) = 0. Thus there exists D = D(6 U 6 2 ) 
such that the j-th component of g satisfies d(g(u)j, 0) < D for every j = 1, 2, • • • . This 
shows that there exists S3 = 63(61,62) such that ||g(w)y|lo { < 63 for all j. In what follows, 
we fix 63 in Lemma 7.4 so that g satisfies HgO^lloi ^ ^3 for all j. Then, Eq. d7.49b is 

applicable to g(u). Note that Eq. d7.48b holds for any R x i(T) without the assumption. 
Lemma 7.5. For a sequence {u ( ~ m) = (u ( ~ m \ u ( ~ m \ ■ ■ ■ )}~ =0 c R x i(T), suppose that 

L q P s u { - q) (7.51) 
as q —* 00, and that there exist constants C > and a sufficiently small b > such that 

d co (u ( ' m \0)<Ce bmT (7.52) 
holds for every m = 0, 1, • • • . Then, {w ( ~ m) }^ =0 satisfies 



uT"" = ip T (u { -' n - X) ) = Lu ( ~ m - l) + g(u { ~'"- l) ), m = 0, 1, 2, • • ■ , (7.53) 



if and only if it satisfies the equation 



u { - m) 



= L- m P c u -J] L k ~ m - l P c g(u ( - k) )+ ^ L k - m - l P s g(u { - k) ), m = 0, 1, 2, • • • . (7.54) 



k=\ k=m+l 



Eq. d7.53b means that {w ( m) }^ =0 = (w (0) , « ( x \ ■ ■ ■ ) is a negative semi-orbit of the mapping 
<p T . Eq. d7.54b is called the Lyapunov -Perron equation ||9l . 

Proof. Suppose that {u ( ~ m) \ satisfies Eq. d7.53b . By iterating Eq. d7.53b . we obtain 

M (0) = p cU (0) +PsU (0) 

= P c u {0) + LP s u ( ~ l) + P s g(u ( ~ l) ) 

= P c u {0) + L 2 P s u ( ~ 2) + LP s g(u { - 2) ) + P s g(u { - 1) ) 

q 

= P c u (0) + L«P s u™ + ^ L k ~ l P s g(u ( - k) ), q = 0, 1, 2, • • • . 

k=i 
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In a similar manner, we obtain 

m q 

= L- m P c u - J) L k " n ~ l P c g{u { ' k) ) + L^ m P s u ( - q) + £ L k - m ~ l P s g(u ( - k) ), (7.55) 



u { - m) 



k=\ k=m+\ 



for q = m, m + 1, • • • and m = 0, 1, 2, • • • . By the assumption, U rm P s u { ~ q) — » as g — > 0. 
Next thing to do is to show that EL +i L k ~ m ~ l P s g(u^~ k) ) converges as g — » oo. Eq. J7.49t 
is applicable to g(w (_m) ) to yield 

J oo (L*-'"- 1 P,g( M ( -* ) ),0) < A^- fl *- w - 1)T d oo (g(M ( - k) ),0) 

< LipC^A^-^-^^C^.O) 

< Up(g)A s Ce~ ia ~ h)kT e a(m+1)T , (7.56) 

which shows that L k ~" l ~ l P s g(u^ k) ) decays exponentially as k — » oo when a > b. Thus 
taking the limit q — » oo in Eq. J7.55t yields Eq. (l7.54t . 

Conversely, suppose that {w ( ~ m) } satisfies Eq. J7.54t . Because of the assumption Eq. J7.52t . 
the series J^=m+i L k ~ m ~ l P sgiu*-'®) exists. Replacing m by m+ 1 and using L for Eq. J7.54t . 
we obtain 

m+l oo 

La ( -" , - 1) = L~'"P c u m - V L k ~ m - l P c g(u ( - k) ) + V L k - m - { P s g{u ( ~ k) ). (7.57) 
Eq. J7.57t is put together with Eq. J7.54t to yield Eq. J7.53t . ■ 



Let X { m) (m = 0, 1, • • • ) be copies of the space R x Define a metric D m on Z ( m) 
to be 

D m (z, z) = e- hntT dM, z), u'eRx i(T), (7.58) 

with a small positive constant b. Let X be a subspace of the product Y\m=o X { ~ m) consisting 
of elements u = (w (0) , u ( ~ l) , w ( ~ 2) , • • • ) such that sup m D m (i/~ m \ 0) is finite. With the metric 
defined by 

D(u, v) = sup D m (u ( ' m) , v [ ~ m) ), (7.59) 

m>0 

X is a complete metric vector space (see Table 2). It is easy to verify the inequality 

D(u, +u) = D(u + u, 0) < D(u, 0) + D(w\ 0). (7.60) 
Motivated by Eq. J7.54t . let us define the map J : X x E c -» X to be 
f /(it, y) = (/<>, y), / ( - 1} ( M , y), / ( " 2) (w, y), • • • ), 

»; oo 

/ ( - m) (a,y) = L~ m y - J] L k ~ m - l P c g(u ( ~ k) ) + J] L k ~ m ~ l P s g(u ( - k) ). 

k=l k=m+l 

If the map /(• ,y) has a fixed point u = q(y) = (q {0 \y), <? (_1) (y), ■ • ■ ), qiy) is a solution of 
the Lyapunov-Perron equation (17.541) with P c uq = y. If g(y) satisfies conditions (17.5 lb and 
(17.521) . Lemma 7.5 shows that q(y) is a negative semi-orbit (that is, it satisfies Eq. J7.53t ) 
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for each y e E c . We will see that this qiy) gives a desired center manifold. At first, let us 
show that J is well-defined. 

Lemma 7.6. J is a map from IxE c into X. 

Proof. Let us show that D m (J^~ m \u,y),0) is bounded uniformly in m = 0, 1, • • • . It 
satisfies 

D m (J ( - m) (u,y),0) = e- bmT dM ( ~" l) (u,y),0) 

m co 

< e~ bmT dUL~'"y, 0) + e^Yj dUL k " n ~ l P c g{u ( ' k) ), 0) + e~ bmr £ ^(L^P^u™), 0) 

Eq. (l7.48l) shows that the first term e~ hTm d ca (L m y, 0) is bounded uniformly in m = 0, 1, • • • . 
Similarly, we obtain 

do O (L*- m - 1 P c g(ii ( -* ) ),0) < AcdooCgCu^)^) < Uv{g)A c dM- k \Q). 
Since u e X, there is a constant C > such that d^iu^'®, 0) < Ce Mr . Therefore, 

m m br _ bT(l-m) 

e- bmT J] dUL k - m ' l P c g(u { - k \ 0) < Up(g)A c Ce- bmT J] e bkT < Lip(g)A c C — — 

k=i k=\ e 

is bounded uniformly in m = 0, 1, • • • . The last term in Eq. fl7.62l) obviously tends to zero 
as m — » oo. This proves that J{u,y) e X. ■ 

Proposition 7.7. If the constants 6i and 6 2 are sufficiently small, J is a contraction map 
on X for each y e E c . 

Proof. For w, w' 6 X, we have 



D m (J ( - m \u, y), J ( -' n \u', y)) = e^dM^iu, y) - /""V, y), 0) 

Ps(g(u { ~ k) ) ~ g(,u' ( ~ k 

(7.63) 



< e- imT £ dUL^PMu^) - g{u' ( - k) )\ 0) + e- femr £ ^(L^PM^) ~ g(u' ( - k) )), 0). 



Eqs. fl7.48D and (17.471) are used to yield 

d x (L k - m - l P c (g(u { - k) )-g{u K - k) )),Q) < Up(g)A c d ca (u i ' k \u' { - k) ). 
Similarly, we obtain 

d^L^Psigiu^) -g( M ' ( " fc) )),0) < Lip(g)A^- a(fe -" ! - 1)T rf 00 ( M ( -^, H ' ( - fc) ). 
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Therefore, we obtain 

D m (J ( ~ m \u,y\J { ' m \u\y)) 

m oo 

< e- bmr ^ Up(g)A c dUu { ~ k) , u K ~ k) ) + e~ hmT Y U V {g)A s e~ a(k ~" l ~ l)T dM~ k \ u' { ~ k) ) 



k=\ k=m+l 
m 



< Up(g)A c Ye b(k ~ myT D k (u^ k \u' ( - k) ) + Up{g)A s Y e^'^ e Kk ~ m ^B k {u { - k \ u'^) 



< Lip(g) 



This yields 



k=\ k=m+l 
m oo \ 

A c Y e b(k -' n)T + A s Y e - a( - k - m - 1)T e Kk - m)T ■ D(u, u'). 

V k=\ k=m+l ) 



D(J(u,y),J(u',y)) = supD m (J ( ~ m) (u,y),J ( ~ m) (u',y)) 

m>0 

We can take 8\ and 6 2 sufficiently small so that Lip(g) becomes sufficiently small and 

/ £ br e br 



Lipfe) + A '*W^)7 ) < 1 ( 7 - 64 ) 

holds. This implies that /( • ,y) is a contraction map on X for each y e E c . ■ 

Remark. The numbers a and b are the same as those in Thm.7.1. The reason we in- 
troduced metrics D m and D is that the center manifold is characterized by the "slow" 



dynamics whose Lyapunov exponent is smaller than b, see Eq. d7.3l) . The above condi- 
tion for Lip(g) shows that if we take b > sufficiently small, Lip(g) (and thus Si and <5 2 ) 
should be small accordingly. Since the open set U in Thm.7.1, in which we can prove the 
existence of the local center manifold, is determined by 6\ and 5 2 , U also becomes small 
as a result. 



(7.65) 



By the contraction principle, J( ■ , y) has a unique fixed point u = qiy) on X: 

' qiy) = (q i0) (y), q^Xy), q { ~ 2) (y), ■■■), q ( ~ m) ■ k -> 

m oo 

q^(y) = L- m y-YL k - m - l Pcg(q ( - k) (y))+ Y L M P,gtf-*<y))- 

k=\ k=m+\ 

In particular, q {0) defines a map from E c into X (0) = R x i(f) given by 

oo 

q i0 \y) = y + Y L k ~ l P*s(q ( ~ k) (y))- (7.66) 

Since q(y) e X, there exists C > such that D m (q ( ' m) (y), 0) = e' hmT dM' m \y\ 0) < C, 
which verifies the condition (17.521) . Further, Eq. (l7.65l) shows that 



L"'P s q^(y) = Y L k - l P s g(q ( ~ k \y)). 



k=m+\ 
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Since this is a convergent series, L m P s q { ~ ,n) (y) —* as m — > oo, which verifies (17.511) . 
Therefore, Lemma 7.5 is applicable to conclude that {<7 (_m) (y)}," =0 is a negative semi-orbit 
for each y. 

Proposition 7.8. For any m = 0, 1 , • • • , 

(i) q ( - m \0) = 0. 

(ii) q ( ~ m) : E c — > is Lipschitz continuous. 

(iii) if llyll^> 2<J!, then q^ m \y) = L~ m y e E c . 

(iv) q ( -' n) : E c -» is a C 1 mapping. In particular, <? (0) :E c ^Rx i(T) is C 1 . 

Proof, (i) Since g(0) = 0, ^ ( -' M) (0) = satisfies Eq. d735b . 

(ii) For y,y e E c , we estimate D m {q^ m) {y), g (_m) (y')). By the same calculation as the proof 
of Prop.7.7, we obtain 

D(q(y),q(y')) < ; -A — -d^y'). (7.67) 



1 - Lip(g) \A C -^ — - + A 



e 



1 1 - e (b - a)7 



This means that q : E c — » X is Lipschitz continuous. In particular, we obtain 

A e bmT 

dU.q { - m) (y),q ( -' n \y')) < -, f T ^ — -djy,j\ 

1 -Lip(g)(A c -r— — -+A S - 



e hT - 1 1 - e< /;, - fl > r 
which proves the proposition. 

(iii) Put y = (y , yi, 0, • • • ), and y\ = £" =0 «/| ///) e E c . Then, the assumption implies 

\\y\\ Ec = \yo\ 2 + \ao\ 2 + \a l \ z + --- + \a n \ 2 >26 1 . 
On the other hand, Lr m y is given by L~ m y = (y , (e~ Tlomr ) x y u O, ■ ■ • ), and 

m m 
7=0 7=0 

Since Aj G V-TR, it turns out that IIL^yH^ = ||v|| £( . > 25\. By the construction of 
the nonlinearity g, g{L~ m y) = if ||L~" I .y||£ c > 25i. Therefore, q { ~ m \y) = Lr m y satisfies 
Eq. d7^5l) . 

(iv) For y, y* G E c and k e R, put 

4 ( ~ m) (y, /, *) = -(^ ( - m) (y + k/) - q { - m \y)). (7.69) 

K 

Then, it satisfies the equation 

1 

§ ( - m) (y,y*,^) = L~ m y* - - YL k - m - x PXg(q { -%) + Kt%,y\K))-g^-%))) 

+- y L i -'"- 1 P. v (^ ) (y) + ^ ) (y,y*,/c))-g(^ ) (y))), (7.70) 



k=m+l 
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for k # 0. If q { m) (y,y*, 0) exists, it should satisfy 



t m) (y,f,o) = L- m f-Y,L k - m - l pAq ( -%))t k) (y,y\o) 



dx 



oo , 



£=m+l 



(7.71) 



Motivated by these equations, we define a map J' : X x E c x E c x R — > X to be 7' = 

L~ m y* - - V L M PjM-*><y) + ku™) - sfo^OO)) 

1 oo 

+- y L M P t (gi^(y) + an™) - g(^(y))), « * 0, 



k=m+ 1 



IT*/ - y L k ~ m - l P c ^-{q ( - k) {y))u { - k) 

h dx 

dg 



k=m+l 



(7.72) 

We can prove that J' is a contraction map on X for each y, y* and /e by the completely same 
way as the proofs of Lemma 7.6 and Proposition 7.7. Hence, there uniquely exists u ( ~ m) = 
q ( ~ m) (y,y*,K) satisfying Eq. dTTOl) and Eq. jTJp . Taking the limit a- -> in Eq. d7770T) 
yields 

lim^ ( - m) (y,/,/c) = L-'Y - y L^»- 1 P c /(^-*>(y))lim§(-«(y,/,ic) 



i = l 



oo j 

+ V L^^/^-^lim^Cy./.ic). (7.73) 



£=m+l 



This implies that lmv->o# ( m) (y,y*,/c) is a solution of Eq. (l7.7 II) . By the uniqueness of a 
solution, we obtain 

\imq ( - m) (y,y*,K) = lim - (q ( ~ m \y + Ky*) - a H,,) (y)) = q ( ~ m \y,y* ,0). {1.1 A) 

From Eq. (l7.71l) . it turns out that q^~ m \y,y*, 0) is linear in y*. Thus we denote it as 

q ( - m) (y,y*,0) = dq { - m \y)f. (7.75) 

Then, dq { ~' n) (y) : E c — » x (_m) defines a linear operator for each y e E c . The remain- 
ing task is to show that dq^ m) : E c x E c — » x ( ~ m) is continuous. This is done in the 
same way as the proof of part (ii) of the proposition. For y',y'* e E c , we estimate 
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dcvicf m) (y,y* , 0),g ( m \y',y'*, 0)). Then, we can show that cf m) is Lipschitz continuous 
in y and y*. The details are left to the reader. This means that dg (_m) (y) gives the deriva- 
tive of q ( ~ m) at y e E c . ■ 

Now we define the center manifold W c of the map (p T by 

W c = {q (0) (y)=y + q(y)\yei c }, (7.76) 

where 

oo 

qiy) = J] L k - l P s g{q(- k \y)) e E^. (7.77) 
k=i 

Proposition 7.9. (i) W c is a dim-E c dimensional C 1 manifold, which is tangent to the 

space E c . In particular, q (0) (y) is expanded as q (0) (y) = y + 0(y 2 ) as y — » 0. 

(ii) V7 f is (p T invariant; that is, if T (W c ) c W c . 

(iii) For any £ € W c , there exists a negative semi-orbit {w ( ~ m) }^ =0 c W c satisfying uq = £o 
and 

dM~ m \Q) < Ce bm \ 
where b > as above and C is a positive constant. 

(iv) if <5i > is sufficiently small, the center manifold W c is included in the strip region 
{z 6 R x W) I cL(P s z, 0)<6 2 } (see FigH2>. 

Proof, (i) Since ^(O) = and q { ~ k) (y) is C 1 , # ( - fe) (y) is expanded as q^ k) (y) ~ 0(y). 
This shows that qiy) ~ 0(y 2 ) because g(z) ~ 0(z 2 ) as z — » (see Sec.7.4). 

(ii) Recall that {q { ~ m Xy)}™ =Q is a negative semi-orbit satisfying Eqs. (l7.51l) and (17.521) . De- 
fine g (1) (y) := ^> T (g (0) (y))- Obviously {<7 ( ~ m+1) (y)}~ =0 is also a negative semi-orbit satisfying 
(17.511) and (17.521) with some C > 0. Then, Lemma 7.5 implies that {<? (_m+1) (y)}^ =0 is a so- 
lution of the Lyapunov-Perron equation (17.541) . By the uniqueness of a solution, there 

exists y e E c such that q { ~ m+Y) (y) = q ( ~ m) (y') for m = 0, 1, • • • . In particular, we obtain 
<Pr(<l {0 \y)) = tpM~ i] (y')) = q <0 \y'), which proves <p T (q (0) (y)) e W c . 

(iii) This is obvious from the definition: if £ = g (0) (y), {q^~ m) (y)}™ =0 is a negative semi- 
orbit included in W c . 

(iv) Prop.7.8 (iii) implies that P s q (0) (y) = if \\y\\ Ec > 25 x . Thus sup & doo(.P s qf- 0) (y),0) 
becomes sufficiently small if <5i is sufficiently small. ■ 

If restricted to a small neighborhood of the origin, W c gives the desired local center 
manifold for Eq. (l6.2l) . 

Proof of Theorem 7.1. If Si > is sufficiently small, W c is included in the region 
{z 6 R x | d^P^, 0) < 8 2 ), on which x(d oo (P s z,0)/6 2 ) = 1. Thus ^-invariant 
manifold W c ' is also invariant under the map <p T given by Eq. d7.44l) . which is a time r 
map of the semiflow of the system (17.301) . Take uq e W c . By Prop.7.9 (iii), there is a 
negative semi-orbit {u^ m) }™ =0 c W c of ^ T satisfying Eqs. d7.51l) and (17.521) . Since <p T is a 
semiflow, we have 

ip T o tp t (u { - m) ) = <p t {ipM- m) )) = &(u ( ~ m+l) ), (7.78) 
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for each t > 0. This means that {<p t (u { m) )}~ =0 is a negative semi-orbit of (p T . Since (p t is C 1 
with respect to the metric dec, there is a positive number C such that 

cU^<y~ m) ), 0) < CdUu^, 0) < CCe bmT . (7.79) 

Further, L m P s (p t (u^ m) ) is estimated as 

L m P s $ t {u { - m) ) = L m+y P s u ( - m) + L m P s g t (u ( ~ m) ). 

Since M ( ~ m) satisfies (17.511) . L' n+l P s u { ~ m) tends to zero as m — » oo. By the same calculation 
as Eq. (l7.56l) . we see that the second term L m P s g t (u ( '~ m) ) also tends to zero as m — » oo. 
This shows that {&(w < ~ m) )}!L satisfies the conditions (17.511) and (17.521) . Therefore, it is a 
solution of the Lyapunov-Perron equation (17.541) . By the uniqueness of a solution, there 

is y 6 E c such that i£ f (w (0) ) = <7 (0) (y) e V7 C , which proves that W° is ^-invariant. 

In Eq. d7.30t . since e is a constant which is independent of t, W c (s) := W c n({e} xi{T)) 
is also ^-invariant for each e. 

On the region U = {z\\\P c z\\ Ec < S u WOI < 1}, ^(ll^lkMi) -ATd^OD = 1 and 
Eq. d7.30t is reduced to the original system (16.21) . Thus W c (s) D U is invariant under the 
semiflow generated by (16.21) . which gives a local center manifold stated in Thm.7.1 with 
U = U n ({e} x 1(T)) and = W c (e) n U. Parts (I) and (II) in Thm.7.1 immediately 
follows from Prop.7.9. It remains to show the part (III) of Thm.7.1. This is proved in the 
same way as Chen, Hale and Tan [9]. In [9], the existence of invariant foliations is proved 
for dynamical systems on Banach spaces. Though our phase space is not a Banach 
space, the distance from the origin Joo(z, 0) plays the same role as a norm. Thus with the 
aid of the estimates (17.48b and (17.491) . we can prove the existence of invariant foliations 
by the same way as O- The details are left to the reader. ■ 

7.5 Reduction to the center manifold 

Let us derive the dynamics on the center manifold and prove the Kuramoto's conjec- 
ture. Recall that for the continuous limit d2.lt of the Kuramoto model, Putting Zj(t, oS) = 

f^ 7T e j ^ e p t (8, af)d8 yields the system of equations (13.21) and (13.31) . Since solutions are in- 
eluded in V^o c F p (Thm.5.10 (iv)), the canonical inclusion is applied to rewrite Eq. d3.2t 
and (13 .31) as equations of the form (16.2b defined on F' p . The order parameter rj is defined 

as T](t) = (Zi,Pq) = (^o I Zi). For this system, we have proved that when < K < K c , the 
trivial solution (de- synchronous state) is asymptotically stable because of the existence of 
resonance poles on the left half plane. In particular, rj(t) — » as t — » oo. When K > K c , 
we have proved that the trivial solution is unstable because of the existence of eigenvalues 
on the right half plane. Thus a bifurcation from the trivial solution may occur at K = K c . 
In Sec. 7.1 to Sec.7.4, we have proved that there exists a smooth local center manifold near 
the origin in YlkLi F'p ^ K ls sufficiently close to K c . Our purpose is to obtain a differential 
equation describing the dynamics on the center manifold to reveal a bifurcation structure 
of the Kuramoto model. 

To derive the dynamics on the center manifold, we impose the following hypothesis 
as is assumed in the Kuramoto's conjecture. 
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(A6) The probability density function g(co) is an even and unimodal function such that 
g"(0) * 0. 

Note that g"(0) < when g(co) is unimodal. Thus (A6) implies g"(0) < 0. Under 
this assumption, there exists only one resonance pole A = on the imaginary axis when 
K = K c . Indeed, the imaginary parts of resonance poles on the imaginary axis are given 
as roots of the second equation of (13.171) . which is calculated as 

( co — y I co 

= lim n ~TT1 7l8(u)dco = hm — — - (g(y + co) - g(y - to)) dco. 

■ v ^ +0 Jr x 2 + (co - yY x-y+0 J X 1 + CO 1 

Since g is even, y = is a root of this equation. Sincegis unimodal, g(y + co)- g(y-co) > 
when y < 0, co > and g(y + co) - g(y - co) < when y > 0, w > 0. Thus y = is a unique 
root. This proves that A = is the unique resonance pole on the imaginary axis when 
K = K c . Further, the resonance pole is a pole of first order. Indeed, 



lim — 

*->+() ay 



I co I Ico cig 

— — - (g(y + co)- g(y - co)) dco = lim — — -—(co)dco. 

y= Jo x 1 + co z x^+q J x 1 + co 1 dco 



Since g is unimodal, we have dg/dco < 0, and the above value is not zero. This means 
that y = is a single root. Thus the center subspace E c is of one dimensional. Let | fi ) be 
the generalized eigenfunction associated with /t = 0. By the definition, |/i ) is given by 

(4>\Mo)= Hm f l -^<P(co)g(co)dco. (7.80) 

■ x ^ +0 Jrx- y-lco 

This is also written as 

I A)) = lim 



^=r-)= lim /( 



where the limit is taken with respect to the weak dual topology on F' p . The main theorem 
in this section, which confirms the Kuramoto's conjecture, is stated as follows: 

Theorem 7.10. For the continuous model (12.11) of the Kuramoto model, suppose that the 
function g(co) satisfies assumptions (Al) to (A3),(A5) and (A6). Then, there exist positive 
constants sq and 5 such that if K c < K < K c + s and if the initial condition h(6) satisfies 

e < ~ lie h(6)d6 <8 (7.82) 

for j = 1, 2, • • • , then the order parameter rj(t) tends to the constant expressed as 



r(t) = \rj(t)\ = J *f ^K-K c + 0(K - K c ), (7.83) 
^ nK^g"(0) 

as t —* oo. In particular, the bifurcation diagram of the order parameter is given as Figj2] 
(a). 
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K 



Fig. 13: Eq. (l7.83l) holds even if the shape of g(oj) is slightly deformed from an even and 
unimodal function. If g"(0) > 0, the unstable synchronous solution appears backward. 



Remark. A few remarks are in order. When g(a>) is a rational function with the assump- 
tion (A6), the theorem also holds without any modifications, though the existence of the 
center manifold is trivial because i(H+) is of finite dimensional. For both cases, we do not 
need the assumption (A4) because (A6) implies that the resonance pole on the imaginary 
axis is a pole of first order. The assumption (A6) is also not essential. The theorem is 
true as long as only one resonance pole gets across the imaginary axis at the origin. For 
example, let us slightly deform an even and unimodal g(co) so that it sags in the center as 
is shown in Fig(T3j If the sag is sufficiently small, only one resonance pole gets across the 
imaginary axis along the real axis, and Eq. (l7.83l) holds. In this case, since g"(0) > 0, |?7(f)| 
in Eq. d7.83l) is positive when K < K c . This solution is unstable as is verified from the 
dynamics on the center manifold given as (I7.99I ). This means that a subcritical bifurcation 
occurs and the bifurcation diagram shown in Figj2](b) is (partially) obtained. In general, 
the more complex a shape of g(co) becomes, the more complex bifurcations occur. The 
result (17.831) coincides with formal discussions and numerical simulations obtained by 
many physicists, see Daido lfT5l and Strogatz Il42l . for example. 

Proof of Thm.7.10. Suppose that an initial condition h(9) satisfies Eq. fl7.82l) . Then, we 
have 



for every j,/3 and n. Thus we can take 6 sufficiently small so that the initial condition 
(|Zi(0, ■)), |Z 2 (0, •)), • • •) for Eq. fl6.2l) is included in the neighborhood U (with respect 
to the metric d^) of the origin given in Thm.7.1. Then, the center manifold theorem is 
applicable. Let us derive the dynamics on the center manifold. 

Since we are interested in a bifurcation at K = K c , put s = K - K c and divide the 




< 



S-\\Po\{ 
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operator T { as 

T^(co) = T 10 <p((o) + |<P I ®Po(<o), (7.84) 

where 

T 10 cf>(co) = V^lW(w) + y (P Q I 0>P o (w). (7.85) 

Then, the operator 7io has a resonance pole at the origin and all other resonance poles are 
on the left half plane. Eq. (l6.2l) is rewritten as 

j t \Z,) = rfolZx) + |<Pol^i>l^o> - f <^o|Zi>|Z 2 >. (7.86) 

To obtain the dynamics on the center manifold, by using the spectral decomposition 
(15.651) . we put 

\Z 1 ) = ^a(t)\^ ) + \Y l ), (7.87) 

where |//o) is defined by Eq. fl7.81l) . | Fi) is included in the complement of E c , and 
where 

a(t) = i<z7^). (7.88) 
Do 

We will derive the dynamics of a. Since (P | fi n ) = 2/K c for n = 0, 1, • • • by the definitions 
of | //„) and the resonance poles, we obtain 



(P \Z 1 ) = a(t) + (P \Y 1 ), (7.89) 

and | P ) is decomposed as 

IA)> = ^lju > + |Io>, (7.90) 
Do 

where | Y ) e E^. By Thm.7.1 (I), on the local center manifold, we can suppose that 

{4>\Y X \ (cf>\Zj)~0(a 2 ,as,s 2 ), (7.91) 

for j = 2, 3, • • • and for every e F p . Let us calculate the expression of the center 
manifold. Substituting Eqs. (l7TH9l) . (f7T87T) into Eq.flO) for j = 2 yields 

^| Z 2 > = T 2 X | Z 2 ) + Z ((a + (PoWi)) (y ar| ju ) + | 7i>) - (a + </><> I Z 3 >) . (7.92) 

We suppose that da/dt ~ 0(a 2 ,ae,e 2 ), which will be justified later. Then, the above 
equation yields 

KK 

T*\Z 2 ) = — j-a 2 \no) + 0(a 3 ,a 2 s,as 2 ,s 3 ). (7.93) 

Lemma 7.11. Define the operator (T£ ) _1 : i(F p ) —* F' to be 

<0|(r 2 x )-^) = ~ lim f l -—cp(aj)Haj)g(co)daj. (7.94) 

2x^+0 J RX - ^p\co 
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Then, 

(T^Xr!)- 1 ^) = {TZT\TZ)W) = \f) (7.95) 
for any | iff) e i(F p ), and it is continuous on i(V). 
Proof. The straightforward calculation shows that 

0(<x>)</f(<x>)g(<^)<^> 



-I lim f 
2 J RX - V^o; 

J(f)(u))iJ/(u))g(u))du) - lim I * (p(oS)ifj(oS)g(co)du). 
R x ^ +0 Jrx - f4w 



Since the limit 



lim f \^^^{co)\]/{co)g{io)dco = {<p ■ iff\fi ) 

exists, the second term in the right hand side above is zero. Thus we obtain 

<0l(r 2 x xr 2 x r 1 ^) = (0l^). 

In the same way, 

{<f>\(T*)-\T*jf) = (<fi\(T*T\-2>/^(oW) 

2V-Fw 



1 r 2V-iw 

= -- lim —-(p{(o)iff(oj)g(oj)doj 

2 Jrx- ^hAoj 



Note that the right hand side of (17.941 ) is also written as —(<f> ■ iff\/j.o)/2, which coincides 
with the coefficient of n | <p ■ iff), up to the constant factor. Thus the proof of continuity of 
(T% X 1 on i(V) is the same as that of the continuity of Ilo on i(V). ■ 

Since (T? X 1 is continuous on i(V), its domain is continuously extended to the closure 



i(V). Since | yu ) e i(V) (see Prop. 5. 20) and it is given as Eq. d7.81l) . (J 1 *) 1 1 fj. ) is calculated 
as 



<0|(r 2 x rW= lim(0 (T^- 1 l —-) = -l lim f \——fta))g(a))da>. 

(7.96) 

Then, Eq. fl7.93l ) provides 

{(f> | Z 2 ) = — -a 2 lim — (h(oj)g(co)doj + <9(a 3 , or 2 e, o-e 2 , e 3 ), (7.97) 

4 ^+°Jr(x- V 3 !^) 2 
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which gives the expression of the center manifold to the | Z 2 ) direction. The projection of 
it to the center subspace is given as 



n |z 2 > = 



4^<Z 2 |jU > • lj"o) 
2D 



lim(z 2 l -=r-) • |ju ) = lim( Z 2 )-|/j >, 



where 



lim( l -— - Z 2 ) = 



KK C _ 2 



4 



im f ■— g{aS)da> + 0{a 3 ,a 2 s, as 2 , s 3 ) 



KK, 



c—2 

a 



(x - V-M 3 

i im r' 

*- +0 Jrjc- V 3 ! 



g"{io)d(x) + 0(a 3 , a 2 s, as 2 , s 3 ) 



KK, 



c—2 



a A ■ ng"(0) + 0(a 3 , a 2 s, as 2 , s 3 ). 



Thus we obtain 

KK 2 

n | Z 2 ) = =^a 2 ■ ng"(0) -\u ) + 0(a 3 , a 2 s, as 2 , s 3 ). 

16D 

Finally, the projection of Eq. fl7.86l) to the center subspace is given by 



(7.98) 



d„ - 



K . - 



^n | Zi> = r* n | z x > + |<A, I z^i p > - '-^{Po 

By using Eqs. flT871) . flT89l . flT90l) and dT98l) . we obtain 
4^li"o> = : y«r i x | A) ) + |(a + (^|T;))i|//o) 



isT / ng"{Q)KKl , , , , , 

--(a + <P I ?i» • - - c <* 2 |//o) + <V, <* 2 £, ae 2 , s 3 ) 
2 \ 16Z) 

e ng"(0)K 4 c 2 2 , , , 

^=a\fio) H = — or|ar| |/io) + 0(ear ,s a,s ,a ), 

2D 32D 



which yields the dynamics on the center manifold as 



d a 
—a = 



dt Do K c 



s + ng " { ^ K " \ a \ 2 ) + 0(sa 2 , s 2 a, s 3 , a 4 ). 
16 



(7.99) 



Since g"(0) < because of the assumption (A6), this equation has a fixed point expressed 
as Eq. fl7.83l) when s = K - K c > 0. Note that the order parameter 77(f) = (Zi,P ) is 
rewritten as 

17(f) = (ZuPo) = (P Q , Z0 = <P I Zi> = y a<A> I A>> + (A) I 7i> = a + 0(a 2 , as, s 2 ). 

(7.100) 
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Thus the dynamics of the order parameter is also given by Eq. (|7.99b . To prove that the 
fixed point (17.831) is asymptotically stable, it is sufficient to show the following. 

Lemma 7.12. D > 0. 

Proof Put 

f(A) = 1 - ^ f l -—g(a>)daj - nK c g(- V=U). 

1 JrA - y-la) 

By the definition of D , 

D =/'(0) = lim^ f ^-g(w)Jw+ V^T^'CO) 



lim f l -^^g'(co)dco+ V^T^ c e'(0). 



Since g(co) is even, 

D = -- ^ lim I — — ,g'((o)d(o = -K c ]im I — — -g'(oj)da>. 

Since g(aj) is unimodal, g'(oj) < when to > 0, which proves that D > 0. ■ 

Since D > 0,K C > 0, g"(0) < 0, the fixed point a = (de-synchronous state) is 
unstable and the fixed point Eq. fl7.83l) (synchronous state) is asymptotically stable when 
s = K - K c > 0. This completes the proof. ■ 



A Resonance poles of the Gaussian distribution 

In this appendix, we give a few simple properties of resonance poles for the Gaussian 
distribution g(co) = e~^~ 12 / yZn. In particular, we prove that it satisfies the assumption 
(A4) (the other assumptions (Al) to (A3), (A5) and (A6) are trivially satisfied). 

The transition point is given by K c = 2 y/2/n. The equation (14.1 II) for obtaining the 
resonance poles is reduced to 




(A.l) 



Let Aq, Ai,-- • be roots of this equation with > Re(/t ) > Re^) > • • • . 

Proposition A.l. If A n is a resonance pole, so is its complex conjugate A n . 

Proof. This is always true when g(oS) is an even function, and immediately follows from 
Eq. fl4~TT1) . ■ 

Proposition A.2. There exist infinitely many resonance poles. As n — » oo, Re(/l„) — » -<x> 
and they approach to the rays arg(z) = 37r/4, 5n/4. 
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Proof. If \A\ is sufficiently large, Eq. dA.ll) is estimated as 

e A2/2 -^— + 0(l/A), (A.2) 

which gives the desired results. ■ 

Proposition A.3. When < K < K c , there exists only one resonance pole on the real 
axis. 

Proof. Put 

fW = M^-£e-*»dx\-j. (A3) 

If A is negative, 

/'(A) = Ae A2/2 1 y| - £e~* ,2 dx\ - 1< 0. (A.4) 

Thus /(/I) is monotonically decreasing on the negative real axis. If < K < K c , then 
/(0) = V^72 - 2/K < 0, and if A < is sufficiently small, f{A) > 0. Thus the 
intermediate-value theorem proves the proposition. ■ 

Proposition A.4. All roots of Eq. dA.ll) are simple roots. Thus the Gaussian distribution 
satisfies the assumption (A4). 

Proof. Put 

f(A) = f l T=r g(a>) + 2ng(- ^A) - | (A.5) 

JrA - y-loj K 

(this is the same function as Eq. dA.3l) ). The derivative is calculated as 

f'(A) = - f l ——g(a>)dio-2jr>Ei g '(->ETX) 

Jr(A- V^Tco) 2 

= - f g'(to)dco-2n^g'(- V-U). 

JrA - V-lw 

Since the Gaussian distribution satisfies g'(oS) = -cog(co), we obtain 

f(A) = f ~]^_ g(a>)dco + 2nAg(-^P\A) 
JrA - y-la> 

= A \ l -^^g(io)dco + 2nAg(- V-U) - 1. 

JrA - y^Ato 

If A n is a resonance pole, it satisfies f(A n ) = 0. By using this, the above is rewritten as 

f'(A n ) = A n j-\. (A.6) 
Since Re(/t„) < 0, f'(A n ) ^ 0. This proves the proposition. ■ 
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B Proof of the existence of solutions of Eq.(E3) 

In this Appendix B, for Eq. (|7.9b with^ given by (17.381) . we prove the existence of solu- 
tions and smoothness of them with respect to initial conditions to complete the proof of 
Prop.7.3. Recall that Eq.dTll is equivalent to the integro-ODE (TTTHl) and (T77T31) . Put 



K 



fi(x, 77, a) = ^(rje-^ + T}e^) -*(^) ' xMVW 
f 3 (x, 77, a) = ^(rje-^ + rje^ x ) ■ \K C + **(^) ^(MODj , 

where a = (otq, ct\, ■ ■ ■ ,a n ) is defined through Eq. (l7.41l) and (17.421) . Then, (17.121) and 
(17.131) are rewritten as a system of integral equations 



fi(x,rj,a) = — =(jie 

K, ^ ^ .ll^zll 



x(t,0;6,to) = tot + 6 + I fi(x(s, 0;9, to), tj(s), a(s))ds, 



f 2 (x(s, t; 6, to), tj(s), a(s))ds 



p t (t, to) = h(x(0, t; 9, to), to) exp 



2x X exp X 



f 2 (x(r, t; 6, to), t](t), a(r))dr 



Tj(t) = (Z l ,P, 

K 



0) = f Z,(t, 
Jr 



Mx(s, t; 9, to), r](s), a(s))ds, (B.l) 



to)g(to)dto, 



aj(t) = -=-{Z x \nj). 
2Dj 

In particular, Eq. (l7.17l) is rewritten as 
Z k (t, to) - 



■In 



e' rike p t (8,to)d6 = 

Jo 



In 



e^ m ^h(d,co)d6 



(B.2) 



+ 



1 /"»2/r r*t r r*0 

2^ J, 6 VrTa ' ( ''° ;e '"j o ex P J /2«r, 0; 0, w), 77(r), a(r))d 



fi(x(s, 0; 0, r/(s), a(s))dsd6. 



Since there is a one to one correspondence between a measure on 5 1 and its Fourier 
coefficients (see Shohat and Tamarkin RUlO . p t is uniquely determined by {Z k }™ =Q . 

Proposition B.l. For any signed measure h(9, to) on S parameterized by to e R satisfying 
Eq. d7.14l) . the system of integral equations (IB.ll) has a unique solution for any t > 0. 

Proof. At first, we construct a local solution by the Picard's iteration method. Once we 



obtain x, 77 and a, then Z l5 Z 2 , • • • are obtained by substituting them into Eq. (IB.2l) . Thus 
we consider to construct x, rj and a by the iteration. Define sequences {x m (t, 0; 6, to)}™ =Q , 
{Z mA (t,to)}™ =0 , {77m(OC =0 and {a m (t)}^ to be 



xo(?, 0; 6, to) = tot + 9, Z Q i(t,to) 
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e^~ w h{9, to)d6, 



(B.3) 



and 



x m+ i(t, 0; 6,co) = x + fi(x m (s, 0; 9, of), rj m (s), a m (s))ds, 



i (t, of) = f e VrTv '» +1 «' 0M h(9, 0))d9 + -!- f 
Jo 27tJ 

I' exp [f 



e -V c Tx ra+ i(t,0;(9,w) x 



h{x m {j, 0; 0, 77 m (r), a m (r))Jr 



Mx m (s, 0; 6>, oj), t^O), a m (s))dsd9, 



L 



7] m (t)= Z mA (t,o))g(oj)do), 



K 



a mJ (t) = -^(Z„ u a m (f) = («m,o(0, • • • ,a m ,„(ty). 

2Dj 

(B.4) 

In what follows, we suppose that 9, to e R. In this case, we can verify by induction 
that x m (t, 0; 9, of) £ R for every m. We will show later that x(t, 0;9, of) has an analytic 
continuation in 9 and oj. 

Take a positive number M\ and put 



R = {(x, 7], a) | \tj - 7] Q \ < Mi} c R x C x C 
When (jc, 77, or) e i?, there exists a positive constant M 2 such that 

\f i (x,? 1 ,a)\<K\r ] \<M 2 , 
for ?' = 1, 2, 3. If (x m -\, Tjm-i, a m -\) e R, \q m - 770 1 is estimated as 



n+l 



(B.5) 



(B.6) 



l/3(^m-i, J/m-i, a m -i)\g(o))dsd9da> 



-2n n t 



\T] m -T} \ < I I 

JrJo 

+ T~ I I f ex P I I^Un-i, J7 m -i,a m -i)|rfT 
^kJrJo Jo Mo 

< \ \ \x m (t, 0; 9, of)-9\- \h(9, oj)\g{of)dojd9 +— I e MlS M 2 g{of)dsd9doj, 

JrJo 27tJ r J Jo 

and |x m - 9\ is estimated as 

\x m -9\< \o)\t + \fi(x m -i,j] 

m—\i &m- 1 

)\ds < (M + M 2 )t. 

Jo 



Thus we obtain 



Vim - Ho\ < 



JrJo 



(M + Af 2 ) • |/z(#, cS)\g(uf)du)d0 ■ t + (e M2t - 1). 



By the assumption (A2) for g(of) and the assumption (F2) for h, the integral above ex- 
ists. Hence, there exists T\ > such that when < t < T\, \q m - 770I < M\ holds. 
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Now the induction shows that (x m ,Tj m ,a m ) e R and thus Eq. (IB.6l) is valid for every 
(x m , ?]„„ a m ), m = 0, 1, • • • on the interval < t < T\. 
When < t < T\ , we obtain 



\xi - x Q \ < l/iOo, ?7o, a )\ds < M 2 t 
|Zii-Zoi| < 



pin 

Jo 



— r r exp r i/ 2 (x , 770,^0)^ 

^ Jo Jo Jo 



+ 



|/ 3 (x , 770, a )\dsd6 



r>2n 

Jo 



2/r w 




e MlS M 2 dsd9 



Jo 



< (M + M 2 )r I \h(9,to)\d9+ — 

= C l (\to\ + M 2 )t + (e Mlt -1), 

where we put C\ = J^\h(6, to)\dQ. There exists a constant M 3 > such that e Mlt - 1 < M 3 f 
on < t < T { . Then, 

|Zi,i-Zo,il < (C 1 |w| + C 1 M 2 + M 3 )f, 

|?7i - ?7ol < ( (CiM + CiM 2 + M 3 )tg(co)dio := M 4 J. 
Jr 

Let /lj = V-Tv 7 (j = 0, • • • , n) be resonance poles on the imaginary axis. By the defini- 
tion, \a\j - a j\ is calculated as 



\a h j - a j\ < 



K 



< 



2\Dj 
K 

K 



|(Z u -Z ,i|^)| 



lim 

t->+0 



f ! 

Jrx + V-lCty - 



JLlim f 
2|D y |^+oJ R 



C\\co\ + C { M 2 + M 3 



-(Z u 0, of) - Z ,i(f, w))g(<y)<ia; 



g{of)do) ■ t := M 5J f. 



2|Z> ; | *-+o J R | x + V^Tfjy 
Putting M = max{M 2 , M 4 , M 5J } yields 

|*i - x l < Mt, I?/! - T/ol < Mr, jor^ - or j| < Mr 



(B.7) 



for < t < T\. Since fi,f 2 , / 3 are Lipschitz continuous, there exists a positive constant L 
such that 

/ n 



\fi(x, T], a) - fi(x', T], a)\ < L 



\x 



- x'\ + \r]-ri'\ + ^Yj\aj- a'j\ 



7=0 



i= 1,2,3, (B.8) 



for (x, i], a) 6 R and (V, rj', a') e 7?. Let us prove by induction that there exists a positive 
constant /V such that inequalities 



N' n L m (n + 3) m +1 

\x m +i — x m \, \rj m+ \ - T] m \, \a m+ i j — a m j\ < M — — — t 

(m + 1)! 



(B.9) 
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hold on < t < Ti for m = 1 , 2, • • • . Indeed, 
Jo 



< 



L \x m — X m -\\ + \i] m — 7/m-ll + / , \&m,j ~ &m-l,j\ 





7=0 

N m - l L'"' l (n + 3)"'- 1 



(is 



< L(n + 3)M 



ml 



■s m ds = M 



N m-l L m( n + 3 y 
(m+ 1)! 



+ m+ 1 



Next, 



l^m+1,1 Z m j | 
-2tt 



< 



+ 



X 

- r 

2tt Jo 



Lfcteiw ^m, «m)I^T • |/ 3 (x m , ?7 m , a m )\dsd9 



+ 



\h{x m ,ri m ,a m )\dsd6 



^ r*hx r-t /-»0 /->0 

—J J exp[ J / 2 (x m , rj m , ajdr] - exp[ J / 2 (x m _i, 77 w _i, Q? OT -iWr] 
^ J" J" exp[ J" |/ 2 (x m _i,77 OT _i,or m _i)|rfT] • |/ 3 (x OT , ?7 m , a m ) - f?,(x m -.uri m -i,a m -x)\dsd6. 



+ 

There exists A/6 > such that 



exp[ J /2(x m , ?7 m , a m )dr} - exp[ J / 2 (x m _i,?7 OT -i,a m -iVT] 

< M 6 I |/ 2 (x OT , ?7 m , a,„) - f 2 (x m -i, i] m -i, a m -\)\dT 
Jo 
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holds when < s < T { . This provides 



- *J • 1^(5, W)|rf0 + — J |-X m+ i - Xj 



J() 



2^J J J M2M6L \ X ™ ~ X ™-^ + fo™ ~ Vm-ll + ^IVmJ ~ 0! m -l,j 

2^X X e 



+ 



< M 



27T \nn-l Tm 



drdsdO 



1 

+ — 

2;r 



(m+ 1)! 

2tt /-v />j 



r +1 |/i(0,w)|d0 + Af 



(m+1)! 




M 2 M 6 L(n + 3)M - — r m drdsde 



3 jojo w! 

27T /-f )\7m-l jm-\( , oyn-l 

+ 3)M - g + J; S m e M ^dsd6 



< M 



^ f f «« 

2ttJo Jo 
N m - l L m (n + 3) 



m! 



1 



2tt 



|/i(6», <u)|</0 + (e M2f - 1) + ^ + g 



M 2 M 6 



,Mif 



(m+ l)! 

Take a number Ni(a>), which is independent of m, so that 



m + 2 



1^(0, o>)|d# + (e M2? - l) + -^t + e M * < Ni((o) 
o m + 2 

on the interval < t < 7\. By the assumption (F3) for /?, Ni((o) is bounded as ct> — » +00. 
Thus we obtain 



JN" 1 ' L m {n + 3)'" 
|z m+u - z^iig^jo* < m — - — — — r + 
r (m+1)! 

By the assumption (A2) for g(u>), N 2 := ^ s Ni{co)g{o))dco exists: 

N m ~ l U\n + 3)'" 

Wm+l - ?7ml < M — 

(m + 1)! 

Similarly, \a m+ \j - a m J is estimated as 

k r 1 



iVi ((jS)g(aS)da). 



r +1 • 7Y 2 . 



|Of m+ ij- &m,j\ — 



-A. lim r. 

2|D ; -|^+oJ R | 



2|D ; -| *-+o Jr^ + V-Uj; - w)| 

N m-l L m( n + 3 y 



N m ~ l r n (n + 1) m 
(m + 1)! 



M 



-r m+1 • /y 3 , 



(m+ 1)! 

with some constant N3 > 0. Now Eq. (IB.9l) is obtained by putting /V = max{ 1 , N 2 , N3}. 
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Eq. (|B.9D shows that sequences {x m }, {rj m } and {a m } converge to solutions x(t, 0; 6, to), 77(f) 
and a(t) uniformly on < t < T\, respectively. Substituting them into (IB .21) provides 
Z\ , Z 2 , • • • . They give a local solution of Eq. dB.ll) on the interval < t < T\ . 

Next thing to do is to show that the solution is extended to any t > 0. For this purpose, 
let us show that x(t, 0; 9, to), 77(f) and a(t) are bounded for any bounded t > if they exist. 

Suppose that there exists f > such that \rj(t)\ > 2 for t > f . Then,;^(|?7(f)|) = 0, and 
thus Eq. (l7.1 II) is reduced to the linear equation dZi/dt = T w Zi. This shows that when 
t > t , 77(f) is given by 77(f) = (e Tlo(t ~ to) Zi(t , to), P ), which is bounded for every f > f due 
to Prop. 5. 13. Thus 77(f) is bounded for every f > 0. In a similar manner, we can verify that 
a (t), ■■■ , or„(f) are bounded for every f > because if |ar ; | > 28\, then^(||P f z|| £ ,/^i) = 
and Eq. (l7. Ill) is reduced to the linear equation. If 77(f) and a(f) exist and are bounded, the 
right hand side of Eq. d7.12l) defines a time-dependent vector field on the compact space 
S l . Thus a solution x(t, 0;6, to) exists. Then, substitution of them into Eq. (IB.2l) shows 
that Z\, Z 2 , ■ ■ ■ are bounded for any bounded f > 0. 

We have proved that if solutions exist, they should be bounded. Now we can extend 
a local solution defined on < f < T x in usual way: Repeating the iteration method by 
using an initial condition given at f = T u we can extend a local solution to a larger interval 
< f < T 2 . Repeating this procedure, we obtain a sequence T 2 < T 3 < ■ ■ ■ such that a 
solution exists on < f < T n . Suppose that 7^ := lim,,^ T n is finite. Then, since the 
solution is bounded at f = T^, we again use the iteration method with the initial condition 
given at f = Too to extend the solution beyond T^. Therefore, the solution can be extended 
for every f > 0. This proves that Eq. dB.ll) has a solution for any f > 0. The uniqueness of 
solutions is proved by a standard argument using the Gronwall inequality and is omitted 
here. ■ 

Proposition B.2. A solution x(t, 0; 6, to) of Eq. d7.12l) has an analytic continuation with 
respect to 8 e C and to eC. 

Proof. Once 77(f) and a(t) are determined, (17. 1 21) is an ODE of x. Since the right hand 
side of d7.12l) is analytic in the parameter to and the dependent variable x, the proposition 
follows from the fundamental theorem of ODE. ■ 

At this time, Prop.7.2 is applicable: solutions of Eq. d7.1 II) exist in T , and thus solu- 
tions of Eq. d7.30l) exist in R x i(T) for any f > 0. 

Proposition B.3. Solutions of Eq. d7.30l) in R x i(T) are continuous in initial conditions 
(| Zi(0)), I Z 2 (0)), • • • ) with respect to the metric d^,. 

One of the modern ways to prove the continuity of solutions of ODEs is to use the implicit 
function theorem on a Banach space. Unfortunately, our space i(T) is not a Banach space 
and the implicit function theorem does not hold in general. Thus we use more "classical" 
way by using the iteration method. 

Proof. Take two elements Z(0) = (Z 1 (0),Z 2 (0), ■ ■ ■ ) and Z*(0) = (Z^(0),Z*(0), • • • ) in T. 
Let (| Zi(f)), I Z 2 (f)>, • • • ) and (| Z[(t)), \ Z' 2 (t)), ■ ■ ■ ) be solutions of (17301) in i(T) with the 
initial conditions 



(| Z,(0)>, I Z 2 (0)>, • • • ), and (| ^(0)) - k\ ZJ(0)> , | Z 2 (0)> - k\ Z*(0)>, • • • ), 
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respectively, where fteR and |Z fc (0)> = j'(Z a (0)),|Z*(0)) = z'(Z*(0)). We want to prove 
that if | Z* k (0)) — » with respect to the weak dual topology uniformly in k = 1, 2, • • • , then 
| Z'At)) — » | Z k (t)) uniformly in k = 1, 2, • • • . The constant k is introduced to prove the next 
Prop.B.4. In this proof, the reader can assume that k = 1. Since Z(0),Z*(0) e T , there 
exist signed measures h and h* such that 

r>2n r>2n 

Z t (0)= I e^ lk6 h{9,cS)d9, Z*(0) = I e^ lk6 h\9,co)d9, (B.IO) 
Jo Jo 

for /c = 1, 2, • • • . Put h'(9, 10) = h(9, co) - Kh*(9, co). Then, Z k (t) and Z' k (t) are Fourier 
coefficients of solutions p and p' of Eq. (|7.9D satisfying the initial conditions h and h', 
respectively, and i{Z k {t)) = \Z k {t)), i{Z' k {t)) = \Z' k {t)). Let (x, rj, a) and {x',r]',a') be 
corresponding variables; x and x' are solutions of the first equation of (IB.ll) . in which 
/z is replaced by h and /?', respectively; rj(f) = (Pq \ Z\), aj(t) = K/ilDj) ■ {Z\ \/J.j), and 
rf, a' are defined in a similar manner. Put 

y = x - x, Y k = Z k — Z' k , £ = rj — rj', yj = ctj - a). (B.l 1) 

Then, they satisfy the integral equations 

y(t,0;9,co) = I (f^x, rj, a) - fi(x - y, t] - a - y))ds, 
Jo 

pin p2n 

Y k (t,co) = e yrikx h(G,cS)d9- e^ lKx - i) h!{Q,oS)dQ 
Jo Jo 

+— I e^~ lkx I exp[ I f 2 (x, 77, a)dr]Mx, 77, a)dsd9 

- n J « ex P[ J M x -y>1 -£><*- 7)dr]f 3 (x -y,r]-Z,a- y)dsd9, 



In 

Z(t) = (Po\~Y l ), 
K 



y j (t) = —(Y l \fi j ). 
2Dj 

(B.12) 

Let us show by the iteration method that these equations have solutions which tend to zero 
as \Y k (0,co)) = *|Z*(0)> -» 0. Define sequences {y m }~_ , {Y m , k }~_ , and {rmX-o 

by 

/-•2;r 

v = 0, Y 0M = Z k (0) - Z' k (0) = KZ* k (0) = k e" rike h*(9, oj)d9, (B.13) 

Jo 
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and 



(/i(x, r], a) - fi(x - y m , 7] - (j m , a - y m )) ds, 



y m +i - I 

Jo 

pin pin 

Y m+U = I e^ kx h(0,u>)dO- \ e" nk(x - y '" +>) h'(6,co)de 
Jo Jo 

pin pt pQ 

+— J e^~ ikx j exp[J f 2 (x,Tj,a)dT\f 3 (x,7],a)dsde 
-^J e^ k{x - y "^J exp[ J f 2 (x - y m , rj-^ m ,a- y m )dT\f 3 (x - y m , rj-^ m ,a- y m )dsd6, 



K 



J m ,j(t) = -=-{Ym,\ 

2Dj 

(B.14) 

Since f„(0 = (A) I K ,i> = k(P q |ZJ(0)>, ^ as |Zf(0)> -> with respect to the 
weak dual topology. Since yoj(t) = K/lDj ■ /f<Z*(0) \fif) = /cII 7 |Z*(0)> and since Ily is 
continuous on i(V), y j — > as | Z*(0)) — » in z'(V). To estimate Eq. (IB.14l) . we use the 
same notation as the proof of Prop.B.l. For a small interval < t < T\, there exists a 
constant L > such that 



(B.15) 
(B.16) 



lyi-yol ^ \f\(.x,r],a)- fx(x-y ,ri-t;Q,a-y Q )\ds 



< l (|v | + iftl + Irol)^ = + Irol) ~ 0{Kt\ 



where \y \ = lro,ol + Iro.il + • • • + |y ,„|. This shows that 

lyi - yd = bil -» 0, 



< 



+ 



+ 



+ 



as | ZJ(0)> -» 0. Next, for any <p e F p , 

\(<P\Y lk -Y , k )\ 

k f f (e^ k(x ~ y>) - e^ ke )(P(to)h*(6,to)g(co)d8dco 
JrJo 

f f |e - e VZT ^-- 1 '' ) | • \(P((o)h(9, io)\g(co)d9dco 
JrJo 

^ XX Ie ^ ~ eV=T " ( " yi)| X exp 'X l/2(x ' ^ a)l i/3(x ' 775 a)l ' ^ oj)i8(aj)dsded( ° 

I/3 (x, 77, a-)| • \(p(u>)\g(co)dsd6dco 

2^X X X^X'^ 2 ^' ?1, a ^ dT \^ X ' T1, a ^~fo( x ~y°> 1 ~ Q,_ ^o)l l0(w)lg(w)d«/6Wa). 
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Let us show that \((f> \ Y\ k - y ,Jt)l ~* as I Z*(0)) — » uniformly in k = 1, 2 • • • , and it is 

of order 0(/<t)- Since yi,£o,7o — » as | Z*(0)) — » and since £o,7o ~ 000, yo ~ 0(Kt), 
it is easy to verify that the second term to the fifth term in the right hand side above tend 
to zero as | Z*(0)) — » uniformly in k = 1, 2, • • • , and they have an order O(Kt). Next, we 

estimate the first term A := f R ffe^ rTk{x ~ y,) - e VrT ^)0(w)/z*(6', co)g(to)d8dco as 

|A| < f f % - yi - 0| • \<p(oj)h*(9, io)g(uj)\d9daj 
JrJo 

< J J £|otf + J ?7, a)ds — y\\ • \(p(to)h*(8, u))g(aj)\d6du> 
JrJo Jo 

<kt\ I \a + M 2 \ ■ \<p{aj)h*{6, cS)g(co)\dedo) + k\y x \ I I \cf>(co)h*(6, co)g{o))\dBdco. (B.17) 
JrJo JrJo 

This and Eq. (IB.15D show A ~ <90O, and thus \{(p\Y u - Y 0<k )\ ~ 0(Kt). On the other 
hand, with the aid of the Fourier expansion of the factor e^ k<JC ~ y ^ in A, we can verify 
that A tends to zero uniformly in k = 1, 2 • • • as (cf> | Z*(0)) tends to zero uniformly for 
j= 1,2---; 

<0 1 Z*(0)> = [ [ e ^ je <f>(io)h*(6, co)g(co)d8dco -> 0, j = 1 , 2, • • • . (B. 1 8) 
JrJo 

These facts prove that there exists a number B = B(k, cp, Z* k (0)) such that 

\(<f>\Y hk -Yo, k )\ < B-Kt, 



and B — > as | Z£(0)) — » uniformly in fc. Therefore, £i - £ an d Tij _ To,y are also of 
order O(Kt) and they tend to zero as |Z*(0)) — » 0. With these estimates, we can show 
by induction that |y„ I+ i - y m \, \{(p \ Y m+Uk - Y mJc )\, |£ m+1 - and |y m+1J - y mJ | have order 
0(Kt n /n\), and they tend to zero as | Z£(0)) — > 0. Thus the sequences converge to solutions 
y, (0 1 F^), £ and y ; - of Eq. OB. 121) . and it turns out that the solution is of order O(k) and it 
tends to zero as | Z* k (0)) — > uniformly in k. The detailed calculation is the same as the 

proof of Prop.B.l. In particular, since (<f>\Z k (t) - Z' k {t)) — » as |Z*(0)) — » uniformly 
in A;, this result means that solutions of Eq. fl7.30l) are continuous in initial conditions with 
respect to the weak dual topology. As was remarked in Sec. 7.4, this also proves that so- 
lutions are continuous with respect to the metric dco- ■ 

The remaining task to complete the proof of Prop. 7. 3 is to show that solutions are 
differentiable with respect to initial conditions. 

Proposition B.4. Solutions of Eq. fl7.30l) in R x i(T) are Frechet differentiable in initial 
conditions with respect to the metric dco. 



Proof. We use the same notation as the proof of Prop. B. 3. Define two solutions (| Z\{t)), \ Z 2 (t)), • • • ) 
and (| Z[(t)), | Z' 2 (t)), ■ ■ ■ ) as the proof of Prop.B.3. We consider the limit k — » for fixed 
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(I Z*(0)), I Z*(0)), • • • )• In the proof, we have proved that the integral equation (IB. 121) has 
a solution, which is of O(k) as k — » 0. Therefore, the following quantities exist: 

ftf,Z(0),Z*(0)) := lim^ = lim ^— y(0,Z(0),Z*(0)) = 0, 

Y k (t,Z(0),P(P)) ■= lim - = lim ^— -A 7,(0,Z(0),Z*(0)) = Z*(0), 
,• f 



f (f, Z(0), Z*(0)) := lim ^ = lim ^ '-, £(0, Z(0), Z*(0)) = (P„ | Z*(0)>, 

c-»o k *-»o a: 

r y (f, Z(0), Z*(0)) := lim 22 = lim ^ (0 , Z(0), Z*(0)) = -^<Z*(0) |^>, 

J k->0 K k—>0 K 1 2D 

(B.19) 

where Z(0) = (Zj(0), Z 2 (0), • • • ) and Z*(0) = (ZJ(0), Z*(0), • • • )• This implies that | Z k {t)) is 
Gateaux differentiable and i(Y k ) = \ Y k ) is the Gateaux derivative. Let us show that i{Y k ) = 
| Y k ) is the Frechet derivative of | Z k (t)) with respect to the initial condition. Dividing by k 
and taking the limit k — » in Eq. dB. 121) yields 

y(*,Z(0),Z*(0)) = f V/ 1 (x,7 7 ,a)(j,|,r) T ^, (B.20) 
Jo 

and 

? t (f,Z(0),Z*(0)) 

nQ.li r>2n 

= I e^~ lkx ^ky(t,Z(0),Z*(0))h(6,a>)d6 + I e^ Zlkx h*(9,cS)d9 
Jo Jo 

+^J o e ^ ^ky(t, Z(0), Z*(0)) J o exp[ J / 2 (x, 77, «)rfr]/ 3 (x, 7/, a)rf«*0 
+— J e J exp[ J / 2 (*, ?7, a)dr\f 3 (x, 77, a) J V/ 2 (x, 77, f , yfdrdsdO 

+ 27r 1 I 6XP [ J /2(X ' ^ a)jT l ' V/3(X ' a) ®' ^ V )Tdsde > ^ B - 21) 

where 



V/,(x, 77, £ 7) = — (x, 77, a)y + — (x, 77, ar)£ + — (x, 77, a)y. 

ox a?7 oa 

Integral equations for f and fj are derived in a similar manner. Then, they define a sys- 
tem of linear integral equations of y, Y k , fj, and initial conditions of them are linear 
in Z*(0) = (Z\(0),Z^(0), ■ • ■). This implies that solutions of these integral equations are 
linear with respect to Z*(0). Therefore, we can write 

Y k (t, Z(0), Z*(0)) = dY k (t, Z(0))Z*(0), (B.22) 
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and 



I y t (r,z(0),z*(0))> = dY k (t,z(0)f\z*(0)), 



(B.23) 



Then, {dY\{t, Z(0)), dY 2 (t, Z(0)), • • • ) defines a linear operator from T into !F, and the dual 
operator (dYi(t, Z(0)) x , dY 2 (t, Z(0)) x , • • •) defines a linear operator from i{T) into j'CF). 
For the integral equations of y, Yk, £, fj above, we again use the iteration method to prove 
that Yk(t, Z(0), Z*(0)) is continuous in a parameter Z(0) and an initial condition Z*(0) with 
respect to the weak dual topology uniformly in k (the detail is the same as the proof of 
Prop.B.3). Therefore, the mapping (dY^t, -) x ,dY 2 (t, ■)*,■■ ■) ■ KT) x i{T) -* i{T) is 
continuous with respect to the metric d^, which means that solutions of Eq. fl7.30l) are C 1 
in initial conditions and (dY l (t,Z(0)) x ,dY 2 (t,Z(0)) x , • • •) gives the Frechet derivative at 

\m) = (\zm),\z 2 io)),---). ■ 
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